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Bromine bombarded with 5-6-Mev protons emits electrons with half-lives 132 seconds and 


55+2 seconds. A 8-spectrograph shows these to be conversion electrons from at least two 
gamma-rays of energies 127 and 187 kev. Physical and chemical tests indicate the activity is 
from a krypton isotope. The long period Kr7**! reported by Snell is found to have a half-life of 
34.5+1 hours and a positron upper limit of about 0.4 Mev. Iodine targets similarly emit 
electrons with a half-life of 75-1 seconds, shown to be conversion electrons from at least two 
gamma-rays of energies 125 and 175 kev, coming from a xenon isotope. Radiation consisting of 
electrons, x-rays, and probably gamma-rays, decaying with a half-life of 342 days is also 
found. Secondary electrons produced in the target material or in Pb absorbers by x-rays from 


the source account for about 3 of the total activity measured by the electroscope. 


BROMINE 


EAD bromide, when bombarded with pro- 
tons of 5-6 Mev energy was found to emit 
electrons in two radioactive periods of half-lives 
1342 seconds and 55+2 seconds (Fig. 1). The 
same activities were found in potassium bromide. 
Since no short periods were found in potassium 
or lead, these activities must be due to some 
process taking place in the bromine. 

To test whether the element responsible was 
Br (¢=35) or Kr (= 36), the potassium bromide 
target was dissolved and boiled to dryness, which 
procedure caused the activity to disappear. A 
proof that bromine is not volatile under similar 
conditions was given by bombarding potassium 
bromide with slow neutrons, forming Br®**®, 
and then boiling the target to dryness. Since the 
activity remained in this case, we believe the 
55-second period to be due to a krypton isotope. 


* Member, the Institute for Advanced Study, Princeton, 
New Jersey. 
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Because of its shortness, no definite information 
could be obtained on the 13-second period, 
although the indication was that it also is 
volatile. 

A further check for the source of the activity 
was made with the apparatus shown in Fig. 2. 
Lead bromide was coated on the nickel foil, and 
bombarded in vacuum. Solid carbon dioxide and 
alcohol slush was then put into the small trap, 
and liquid air in the larger one, while krypton 
gas flowed into the apparatus. The foil was then 
heated to vaporize the lead bromide and release 
any occluded gas. Free bromine should be frozen 
out by the carbon dioxide, and some of the 
krypton by the liquid air. The activity on the 
liquid-air trap was then measured and found to 
be of about 55-seconds half-life. As a check 
experiment the carbon dioxide trap was left 
empty and the CO,-alcohol mixture put in the 
liquid-air trap. A bombardment and heating of a 
new lead bromide target then caused less activity 


CREUTZ, DELSASSO, SUTTON, WHITE, BARKAS 


1342 SEC 


55+ 2 SEC. 


ELECTROSCOPE DiviSIONS / SECOND 


a 


TOTAL ACTIVITY 
FROM BROMINE 


Fic. 1. Short Periods from 


tons. 


100 200 300 400 500 #4600 


to appear on the liquid-air trap than in the 
previous experiment, with little or no 55-second 
period measurable. 

These experiments make it seem probable that 
the electrons are from krypton. The stable 
isotopes of bromine and krypton are shown in 
Fig. 3. Krypton could be produced in the 
following ways: 


(1) Br78(p, 2) 
(2) Br7%8(p, -y)Kr8® 8, 


Because the occurrence of negative electrons 
suggested an internally converted gamma-ray, a 
B-ray spectrograph was fitted with a vacuum lock 
for introducing targets repeatedly until sufficient 
effect was built up on the photographic film. 
This lock is described in a note soon to appear in 
Review of Scientific Instruments. One thousand 
and one bombardments were made to obtain the 
spectrum. Bombardment time was about one 
minute, while fifteen seconds elapsed between 
removing the target from the hydrogen-filled 
bombarding chamber of the cyclotron and 


SEC. 


placing it in the spectrograph. Since the intengi- 
ties of the two periods were about the same 
immediately after bombardment, the longer 
should emit approximately four times as many 
electrons as the shorter. About one-half of the 
short period electrons were lost while the target 
was being placed in the spectrograph, while only 
one-half of the long period electrons were used 
during the one-minute exposure in the spectro- 
graph. Therefore the lines due to the 55-second 
period should be about four times as strong as 


Fic. 2. Apparatus used to freeze out the 
radioactive krypton. 
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those due to the 13-second period. For this 
reason the stronger two lines may be tentatively 
assigned to the long period, and the weaker 
lines to the 13-second period. In Fig. 3, which 
reproduces the spectrogram, three lines are 
visible. The very weak fourth line can be seen 
but poorly on the original. The spectrograph was 
calibrated with the active deposit of thorium, 
and on various trials the magnetic field (furnished 
by an electromagnet and storage batteries) could 
be reproduced to within 3 percent. 

The energies of the electron lines are 0.113, 
0.125, 0.173, and 0.185 Mev, respectively. The 
separation of the K and Zin x-ray absorption 
edges in bromine is 0.0116 Mev, and that in 
krypton is 0.0124 Mev, while the separation of 
the two lower energy electron lines is observed 
to be 0.012, and of the higher two, 0.012 Mev. 
Because the uncertainty in measuring the lines 
is about 1 percent, this cannot be considered 
definite evidence that the lines are from krypton, 
but if so, the energies of the gamma-rays are: 


y (strong) =0.125+-0.002 (=Kr Zin absorption 
edge) 
=0.127 
=0.113+0.014 (=Kr K absorption 
edge) 
=0.127 Mev. 
y (weak) = 0.185+0.002 = 0.187 
=0.173+0.014=0.187 Mev. 


RADIOACTIVITY OF BROMINE AND IODINE 


TABLE I. Stable isotopes. 
Bromine and krypton 


Sa 
Xe 
MEV 


Fic. 3. Upper photograph: conversion electrons from 
krypton. Lower photograph: conversion electrons from 
xenon, 


The uncertainty in these values is about 1.5 
percent. If the gamma-rays are from bromine, 
their energies will be lowered by 1 percent, since 
the K and L electron binding is smaller than in 
krypton. 

Another check on the element was made’ by 
studying the x-rays accompanying the 55-second 
period. No absorbers are available to distinguish 
Br characteristic x-rays from those of Kr by 
use of the K-absorption edges, so it was necessary 
to use the L edges of Au and Pt. A thin layer of 
lead bromide was melted onto 0.2-mil platinum, 
and bombarded for two minutes. The 13-second 
period was allowed to decay for 45 seconds, 
after which 3 to 4 electroscope readings were 


83 


84 


11.53 


57.10 


Iodine and xenon 


129 


130 


0.094 0.088 1.90 26.23 


4.07 21.17 | 26.96 10.54 $8.95 


483 
Kr 
from 
rensi- 
same 
onger 
Nany 
f the 
| 
used Poy | 78 79 80 81 82 = = 85 86 
cond 
as | | 50.6 49.4 
| 0.35 | 2.01 11.53 | — 17.47 
| 
| 

- ~ 124 | 125 | 126 | 127 | 128 | mm | mm | 131 | 132 | 133 | 134 | 135 | 136 

sal | 100 


CREUTZ, 


DELSASSO, SUTTON, WHITE, BARKAS 


Kr7®8 


2 
Pe 34.5 + 1 HOURS 
— 
> 
r=) 
tol Fic. 4. Positron decay of 
a 
° 
w 
02 
01 l lL l l 
° 24 48 72 120 144 
HOURS 
3.0 04 
2.0 
TOTAL Ac Se - 
ACTIVITY ABSORPTION OF 
75+ 1SEC FROM 02+ 75-SECOND 
1ODINE > RADIATION 
FROM IODINE 
sb 
2 
.005 ° 
tre l 
° 1000 2000 MG 
8 cme 
5 Fic. 6. Absorption of 75-second radiation from 
we showing the x-rays are characteristic of xenon. Total 
wu activity was about 1 division/second. Absorption coeffi- 
02 | ison, X-Rays in ory and Experiment, are as follows: 
° 100 ©6200 «6300 400 500 600 SEC. For iodine characteristic K radiation: in Sn, «/p=9; Ag 


Fic. 5. Short period from iodine bombarded by protons. 


made with the absorber in place. The absorber 
was then removed, and 6 to 8 readings taken. 
The data suggest that thex-raysarefrom krypton. 

Kr?! postulated in reaction (1) above is not 
stable. It has been produced by Snell! by deuteron 
bombardment with a half-life of the order of 
magnitude of 18 hours. We have bombarded lead 
bromide and also sodium bromide with protons 
for about 30 hours and obtained a half-life of 
34.5+1 hours which magnetic analysis shows to 
consist of positrons. See Fig. 4. Examination of 
these positrons in a cloud chamber reveals an 
upper energy limit of 0.4 Mev. Scattering in the 


1A. H. Snell, Phys. Rev. 52, 1007 (1937). 
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gas is rather bad at this energy, however, so the 
uncertainty in this value is about 0.1 Mev. 
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IODINE 


The procedure with iodine closely parallels 
that with bromine. The stable isotopes, and 
those of xenon, are shown in Table I. When lead 
jodide, potassium iodide, and lead on which 
iodine had been melted and recrystallized were 
bombarded with protons, an activity of 75+1 
seconds was observed (Fig. 5). This activity was 
volatile, while I"** formed by slow neutron cap- 
ture was not volatile under similar conditions. 

Two §-ray spectrograms were made of the 
electrons, one using 801 lead iodide targets, and 
one using 1301 potassium iodide targets. Figure 3 
shows that obtained from lead iodide which 
differs in no important respect from the other. 
Three lines are visible, at 0.0914, 0.140, and 0.170 
Mev, respectively. The K—Lin separation in 
xenon is 0.029 Mev, while that in iodine is 
0.0278 Mev. Here, as in the case of bromine, the 
error is too great to assign definitely the observed 
separation (0.170—0.140=0.030) of the upper 
two lines to either element. The lowest energy 
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Fic. 8. Absorption of 34-day period from I7+H}, 
Curves A and B show absorption in aluminum. Curve C 
was made with lead absorbers but with 31 mg/cm? of 
aluminum between the absorber and the electroscope to 
remove the soft secondary radiation produced in the lead. 
Total activity was about 2 divisions/minute. 
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Fic. 9. Absorption of 34-day period from I"*’+H". The 
filled-in circles show the absorption by lead alone; the 
open circles, by aluminum alone. The points marked x 
were obtained by putting various thicknesses of aluminum 
between 318 ane? of lead and the electroscope. The 
points marked + show the activity observed through lead 
absorbers, but with 31 mg/cm? of aluminum absorbing the 
secondary electrons produced in the lead. 


line appears to be single, although one several 
times weaker could be lost in the background. 
If this is assumed to be a K line from xenon, the 
gamma-ray energy is: 


0.0915+0.034 (Xe K-absorption edge) 
=0.125 Mev. 


The higher energy gamma-ray is, similarly, 
0.140+-0.034=0.174 


or 0.170+0.005 (Xe Lis:-absorption edge) 
=0.175 Mev. 


The uncertainty in these values is about 1.5 
percent. A chemical investigation showed that 
the amount of bromine in the lead iodide, or 
iodine in the lead bromide, if any at all, was 
much too small to account for the coincidence in 
energy of the two gamma-rays (Br, 127 kev and 
I, 125 kev). 

An attempt was made to freeze out the active 
material with liquid air, but the results were 
less positive than in the case of bromine, which 
is perhaps not surprising, as the vapor pressure 
of xenon at liquid-air temperature is about 13 
centimeters, and only a small amount was used 
as a carrier. However, the volatility of this 
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activity makes it seem that a xenon isotope is 
responsible. Furthermore, the x-rays accompany- 
ing the internal conversion were shown to be 
characteristic of xenon. See Fig. 6. A more com- 
plete study of the absorption in silver showed the 
ionization due to the x-rays was about equal to 
that of the unconverted gamma-rays. 
The possible ways of producing xenon are the 

following : 

(3) I*"(p, m)Xe¥" 

(4) 


Xe"8 is stable, and the gamma-rays may arise 
from an excited state in it. Xe? is not found in 
nature, nor has it been reported as a radioactive 
isotope. We have bombarded lead iodide and 
also sodium iodide for about 30 hours and ob- 
tained an activity of 34+2 days half-life (see 
Fig. 7), consisting of electrons, x-rays and 
probably some gamma-rays. The range of the 
electrons is about 350 mg/cm? of aluminum. 
By means of a counter kindly supplied by Mr. 
J. G. Fox, they were found to be negative. 
A softer (negative) group also appears in the 
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complete absorption curve (Fig. 8), with range 
about 35 mg/cm? of aluminum. These, however 

are undoubtedly secondaries produced by the 
harder radiation in the material of the source 
and the absorber, since the same group is found 
by absorbing in aluminum the radiation that 
comes through various thicknesses of lead. The 
greater number of these secondaries produced jn 
lead than in aluminum makes the total radiation 
received through a given mass of lead greater 
than that through the same mass of aluminum, 
until the total absorber thickness is about 0,9 
gram/cm?. This effect is shown in Fig. 9. 

A study of the x-radiation with critical ab. 
sorbers in an attempt to ascertain whether the 
x-ray is from iodine or xenon gave an incon- 
clusive result. From the absorption coefficient jn 
Pb, its energy is found to be about 40-45 key, 

Although the number of nuclei with the 75. 
second half-life made in a given bombardment is 
of the same order of magnitude as the number 
decaying with the 34-day period, no genetic 
relationship between the two processes has as 
yet been established. 
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Internal Conversion of Gamma-Radiation in the Z Shell 


M. H. Hess 
Department of Physics, Duke University, Durham, North Carolina, 


AND 


ELpRED NELSON, 
Department of Physics, University of California, Berkeley, California 
(Received May 31, 1940) 


The nonrelativistic calculations of the internal conversion of electric multipole radiation 
which have been made for the K shell are extended to the L shell. Numerical values of the 
conversion coefficient for both the K and L shells are given and curves showing the ratio of the K 
conversion to the LZ conversion for different energies and atomic number are presented. This 
latter ratio is quite sensitive to the multipole order and varies between ~0.1 and 10 in the range 
in which the calculations are valid. A simple relativistic formula for the conversion of magnetic 
multipole radiation is given. In Section 5 we summarize selection rules and give applications of 


the formulas to experiments. 


1. INTRODUCTION 


HE increasing number of cases in the lighter elements in which the internal conversion of 
gamma-radiation has been observed makes a study of this phenomenon more important. Cal- 
culations of the ratio of the number of conversion electrons to the number of gamma-quanta N,/N, 
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CONVERSION OF GAMMA-RADIATION 


have already been carried through for electric and magnetic multipole radiation for electrons of the 
K shell.! In this paper we shall extend these calculations with the same restrictions on their validity 
to electrons of the L shell. This has seemed worth while because experimentally it is much easier to 
obtain the ratio of the conversion in the K shell to that in the Z shell Nx/Nz than to measure either 
Nx/N,or N ./N, separately. Consequently, with calculations of both Nx/N,and N,/N, available, 
an easier and more accurate comparison of experiment and theory can be made. In order to facilitate 
this comparison we have made fairly extensive numerical calculations of both quantities Nx/Nq 
and N,/N, and of the ratio Nx/N1. 

First we present the calculations for conversion of electric multipole radiation; in Section 4 we 
consider the conversion of magnetic multipole radiation. 


2. CONVERSION OF ELECTRIC MULTIPOLE RADIATION IN THE L SHELL 


In many of the artificially radioactive gamma-emitters nonrelativistic calculations of the con- 
version of electric multipole radiation are appropriate ; namely, where the energy of the gamma-ray 
is small so that for the ejected electron vc and where Z is less than ~50. We present such calcu- 
lations following rather closely the method of Dancoff and Morrison. Therefore we shall give only 
the most important steps and the results, referring to their paper for details. 

We are concerned with the ratio N./N, where N, is the number of electrons per second ejected 
from a given state or shell by the gamma-radiation from the nucleus and N, is the number of gamma- 
quanta per second which escape from the atom. In the nonrelativistic approximation in which the 
wave-length of the gamma-ray is large compared with atomic dimensions, we have? 


(1) 


=2n — 


Ny 1+1L 24! 


N. L f(2l) 


f 


for a single electron initially in the state Yo. For simplicity in writing we use the natural system of 
units in which the units of mass, action and velocity are, respectively, m, h and c. Thus in Eq. (1), 
W is the energy of the gamma-ray in units mc? and 27/W is the wave-length of the gamma-ray in 
units of the Compton wave-length. Equation (1) holds for an electric 2'-pole, 1 being the change of 
angular momentum of the nucleus in the gamma-transition. We use a to represent the fine-structure 
constant ~1/137, Y," for the normalized surface harmonic and y,; for the final electronic state in 
the continuum, normalized to unit flux at infinity. We have 


+1+in) | 
(21’+1)! 


+1 + in; +2; 2ipr), 


where p is the momentum of the ejected electron (in units mc) and n= y/p=aZ/p. 
For the L shell we have the 2s-function 


Vo= (y*/8x) —yr/2)e-7"? 


and the three 2p-functions 


The angular integrations in Eq. (1) require /’=1, m’=m if yo is the 2s-function and /’ =11, m’=m, 
m+1 for the 2p-functions. The radial integrations can be performed in series and we find for two 


'M. H. Hebb and G. E. Uhlenbeck, Physica 7, 605 (1938); S. M. Dancoff and P. Morrison, Phys. Rev. 55, 122 (1939). 

* Equation (1) can be obtained most easily by taking for the perturbation acting on the electron the potential of a 
static multipole V=const. r~'~! Y," multiplied by e~*™*. It is also in agreement with the method of Dancoff and Morrison 
as can be shown by transforming the potentials from their primed gauge for which div A’ = ¢’=0 to the conventional 
= In the conventional gauge the only important potential in the nonrelativistic approximation is the scalar po- 
ial ¢. 
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1 Y 
mane” 
where p=(y/2)+ip and F,'=F(q,l1+1+in; 214+2;4/(2—in)). 

Similarly for the six 2p-electrons, we find two terms. The first of these for which /’=/—1 is 
N. of 

N, 2 
and the second where /’=/1+1 is 


The sum of the expressions (2), (3) and (4) then gives the quantity N1/N, for the eight electrons 
of the complete L shell. 
To evaluate the hypergeometric functions we apply the Gauss relations between contiguous hyper. 
geometric functions* and the formula for the analytic continuation of the hypergeometric series. 


Then 


1 2 
J 


ial 291)” 


2041 (44+n2)2 | 
TABLE I. Values of logis [y?'#(N./N,)] for electric 2'-pole radiation. W=energy of gamma-ray in units mc, y=Z/131, 


N1 
=128r 
Ng 


A? » (3) 


K SHELL L SHELL ! 

0.125 1.67 2.11 3.83 4.95 5.71 
0.25 3.80 0.69 1.98 3.02 3.89 
0.5 3.73 3.67 3.20 3.47 2.58 3.94 J.07 0.14 0.96 1.65 

3.25 3.33 3.27 3.17 3.06 3.42 3.22 I. 1.72 0.28 
1.0 3.90 3.97 3.93 3.86 3.79 3.05 3.62 2.30 3.84 1.29 

3.62 3.65 3.59 3.50 3.41 4.76 3.18 3.72 3.16 35) 
1.5 3.39 3.38 3.28 3.16 3.04 4.52 4.82 3.25 3.60 3.89 

3.19 3.14 3.01 4.85 4.70 4.32 4.52 4.85 3.13 3.35 
2.0 3.02 4.93 4.76 4.57 4.38 4.14 4.27 4.51 4.72 188 

4.87 4.75 4.54 4.32 4.09 5.98 4.05 4,21 4.37 44] 
2.5 4.73 4.57 4.34 4.09 5.83 5.84 5.85 5.95 4.05 iil 

4.60 4.42 4.15 5.87 5.58 5.71 5.67 5.71 5.76 B78 
3.0 4.49 4.27 5.98 5.67 5.36 5.60 5.51 5.50 5.50 9-547 

4.38 4.14 5.82 5.48 5.14 5.49 5.36 5.30 5.26 5.0 
3.5 4.28 4.02 5.67 5.31 6.95 5.39 5.23 5.12 5.04 3-894 
4.0 4.10 5.79 5.40 6.99 6.58 5.21 6.98 6.64 648 
4.5 5.94 5.59 5.16 6.71 6.26 5.05 6.77 6.52 6.30 =: B08 
5.0 5.80 5.41 5.94 6.45 7.96 5.90 6.58 6.27 7.99 T77 
5.5 5.67 5.24 6.74 6.22 7.69 6.77 6.40 6.04 7.71 7.39 
6.0 5.55 5.09 6.56 5.00 7.44 6.65 6.24 7.84 746 7.09 
6.5 5.44 6.95 6.39 7.80 7.21 6.54 6.10 7.65 7.23° 882 
7.0 5.34 §.82 6.23 7.62 7.00 6.44 7.96 7.48 7.02 8.57 
7.5 5.24 6.70 6.08 7.44 5.80 6.34 7.84 7.32 5.82 834 
8.0 5.16 6.59 7.94 7.28 8.61 6.25 7.72 7.17 5.64 8.2 
8.5 5.07 5.48 7.81 7.13 8.43 6.17 7.61 7.03 5.47 4.92 
9.0 6.99 6.38 7.69 8.98 8.27 6.09 7.51 8.90 8.31 9.73 
9.5 6.92 6.28 7.57 8.84 8.11 6.01 7Al 8.77 8.15 9.55 
10.0 6.84 6.19 7.46 8.71 9.96 7.94 7.31 8.65 5.01 9.38 


3 Gauss, Werke 3, 130. 


4 Whittaker and Watson, Modern Analysis(Cambridge, 4th ed.), Sec. 14.51. 
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Fic. 1. Curves for Nx/Nz as a function of Z?/F. 


A, =2(1+1)[44+ tan — B,, 
Ast = + 152?+311+16) tan Bt, 
The B’s are conveniently specified by giving the recurrence formulas which they satisfy. 


[4(3/2-+91+4) 15/?+ 311+ 16)? = [4(3/2+15/+16) 
| 14242 
(2/+1)! (4-+n?)! 
(l—2) i)?+n?]---[1 +n*] 
(21—1)! 


= (1+-2)B,-— 
The first few values are 
B,=0, B,*=0, 
B.=j, =§, 
Bs-=}. 
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3. NUMERICAL RESULTS 


In Table I are shown the results of the calculations. The values for the K shell given in the left 
half of the table were calculated from the formulas previously obtained by Hebb and Uhlenbeck 
and by Dancoff and Morrison. The right half of the table refers to the L shell and was compute 
with the aid of Eq. (5). We have tabulated the quantity logio [y?'**(N./N) ] since this depends on 
W and Z only in the combination W/y*. The energy of the gamma-ray is related to the momentum 


of the electron and the parameter n by 


W/y?=(4+n*)/8n° 


for the K shell and by 


for the L shell. If EZ represents the energy of the gamma-ray in kev, then W=E/511 and W/,! 
= 36.7E/Z*. 

The influence of screening is easily taken into account in Eq. (5) and the table. One has simply 
to replace Z by Z.1s:=Z—o where o is a screening constant. According to Slater® one should take 
oa =0.30 for the K shell and =4.15 for the L shell. With neglect of screening it is possible to express 
Nx/WN_ in terms of one parameter W/7? or n, but when screening is taken into account this is no 
longer true, and Nx/N,z depends essentially on both W/7? and Z. However, the dependence on Z 
is not strong, and it is still feasible to construct approximate curves giving Nx/Nz directly. These 
are shown in Fig. 1 plotted against Z?/E. They are strictly correct only for Z=35; for higher Z 
the curves should be lowered somewhat and for smaller Z raised compared to the positions shown, 


With 25<Z<50 the departures are ~ 10 to 20 percent. 


4. CONVERSION OF MAGNETIC MULTIPOLE RADIATION IN THE L SHELL 


Since the conversion of magnetic multipole radiation by s electrons depends essentially on the 
electron spin,! its calculation must be relativistic. In the lighter elements for gamma-ray energies 
not too near the threshold, the binding of the electron will not play an important part and a good 
estimate of the magnetic conversion may be obtained by using the Born approximation in evalu- 
ating the matrix elements. We have calculated a simple relativistic formula giving the conversion 
of magnetic multipole radiation in the LZ shell similar to that obtained by Dancoff and Morrison for 


the K shell. 
The number of electrons ejected per second from the L shell into a solid angle dQ is 


2re? 
2, 


where 
2 


A=0 2, spin 


f dev 


a is the vector whose components are the first three Dirac velocity matrices, A’ is the vector potential 
of a magnetic multipole in the gauge div A’=0, and yy is normalized to unit energy. The sums are 
taken over the spin and orbital states of the L shell. To evaluate the matrix elements we use the Born 
approximation, neglecting the influence of the Coulomb field and taking a plane wave for the space- 
dependent part of ¥;. Doing the sums and integrals and dividing by the rate of radiation, we have 


for eight electrons in the L shell: 

N, 4 2i+1 W+2 
$j. ~ 85. C. Sk Slater, Phys. Rev. 36, 57 (1930). 
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The first term in Eq. (6) comes from the 2s electrons and is just one-eighth of the conversion for the 
1s electrons of the K shell. As pointed out by Dancoff and Morrison, conversion of magnetic radiation 
by s electrons is essentially relativistic. On the other hand, this is not so for p electrons for which 
nonrelativistic calculations give a finite conversion. The second term in Eq. (6) is the contribution 
from the 2p electrons. The ratio of the magnetic K conversion to the magnetic L conversion is 
readily obtained, since the formula for the K shell' 


/ 
w\ w 


is contained as a factor in Eq. (6). Like the corresponding ratio for electric multipole radiation, it is 
sensitive to the multipole order, but for a given case it is definitely larger than that for the electric 


multipoles. 
5. APPLICATION TO EXPERIMENTS 


A gamma-ray of energy W is emitted in a nuclear transition between states of total angular 
momentum J and J’ and energy difference W. From two constants of the motion, the parity and the 
total angular momentum of the system nucleus+gamma-ray, we obtain the selection rules.! The 

ma-ray will be predominantly multipole radiation of the lower order / in the range | J—J’| =/ 
<|J+J’| permitted by the parity and not reduced by special arguments of symmetry. If it is 
permitted by the parity, we need consider only electric multipole radiation of order 1= | J—J'| ; 
if it is not permitted by the parity, the gamma-ray will, in general, be a mixture of magnetic 2!-pole 
and electric 2'*!-pole radiation. It would seem that the occurrence of these two cases should be 
equally probable. 

The type of radiation involved may be classified by fitting the experimentally determined values 
of the total internal conversion coefficient, the K to L conversion ratio, and the lifetime, each of 
which is sensitive to /. Since the dependence on / of the K to L ratio for electric and magnetic multi- 
poles is different, it provides a method of distinguishing between them. A mixture of the two will 
appear in a larger K to L ratio than is consistent with the lifetime for electric multipole radiation. 

In two cases sufficient information is available to assign the conversion to electric multipole 
radiation. In Zn®’ the observed total internal conversion coefficient of ~0.5 to 1.0,* the K to L ratio 
of 8,’ and the short lifetime are in good agreement with the calculated values for electric quadrupole 
radiation. The 6.4-hour Cd activity has a K to L ratio of one. This could be fitted to electric multipole 
radiation for ]=4 if the conversion electrons came from a short-lived product of the disintegration. 
Experiments have shown that the gamma-ray and the conversion electrons actually come from a 
40-sec. isomer of Ag.* The 170-kev internally converted gamma-ray of the 1.2-day Te has a K to L 
ratio of 2.* This can be interpreted consistently if the gamma-ray is } magnetic multipole radiation, 
l=4, and 3 electric multipole radiation, /=5. 

From the preceding applications it is clear that the ratio of the K conversion to the L conversion 
Nx/N, will be useful in determining the multipole order of gamma-rays of energy «mce?. For gamma- 
rays of energy mc? the K to L ratio for both electric and magnetic multipoles approaches the 
same limit ~10 for all multipole orders. 

The authors wish to thank Professor J. R. Oppenheimer, Professor G. E. Uhlenbeck, and Dr. 
P. Morrison for their interest and suggestions. 


* The authors are indebted to Dr. A. C. Helmholz for this information. 
7G. E. Valley and R. L. McCreary, Phys. Rev. 56, 863 (1939). 
*L. W. Alvarez, A. C. Helmholz and E. Nelson, Phys. Rev. 57, 660 (1940). 
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Theoretical Maxima and Minima in the Anomalous Scattering of Alpha-Particles 


PauL SELIGMANN 
Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois 
; (Received June 8, 1940) 


The work of Wenzel on the determination of theoretical maxima for the anomalous scattering 
of a-particles by light nuclei under the assumptions of the two-body problem has been general- 
ized to the case of an finite number of / values and extended to include minima. The results have 
a more general validity if the scattering nucleus has zero spin. Possible applications to analyses 
of experimental data are discussed. Treatment of the existing data for a-particles scattered by 
carbon and oxygen indicates that / values at least as large as two are effective in producing the 


anomalies in both cases. This conclusion is verified for carbon by a phase-shift analysis. 


I. INTRODUCTION 


XPERIMENTS on the anomalous scatter- 
ing of alpha-particles by light elements 
heavier than He‘ have been performed by 
Rutherford and Chadwick,' Riezler,? Brubaker,’ 
and Devons.‘ Although no complete phase- 
shift analyses have been made, the data of 
Riezler have been studied by Wenzel,> who 
calculated maximum values for the anomaly in 
the scattering due to a single value of the orbital 
angular momentum of the incident particle. 

In the present paper, the work of Wenzel is 
extended to include the derivation of expres- 
sions for maximum anomalies in the scattering 
due to combinations of two or more values of the 
orbital angular momentum of the incident 


R=|1+ 
n 


where 6 is the angle of scattering in the center of 
gravity system, »=2Ze?/hv with v the velocity 
of the incident a-particles, o,x=arg T'(/+1+7n), 
and the P; are zonal harmonics. The K; are 
here treated as arbitrary parameters charac- 
teristic of the nuclear interaction, dependent 
on the energy but independent of the scattering 
angle. As is usual in empirical phase-shift 
analyses, the sum in the expression will be taken 
to consist of only a few terms. 


1E. Rutherford and J. Chadwick, Phil. Mag. 50, 889 


(1925). 

2W. Riezler, Proc. Roy. Soc. 134, 154 (1931); Ann. d. 
Physik 23, 198 (1935). 

3G. Brubaker, Phys. Rev. 54, 1011 (1938). 

4S. Devons, Proc. Roy. Soc. 172, 127 (1939). 

5 P, Wenzel, Zeits. f. Physik 90, 754 (1934). 
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particle as well as expressions for the minimum 
anomalies under the same conditions. The ex. 
pressions are then applied to the computation of 
maxima and minima for the cases of carbon and 
oxygen. Comparisons with the experimental 
data are carried through. It is found that values 
of / at least as large as two are needed to account 
for the observed results. This is verified in the 
case of carbon by an actual phase-shift analysis 
on the Riezler data. 


II. THEORETICAL MAXIMA AND MINIMA 


The analysis to follow is based on the assump- 
tion that the scattering anomaly R for alpha- 
particles (ratio of the scattering cross section to 
that for a pure Coulomb field) can be written 


in the form :° 


The above form for R is the familiar one which 
can be deduced rigorously when the interaction 
is treated as a two-body problem with the forces 
derivable from a potential dependent only on the 
separation of the a-particle from the nucleus 
and behaving at infinity like a Coulomb field. 
The K; are then the usual phase shifts over 
Coulomb for the radial wave functions at in- 
finity. However, there is reason to believe that 
an expression of the form:(1) has more general 
validity, with, of course, a modified interpreta- 
tion of the K;. Thus, according to Bethe,® the 
following assumptions are at least necessary and 
possibly sufficient for its validity : 

*H. A. Bethe, Rev. Mod. Phys. 9, 69 (1937), §73. 
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(a) There is no competing reaction, i.e., there is no other reaction such as disintegration of the 
nucleus which can occur with a probability comparable with that to be described. 
(b) Orbital angular momentum is conserved. This condition is most easily satisfied in cases 
where the scatterirz nucleus has zero spin, but might be satisfied in other cases as well if the 
bility of transformation of orbital angular momentum of the a-particle into spin angular 
momentum of the nucleus is negligible. 
In order to bring out the dependence on the K:, Eq. (1) may be rewritten in the form 


R= S, exp [ir S, exp [i(a:+2K,)]|*=|A exp S,exp [i(a:+2K),) ]|?, 


where (214+-1) sin? 36 


| Pi(cos 6)|, 
when P,(cos 6) >0 


=» log sin? 
log sin? 36+2(o when P,(cos 6) <0 


and exp [ic;] 


with A and 6 real. 
An equivalent and more convenient form for R is 


R=|A4+2 S, exp (4) 


The maxima of R as a function of any tates number of K; at a fixed value of @ are now er in- 
vestigated, since Eq. (4) expresses R as the square of the length of the sum of a number of vectors in 
the complex plane, and the maximum will be attained when these vectors are all colinear. The de- 
sired condition for a maximum is therefore simply that 


a:+2K,—6=0 (5) 


for each value of / involved in the analysis. 

Since, for the analysis of experimental data, the K; are considered arbitrary parameters, a set of 
K, values can always be found which will satisfy the set of conditions represented by (5). The maxi- 
mum of (4) is then 


R™===|A +2 Si|?=A?+2A Std SP+2 SiSv. (6) 
The quantity A, evaluated from (3) is 

A= (+2 S?—2 cos a,+2 Si:Sv cos (a:—ay))*. (7) 
By the substitution of (7) into (6), and further reduction, one obtains 


1+D+2(2 Si)(1+D— Si)?)4, (8a) 
where 


Equations (8) then represent the maximum: possible scattering for any case for which Eq. (1) 
is valid. 

The minimum scattering may be calculated in an analogous fashion. Examination of Eq. (4) shows 
that the minimum will be zero unless one of the vectors is greater than the sum of all the others. 
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When a minimum other than zero is possible it is convenient to consider two cases separately: 


Case i. S;. Then 


S,|? 

Si (9a 
and substitution of A from (7) yields 
S)(1+D— Si)*)'. (10) 
Case ii. S;>A+ Sp. Then 
Vel 
Sy |?*. 


In this case it is found more convenient to compute the minima from the magnitudes of the vectors 


concerned rather than to substitute (7) into (9b). 
A procedure analogous to the above may be readily carried out for symmetric cases of scattering 
such as the scattering of a-particles by helium. The expression corresponding to (1) is? 


in log cos? 6 


sin? @ cos? 6 leven 


2 
1 1 ( log tan 6) 


e~ in log sin? @ 


sin‘ @ cos‘ sin? 6 cos? 6 


which may be reduced to 
| 1+tan? be” log tan? @__ T T | 2 


leven leven 


1+tan‘ 6+tan? 6-7 cos (7 log tan? 6) 


where 


2(21+-1) sin? @ 
| P.(cos 26) |, 


when P;(cos 26) >0 
+3n/2 when P;(cos 26) <0. 


The derivation of an expression for the maximum may be carried through as before. The final result is 


1+D,—tan? 6+2tan?@ © T)(1+D.—( T,—tan? 6)?}! 


l even even even 


1+tan‘ 0+tan? 6-7 cos log tan? 6) 


Rmx = 


where 
D,=2 tan? 6 cos log tan? 6)+2 tan‘ Ti cos 


even even 


+ TiTv cos? }(8:—By) —2 tan? @ 7; cos* log tan? 


U>I even even 


In the case of the minimum, it is necessary to find the quantity 
A,=[1+D,—( T:—tan* 6)*}! 


even 


and compare it as before with the values of 7:. 
7H. A. Bethe, reference 6, §74. 
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Expressions for maxima and minima may be 
applied to experimental data in various ways, 
the most obvious of which is the comparison of 
theoretical maxima with those experimental 
data which give the variation of the scattering 
with energy at constant angle. This was the 
type of application made by Wenzel. The pro- 
cedure consists in making assumptions about the 
] values effective in a specific case and then com- 
paring the maxima or minima calculated on this 
basis with the experimental data. If the com- 
parison shows that the experimental curve 
crosses the theoretical extremum curve, then 
the particular set of / values chosen is ruled out. 
As will be shown below, if it so happens that 
the experimental energy variation curve is 
tangent to a maximum or minimum curve, then 
the point of tangency will determine uniquely 
a set of phase shifts. The validity of such a set 
can be tested if the experimental angular distri- 
bution curve for that particular energy is known, 
since a theoretical angular distribution curve 
can be calculated from the phase shifts and then 
compared with the data. Such cases of tangency, 
although looked for, were not found in the ap- 
plications to be discussed below. 

It is also of interest to examine the behavior 
of the extremum curves as a function of angle. 
For the case of a single phase shift and a given 
constant energy, the maximum curve plotted 
as a function of angle will be an envelope of all 
possible angular distribution curves for varying 
values of the single K;. From (5), the value 
of the phase shift is given by 


Ki=}(6—a)). (11) 
For the maximum R, the quantity 6 as deter- 
mined from (3) is 
sin ai 


tan 6 (12) 
Si cos a,;—1 


When (12) is substituted into (11), one obtains 
sin ai 


tan 2K, (13) 
Si—cos a) 


This expression establishes a correspondence 
between the K; values and the points on the 
maximum curve. For every value of @ there can 


be calculated a value of a: and also one of S; 
and hence of Kz, since the K; is usually restricted 
to the range 0 to z or —2/2 to +2/2. The same 
conclusions are readily seen to hold for minimum 
curves, the slight modification of (5) for this 
case resulting in corresponding changes in 
(13) and (14). 

These properties may be of use in simplifying 
phase-shift analyses of experimental data. If an 
experimental angular distribution curve crosses 
an extremum curve, the corresponding set of K, 
values is automatically excluded. Furthermore, 
for a single phase shift, if an angular distribu- 
tion curve is tangent to an extremum curve a 
value of the phase shift is uniquely determined 
by Eq. (15) and its ability to account for the 
observations can be tested by calculation of the 
theoretical angular distribution curve. Similarly, 
the case of two phase shifts can be handled by 
the determination of any points of tangency 
between the observed curve and any member of 
the one-parameter family of extremum curves 
corresponding to the two values of / under 
consideration. 


III. APPLICATIONS 


The cases treated will be those most nearly 
satisfying the conditions under which Eq. (1) 
is presumably valid. Of the light elements heavier 
than He‘, only carbon and oxygen satisfy the 
conditions of zero spin and nondisintegrability 
by alpha-particles. 


A. Oxygen 

The best available data seem to be those of 
Brubaker* and Devons.‘ These consist of energy- 
dependence curves at six angles. The required 
maxima and minima for all combinations of / 
values less than four were calculated and are 
shown in Fig. 1. Examination of these curves 
shows that there is no set of / values which will 
yield consistent maximum and minimum curves 
for the data, because of the rise above all the 
theoretical maxima in the experimental curves 
of Brubaker for higher energies at the largest 
angles considered. It is difficult to account for 
this phenomenon. A plausible explanation can be 
given, however, without assuming values of / 
greater than three, or that Eq. (1) is not valid. 
The theoretical curves may vary rapidly with 


(9a) 

(10) ii 
(9b) 
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@=53° 
6=64 


© = 65 30! 
6=78°39' 
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2 


ENERGY IN MEV 


mental curve at 6=90° is 


Fic. 1. Scattering of a-particles by oxygen. The heavy curves are the experimental values; the lighter ones are the 
theoretical maxima and minima, the numerals at the side indicating the / values for which they were computed. © denotes 
the scattering angle in the laboratory system of coordinates; 0, the corresponding angle in the center of gravity system. 
at of Devons; the rest are those of Brubaker. The minimum curves for /=0,2 
and /=0,3 at @=90° have not been plotted because they lie too close to zero to be observable on the vertical 


angle® and the angular spread in Brubaker’s 
apparatus was large. It is quite possible that a 
rapid rise in the scattering just above the angles 
in question would make the type of average 


§ G. Beck in Handbuch fiir Radiologie, Vol. 6 (1), Chapter 
IV, §25 (1933). 


@= 104° 
= 118° 


hb 
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ENERGY IN MEV 


value experimentally obtained much higher 
than the actual scattering at the mean angle 
given in the data. This conclusion is supported 
by the Devons data where better angular 
conditions existed. No such rise occurs, despite 
the fact that the mean angle given is very close 
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to one of the Brubaker curves being discussed. 
Another marked difference between the two 
sets of experimental curves is the fact that what 
appears in Brubaker’s work as an irregularity, 
as a maximum in the Devons data. 

It is quite possible that this “singularity” at 6 
Mev, which was ascribed by Brubaker to a 
resonance level, is in reality a maximum which 


has been covered up by the much larger scatter- . 


ing at slightly higher energies and angles within 
the angular range of the detecting apparatus. 
Under this assumption, the Brubaker and 


Devons data become consistent in showing that’ 


the scattering at this energy is due to a d reso- 
nance, for there is no combination of the values 
0, 1, and 3 for J, which will give maximum values 
higher than the experimental and minimum 
values lower than the experimental for all 


angles at this energy. 
B. Carbon 


The data available in this case are those of 
Riezler? and Devons.‘ They consist of energy 
dependence curves at two angles and one angular 
distribution curve. The maxima and minima for 
the energy variation are shown in Fig. 2. It will 
be observed that the curves differ from those of 
Wenzel,® the discrepancy being due on the one 
hand to the fact that Wenzel’s curves were 
calculated for 180°, and on the other to an error 
in Wenzel’s expression.® It is seen, in a manner 
similar to that for oxygen, that the resonance 
in the Riezler data is consistent with an / value 
of two rather than of one as postulated by 
Wenzel. 

This conclusion was verified by a phase-shift 
analysis based on Riezler’s angular distribution 
curve, which is at an energy slightly higher than 
the resonance maximum. It was found that 
Ky and K, together could not account for the 
data and that K» at least had to be used. The 
procedure adopted was as follows: 

Three points on the angular distribution curve 
determine, in principle, the three phases Ko, K, 


*A comparison of the first two terms of Wenzel’s 
expression (12) with the corresponding terms derived 
from (1) of this paper shows that the relationship between 
the parameters used is K:=y:—o:+4r. A comparison of 
the third term of Wenzel’s expression with the correspond- 
ing term derived from (1) shows that, in Wenzel's expres- 


sion, 2(¢9—0;) in the argument of the sine should be 
replaced 


209 
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Kz as the point of intersection in (Ko, Ki, K:) 
space of the three surfaces 


fi(Ko, Ki, K2)=R,, 
fe(Ko, Ki, K2)=Rz2, 
fs(Ko, Ki, K2)=Rs, 


where fi, f2, and f; are the functional forms 
which Eq. (1) takes for the values 6@;, 02, and 63 
corresponding to the experimental values Rj, Re, 
and R;. In the determination of the intersection 
of these three surfaces analytic methods were 
excluded because the complexity of the func- 
tional forms eliminated the possibility of an 
explicit solution ; graphical and numerical meth- 
ods were therefore resorted to. Tables were first 
constructed for values of f; and fe as functions of 
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@=90° 
6=109° 
ck. 
— | 
4, 
d 4 
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Fic. 2. yy of 4 by carbon. The notation 
is the same as in Fig. 1. The upper experimental curve is 
that of Riezler; the lower, that of Devons. On the lower 
plot, only a few minimum curves have been shown; the 
rest lie in the same region and very close to those given. 
On the upper plot, the maximum curves not shown all lie 
above the range of the graph. 
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Ko, Ki, and K, at intervals of 7/12 for the param- 
eters in the range 0 to 7. For each Ko=n7/12, 
(n=1, 2, ---, 12), the curves fi(m7/12, Ki, Ke) 
fe(mr/12, Ki, K2)=R: were then con- 
structed and intersections of curves for the same 
Ko were determined by Newton’s method. The 
points of intersection determine the curve in 
(Ko, Ki, K2) space on which the intersection of 
the three surfaces must lie, if, indeed, such an 
intersection exists at all. This point is deter- 
mined as that point for which fs=R3, and was 
found, again, by Newton’s method. The values 
used were: 6,;= 142° 26’, Ri=1.10; 02=166° 33’, 
R2=11.5; 63=152° 23’, R3=3.62. From these 
values, the procedure described yielded a set of 
values Ko=105°, K,=150°, and K2=94°. The 
curve calculated from these values gives ap- 
proximate agreement with the Riezler data. 


As far as the Devons data are concerned, jt 
will be seen that the low values for the scattering 
at low energies are consistent with any values of 
the orbital angular momentum with 7 less than 
or equal to three, and data at other angles are 
needed in order to settle the nature of this 
phenomenon. The experimental maximum at a 
higher energy has about the same property, but 
if the theoretical curve found above is prolonged 
to 6=109° 28’ (corresponding to 90° in the 
laboratory system), a very large discrepancy 
occurs. The question of whether this discrepancy 
is due to experimental error or to weaknesses jn 
the theory can only be settled by the publication 
of more accurate data. 

In conclusion, the writer wishes to express, his 
thanks to Professor F. C. Hoyt and Professor 
Carl Eckart for suggestions in this work. 
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The Distribution in Angle of Alpha-Particles from Lithium Bombarded with Protons 


Victor J. Younc, A. ELLett anp G. J. PLAIN 
State University of Iowa, Iowa City, Iowa 


HE distribution in angle of a-particles 
from the Li?7+H!—2He? reaction has been 
studied by Kirchner,! by Giarratana and Bren- 
necke? and by Neuert,’ all of whom report the 
distribution to be spherically symmetric. How- 
ever, the data reported by these observers have 
in every case been taken at bombarding energies 
of 270 kev or less, and with relatively large 
statistical errors. We find that even at 100 kev 
there is a definite but small departure from 
spherical symmetry, while at energies of 350 kev 
the intensity in directions near that of the 
bombarding beam exceeds that in a direction at 
90° to the beam by a factor of the order of 2. 


1F. Kirchner, Physik. Zeits. 34, 785 (1933). 
as eo and C. G. Brennecke, Phys. Rev. 49, 35 
+H, Neuert, Ann. d. Physik 36, 437 (1939). 


(Received June 28, 1940) 


The distribution in angle of alpha-particles from Li’+ H! is found to be of the form 1+A cos? 6. 
The variation of A with bombarding energy is characterized by a rapid increase in the range 250 
to 325 kev presumably indicating resonance toa broad excited level of Be’. 


This is not surprising as it appears from other 
considerations that the reaction proceeds by 
capture.‘ 

APPARATUS 


The apparatus is shown in Fig. 1. A lithium 
target formed by evaporation from the furnace 
upon one of the surfaces of the hexagonal target 
support is bombarded by a beam of magnetically 
resolved protons or H,* ions defined by the 
apertures. The two upper apertures are }” 
circular holes, while the lowest is reduced to 
3/64” in the plane of the paper. The target 
surface being tipped at an angle of 22.5° to the 
beam results in an approximately rectangular 
target spot about 3X3 mm. The sensitive volume 


‘M. Ostrofsky, G. Breit and D. P. Johnson, Phys. Rev. 
49, 22 (1936). 


°H. 
(1938), 


it 
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in front of any one of the eight collector plates 
can see this target spot through a 4” clear 
aperture made vacuum tight with 0.0015” 
aluminum foil. The hole over which the foil is 
supported is }’” in diameter so that particles 
emerging through the 3” defining apertures do 
not strike the foil in the more steeply sloping 
region near its edge. 

For counting, two independent linear ampli- 
fiers with associated scales-of-eight’ and mechan- 
ical counters were used. The data obtained are 
the number of particles observed within a given 
time interval at any chosen window and the 
number observed in the same interval at the 
three forward windows. These three were con- 
nected in parallel to the input grid of one ampli- 
fier and the chosen window to the other, unless 
the latter was one of the forward three, in which 
case only two were in parallel and the total 
count on the three was obtained by adding the 
scores of the two counters. In any case the final 
result is the ratio of the number of particles 
counted at a given window to the number at 


Fic. 1. Cross section of target chamber 
and defining apertures. 


*H. Lifschutz and J. L. Lawson, Rev. Sci. Inst. 9, 83 


(1938) 


ALPHA-PARTICLES FROM LITHIUM 


YIELD WITH 
ENERGY AND ANGLE | | 


H' 2 He* | 


| 
12 
10 UV 
° 
| 
fel 
100 
| | | | 


cos* @ 


T 
6 


Fic. 2. Distribution in angle of alpha-particles from 
Li’+H*~+2Het for 8 bombarding energies from 100 to 
400 kev. Relative numbers of alpha-particles per unit 
solid angle vs. square of cosine of angle between direction 
of emission and direction of bombarding beam in center of 
mass coordinate system. 


the three forward windows. The latter number 
constitutes a measure of the number of dis- 
integrations occurring in the interval in question. 
The use of such a monitoring arrangement 
eliminates obvious errors associated with meas- 
urement of bombarding current and the assump- 
tion of a yield proportional to the current. 

Bombarding energies are taken from a magnet 
current vs. bombarding energy curve for which 
the 330-kev fluorine y-ray resonance served as a 
fixed point, other points being interpolated by 
the usual procedure of bringing various m/e 
spots onto the target. This curve was checked 
using a resistance voltmeter. (450 10-megohm 
IRC resistors.) 

Thickness of targets was estimated at from 2 
to 40 kev stopping power, the thicker targets 
being used at low energies in order to secure 
increased yield. The estimate of target thickness 
was made by determining the yield from the 
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KEV 
Fic. 3. The variation with bombarding energy of A in 
the expression 1+ A cos* @ which is found to represent the 
distribution in angle of the alpha-particles from Li’ 
+H'!—2He't. 


target in question at a particular bombarding 
energy, then the yield from a thick target both 
at this and at a somewhat lower energy so 
chosen that the decrease in thick target yield 
was approximately equal to the yield of the thin 
target. This gives only a rough value of target 
thickness, but probably good enough for the 
present purpose which is merely to make sure 


that the use of even thinner targets would not 
give sensibly different data. Targets were re. 
placed after bombardment for less than half the 
time which experience showed to produce appre. 
ciable ‘‘fatigue.”’ 


RESULTS 


The relative yields at various angles, trans. 
formed to the center of mass coordinate system, 
normalized to unity at 90°, and plotted against 
cos* @ are shown in Fig. 2 for several bombarding 
energies. The number of particles counted at any 
window is 3000 or more for the window in ques. 
tion and 10,000 or more on the monitor. Straight © 
lines, 1+A cos? have been drawn through 
these points by inspection and the values of 4 
so determined are plotted against bombarding 
energy in Fig. 3. Vertical bars are estimated 
limits of error arrived at by drawing through the 
observed points of Fig. 2 straight lines of greatest 
and least admissible slope, admissibility being 
determined merely by the best judgment of the 
authors. 


DISCUSSION 


The appearance of Fig. 3 obviously suggests 
resonance to a rather broad excited state of Be’, 
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The Distribution in Angle of the Long Range Alpha-Particles from Fluorine 
Bombarded with Protons 


W. B. McLean, A. ELLETT AND JAMEs A. Jacops 
State University of Iowa, Iowa City, Iowa 
(Received June 28, 1940) 


The distribution in angle of the long alpha-particles from the reaction F!*+ H*+O"* + Het has 
been determined at bombarding energies of 330, 375 and 435 kev. This distribution shows a very 
strong fore and aft asymmetry which varies only slightly with the bombarding voltage. The 
distribution is well represented by the expression 1+0.66 cos 6+ 0.25 cos* 6+0.41 cos? @. 


HE present paper is concerned with the 
distribution in angle of the long range 
alpha-particles from fluorine bombarded with 
protons. The apparatus shown in Fig. 1 differs 
essentially from that described in the previous! 


1 Young, Ellett and Plain, Phys. Rev. this issue. 


paper only in the form of target support and in 
the dimensions of the defining apertures, which 
are as given in Fig. 1. Targets were prepared by 
electrolysis of hydrofluoric acid onto tantalum. 
Care was used to avoid boron contamination 
arising from the solution of Pyrex glass in the 
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ional vi f the target chamber and ait , , 
The canker is statistical error. Data secured with a thicker 


celerating tube by means of the greased joint shown at target and larger number of counts (between 2000 q 
the top of the figure. The target spot is approximately a 
rectangle 


hydrofluoric acid, as the boron reaction has at 
low energies a much larger yield than the v 
fluorine reaction. Observations were made at | 

330, 375 and 435 kev by the method outlined in wi 
the preceding paper, and the results are presented tH 
in Fig. 2. 


RESULTS 


The distribution shows a very large fore and 
aft asymmetry, with a strong preference for the 
forward direction. It appears to change very 
slightly with bombarding energy. It is true that ; 
the points taken at 330 kev and 435 kev do not ny 
coincide exactly, but their differences are hardly 
significant, as the number of counts was only 
about 500 on a particular point, giving a statistical 
probable error of 3 percent. These curves were 
taken with a relatively thin target (apparent 
width of the gamma-ray resonance was about 40 
kev) and since the distribution appeared to 
change but little with energy we decided that the 4 
use of a somewhat thicker target would be 
legitimate and would permit the reduction of the 


and 3000 per point) at a bombarding energy of 


RATIO TO 90° 


© 375 KEV 
O 435 KEv 


& 330 KEV 
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Fig. 2. The distribution in angle in the center of mass system of coordinates of the long range alpha- 
the reaction F!°+H+O!*+ Het. The triangles are for a bombarding energy of 330 kev, the circles for 375 kev, and the 
squares for 435 kev. The broken line represents the expression 0.985+0.656 cos @+-0.25 cos* @+0.40 cos* @ and the solid 
line the expression 1+0.65 cos 6+0.11 cos? 6+0.44 cos? 6+0.24 cos* 6. 
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Fic. 3. The ratio of the yield of alpha-particles at three 
backward windows (105°, 120°, and 135°) to three forward 
windows (45°, 60°, and 75°) as a function of energy. This 
curve is plotted in the laboratory system of coordinates. 


375 kev have been fitted with an empirical 
equation by the method of least squares. The 
best curve of the fourth degree is found to be: 


1+0.650 cos 6+0.11 cos? 6 
+0.44 cos* 6+0.24 cos‘ 6 


and is shown in Fig. 2 by the solid curve. The 
best third degree curve is found to be: 


0.985 +0.656 cos 6+0.25 cos? 6+0.40 


and is shown by the broken curve in Fig. 2. It is 
clear that the inclusion of the fourth-degree term 
leads to no significant improvement in fit. 
We‘have not been able to make similar ob- 
servations at lower energies because of the low 
yield. However, we determined the ratio of the 
number of alpha-particles backward to forward 
(windows at 45°, 60°, and 75° were connected to 
one amplifier and those at 105°, 120°, and 135° to 
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the other) down to 280 kev. This ratio, plotted ag 
a function of bombarding energy, is shown jp 
Fig. 3. 


DISCUSSION 


In view of the rather unexpected nature of this 
distribution we were careful to check possibje 
sources of error in the experimental arrangement. 
The foils were measured by means of a standard 
polonium source and a shallow ionization chamber 


and found to have a stopping power of 23. 


centimeters air equivalent. The outer edges of 
the collector plates were at a distance of 15 
centimeters from the foils so that all particles of 
range greater than 3.8 centimeters should traverse 
the full depth of the ion chamber. The pulses ag 
seen on the screen of an oscilloscope were of 
uniform height, indicating that all of the alpha. 
particles were traversing the full depth of the 
sensitive volume of the ionization chamber. Asa 
further check the alpha-particles from B+) 
were observed using the same target chamber 
and foils. Since these particles, which have a 
range about one centimeter less than that of the 
fluorine alpha-particles, were counted at all 
windows it was certain that none of the fluorine 
alpha-particles was being lost. The distribution 
of the alpha-particles from Li’+H! at bom. 
barding voltages below 100 kev was also observed 
with this apparatus and found to be yery nearly 
spherically symmetric in agreement with the 


results of the preceding paper. 
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On the Angular Distribution in the Reaction F" (9, «) 


VOLUME 58 


E, GEeRjuoy 


Department of Physics, University of California, Berkeley, California 


The angular distribution of the long range a-particles shows marked lack of spherical sym- 
metry, with asymmetry in the forward and backward directions. The distribution is constant 


with bombarding energy in the range 330 to 435 kev and contains large cos @ and cos? @ terms 
with somewhat smaller and less certain amounts of the third and fourth powers. It is shown 
that the energy dependence implies broad levels and constant penetrabilities; the latter is 
found to agree with calculation. The experimental distributions can most simply be accounted 
for by the assumption that the chief contributions are from three levels of total angular mo- 
mentum J=0,1,2. This leads to the experimental curves without arbitrary assumptions as to 
the values of the line breadths or detailed assumptions about the nuclear coupling scheme. 


Powers of cos @ higher than the fourth in the angular distribution are not to be expected. 


HE data obtained by Ellett and his co- 

. workers on the angular distribution of the 
long range a-particles from the proton bombard- 
ment of F!* show a marked lack of spherical 
symmetry with asymmetry in the forward and 
backward directions. A distribution containing at 
least the first and second powers of cos @ is 
required to fit the data at all well, and somewhat 
smaller and less certain amounts of cos* @ and 
cos‘ @ also seem to be present. More explicitly, 
the observed angular distribution is well repre- 
sented by a fourth-degree equation of the form: 


1(@)=1+0.65 cos 6+0.11 cos? 6 
+0.44 cos* 6+0.24 cos‘ @ (1) 


or probably equally well by: 
1(@) =1+0.66 cos @+0.25 cos? 6+0.41 cos* (2) 


These curves are calculated to fit the data for 
the bombarding energy of 375 kev but they fit 
the distributions for bombarding energies of 330 
and 435 kev almost as well, and in fact the 
quadratic equation : 


I(6)=1+0.81 cos 6+0.13 cos? 6 (3) 


fits the data for the three bombarding energies 
to within less than ten percent. 

On first sight these curves, drawn for the 
coordinate system in which the center of mass is 
at rest, present some surprising features, namely 
the pronounced asymmetry in the forward and 
backward directions, the large coefficients of the 
cos 6 and higher terms, and the marked con- 
stancy with bombarding energy. It is the purpose 
of this paper to show that none of these features 


503 


is inconsistent with our present theory of nuclear 
reactions, and that it is possible to fit the experi- 
mental curves by making reasonable assumptions 
about the unknown quantities such as line 
breadths which occur in the theory.* 

For any initial state of the system, the angular 
distribution of the final a-particle wave function, 
in the coordinate system in which the center of 
mass is at rest is: 


where the Y," are the spherical harmonics and 
the A," are complex coefficients which are 
determined by the initial state i and by the 
selection rules conserving the parity, total angular 
momentum J, and total magnetic quantum 
number m, as well as by the probabilities of 
capture and emission of the respective particles. 
The total angular momentum of F* in its 
ground state is 4. The resultant angular dis- 
tribution, obtained by summing (4) over the 
statistically independent initial states 1, may then 
be found by summing over the possible different 
orientations of the spin S, of the incoming 
proton and the total angular momentum J of 
the F'* nucleus. We shall adopt the procedure of 
summing (4) over the four independent eigen- 
functions of the operator T=J+S,, the sum of 
the spin of the proton and the total angular 
momentum of the nucleus, indexing these states 

*To calculate these unknown quantities explicitly, a 
detailed knowledge of the coupling scheme in the nucleus 
would be required. We shall fit these quantities to the 
experimental curves, and our calculation will then be 


carried through without any assumptions about the 
coupling. 
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by 4. Furthermore, taking the parity of the 
ground state of F!® as even, because JT is at 
most 1, and the final spin is 0, the disintegration 
into an a-particle of orbital angular momentum 
L=!I can only arise from the capture of a proton 
of equal Z to form a compound state of Ne?* 
of J=/. 

It is, in principle, possible for any number of 
levels of J=/ of Ne*** to contribute to A,” in 
(4) and it is clear that a complicated picture of 
the reaction of this sort could readily lead to the 
distributions (1), (2) or (3). We shall assume 
simply that the process is one of resonance with 
two, or at most three, broad overlapping levels 
of the compound nucleus Ne®* of total angular 
momentum J=0, 1, 2, and therefore write: 

(T pimil aim 


In (5) ax” is the amount of ¥(7ilm) in the 
initial wave function, where ¥(7LJm) is an 
eigenfunction of J expressed in terms of the 
initial eigenfunctions of T and L. T; and E,; are 
the total width and resonance energy of the 
level J=/ of the compound nucleus, I'pin; and 
Taim the decay constants from this state to 
vu” and the final Y,", respectively. The wim; are 
determinate, though unknown phases, and finally 
we may write I'pyimi aS Pprypimi the product of 
P 1, the penetrability of the barrier for the state 
whose wave function is y¥,;" and the “partial 
width without barrier,”’ 7 pimi, which may also be 
regarded as the probability of forming the 
particular compound state once the particle has 
penetrated the nucleus.f{ It follows from (4) and 
(5) that the constancy of the distribution with 
energy will be highly unlikely unless the penetra- 
bilities for the various /’s are constant with 
energy. This necessary condition is satisfied as 
seen from Table I which summarizes the calcu- 
lated penetrabilities for the three different 
t The a,” are calculated with the aid of the Clebsch- 
Gordan series and depend on the order in which the 
eigenfunctions of the three operators J, S, and L are 
bine first J and S, 


combined. We have here chosen to com 

but a different initial pair would lead to somewhat different 
numerical coefficients in (6) below and consequently 
require a different relation (7) to fit the observed distribu- 
tion. The physical results are, of course, the same and 
this is simply a consequence of the fact that the different 
orders of combining the operators correspond to a different 
choice of the set of independent states of the initial 


system. 
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TABLE I, Calculated penetrabilities of F, 


— 

E RxXx103 Ppi-t/Ppi-o Ppi~t/Pying 

330 kev 4.00 0.106 0.002 
375 4.00 0.110 0.003 
435 4.00 0.112 0.003 
330 3.52 0.092 0.002 
375 3.52 0.094 0.002 
435 3.52 0.096 0.002 


energies and for two different values of the 
radius of 

In addition, in order to explain’ particularly 
the constancy of the interference terms in (1), (2) 
and (3) with energy, it is necessary that the 
phases of the A,;" remain essentially constant in 
the energy range considered, or at most vary go 
that the phase differences remain constant. It js 
this requirement which makes necessary the 
assumption that the levels are broad. We wish 
to stress that this requirement and the absence 


of marked resonance structure in the yield curye - 


for this reaction at these energies, means that 
the assumption of resonance with three levels 
implies nothing more than that these levels 
contribute larger terms to the cross section than 
other levels. The assumption is by no means 
necessary and is justified only by the fact that 
it leads to the observed distribution without 
arbitrary assumptions about the values of the 
line breadths. In addition, in order to simplify 
the calculation, we shall assume at once that the 
level widths I; in (5) are approximately equal. 
We now note the following: The Clebsch- 
Gordan series gives a;;°=0 for 7 corresponding to 
the initial state 7=1, m=0. This, coupled with 
the selection rules, will mean that we obtain 
interference between P and D waves in the final 
distribution for incident states T=1, m=1 and 
—1, and interference between all three S, P, D 
waves for the 7=0 state. Or, alternatively, 
that there is interference between P and D 
when the z components of the angular momenta 
of the incident proton and F’* nucleus are 
parallel, and interference between the three 
waves when the z components are antiparallel. 
Moreover, we can now drop the i subscript in (5) 
and index the terms by their m values, m=0 
corresponding to an initial T=0 state and the 
other m values to initial T=1 states. The Tam 
are independent of m while the ypim and @im 
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depend only on whether the initial state was of 
T=1 or 0. With our assumptions the denomi- 
nator of (5) is approximately the same for all 
three levels and may be canceled from (4). The 
Tem will form the largest part of T; and our 
assumption of approximate equality for the T's 
implies approximate equality for the T'aim’s as 
well, and these may also be canceled. Now insert 
in (5) the calculated values of the a and the 
penetrabilities for E=330 kev, R=4X10-, put 
wo=0, and equate the m=—1 to the m=1 
terms. If we substitute and sum over m in (4), 
the resultant angular distribution is: 


1(0) 0.477711 +0.01 2520 
Cos weo 
+cos cos 
—0.220(y:0720)! cos (w20—w10) 
+0.381 (yi1721)! cos (wei —w11) 
+cos? 
Cos woo 
+cos* 6[0.660(y10720)? cos (w2o—w10) 
— 0.381 (y11721)! cos (we1—w11) J 

+cos* 6[0.1125y20—0.075y21]. (6) 


All unnecessary subscripts have been omitted. 

To fit the simple quadratic distribution (3) we 
need only assume the existence of two over- 
lapping levels of total angular momentum J=0 
and 1. Dropping the terms in 20 and ya: in (6) 
we find that with 


Y11=10= 00/3, (7) 
the resultant angular distribution will be 
I(6) =1+0.97 cos w39 cos 6+0.13 cos? 6. 


The relation (7) is satisfactory because it is not 
improbable and is what might be expected from a 
nucleus in which Russell-Saunders coupling is not 
valid and the states of equal J are well mixed. 
Clearly (6) can be made to fit the observed 
distribution if we assume other sets of over- 
lapping levels. These, however, will in general 
require rather more arbitrary assumptions about 
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the values of the yim’s. In particular, the assump- 
tion of S and D resonance will lead to a distribu- 
tion with a negligible cos @ term, and the 
assumption of P and D resonance requires 
juggling of the phases in order to make the cos @ 
term larger than the cos* @ term. In addition, 
with either of these assumptions, as well as for 
the assumption of the overlapping of all three 
levels, or S, P, D resonance, the cos‘ @ term 
cannot be fitted without always assuming the 
‘Yom’s rather larger than the other yim's. 

Because of the numerous assumptions in- 
volved, it is not possible to say any more than 
that the results of this paper suggest the follow- 
ing conclusions: The possibility of fitting the 
curve so well with the distribution (3) and of 
deriving this distribution without making arbi- 
trary assumptions as to the values of the level 
widths seems to indicate that the observed 
distribution arises chiefly from a resonance, in 
the sense indicated, with two broad overlapping 
levels of J/=0 and 1. The difficulties encountered 
in attempting to fit the cos‘ @ term might indi- 
cate that this term and the cos*@ term are 
actually somewhat smaller than observed, but 
it is quite possible that these terms arise from 
variation in the level widths among these three 
levels or even in a more complicated way from 
the contributions of other more distant levels. 
It should be noted that the possibility of large 
P scattering arises from the fairly large value of 
the penetrability, combined with the normaliza- 
tion factors of the spherical harmonics and the 
weights corresponding to the more numerous 
possible orientations of the particles of higher 
orbital angular momentum in the incident beam. 
These factors are not quite sufficient to com- 
pensate for the decrease of penetrability for the 
D wave, and the increasingly rapid rise of the 
potential barrier for waves of higher L make it 
unnecessary to consider these waves of higher L 
and similarly they should not appear noticeably 
in the angular distribution. 

I wish to express my thanks to Professor J. R. 
Oppenheimer and Dr. L. I. Schiff for suggesting 
this problem and for much encouragement and 
helpful discussion. I should also like to thank 
Dr. Ellett for sending me his data in advance of 
publication. 
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PHYSICAL REVIEW 


The Interpretation of Resonances in Nuclear Reactions 


G. Breit 
University of Wisconsin, Madison Wisconsin 
(Received July 15, 1940) 


The theory of resonances in the one-body problem model 
is systematized. Energy shifts in the apparent position of 
levels are found. The apparent position of levels in scatter- 
ing and in y-ray emission is expected to be the same to 
within ~T*/E in the usual notation if the only repulsive 
barrier is the centrifugal one. In the absence of strong 
repulsive barriers there may be additional shifts of the 
order I’. Interference between levels is shown to occur for 
the scattering cross section in the numerator and the de- 
nominator of a fraction in accordance with Eqs. (3.3), 
(3.6). The artificiality of usual formulas with interference 
for cross sections is discussed in this connection. The 
construction of Green’s function for two-dimensional 
separable problems is described. The construction works 
also for some many-dimensional problems. The reduction 
of the number of essential dimensions is outlined by a 
reformulation of the problem in one less dimension. Ap- 
plications are made to the solution of scattering problems 
which represent schematically idealized nuclear problems. 
The single particle potential barriers enter the solutions 
through the regular and irregular solutions f, g of the 
radial equation. It is seen from the solutions that the 


yield may depend on the f and g for various excited states 
of the residual nucleus. A simplification occurs if the inter. 
action is highly repulsive and is localized in the two- 
dimensional space of the distances of the particles from 
the center. In this case the regular functions f for the 
incident and final state enter linearly in the wave function 
as in Eq. (15.5). If the interaction region is made large and 
if the incident state happens to be especially important 
in the expansion of Green’s function then the factor g of 
the incident wave occurs in the formula for the cross 
section as in Eq. (16.4). A resonance model is worked out 
by assuming the interaction to be highly attractive within 
a circle and repulsive in a ring surrounding the circle in the 
two-dimensional diagram. The damping constants in the 
resultant formula contain f? as in Eq. (18.8). Here also 
the single dependence on f? is a consequence of having 
localized the interaction to a narrowly defined radius for 
the incident particle. In the general case each damping 
constant involves the g of different states. The examples 
show that quantitative applications with simplified forms 
of damping constants and with interference types of dis- 
persion formulas have only a limited validity. 


INTRODUCTION 


ANY nuclear reactions show excitation 
curves with pronounced peaks which sug- 
gest that there is a resonance of the nuclear 
system to certain energies. The older discussions 
of such resonances have been concerned mainly 
with “one-body” models in which the actual 
system was schematically represented as a 
system consisting of the incident particle moving 
in the field of the bombarded nucleus. The large 
number of resonances discovered in the bombard- 
ment of nuclei by slow neutrons as well as other 
reactions have shown the inapplicability of the 
“one-body” picture in the majority of cases and 
it has become clear that many nuclear particles 
are usually involved in a resonance state. 
The general physical aspects of the problem 
are clear and well known.'-5 The bombarded 


1N. Bohr, Nature 137, 344 (1936). 

2G. Breit and E. Wigner, Phys. Rev. 49, 519 (1936). 
Referred to as BW. 

2H. A. Bethe and G. Placzek, Phys. Rev. 51, 450 (1937). 

‘F. Kalckar, J. R. Oppenheimer and R. Serber, Phys. 
Rev. 52, 273 (1937). 

5H. A. Bethe, Rev. Mod. Phys. 9, 71 (1937). 


nucleus together with the incident particle can 
form a state which has an appreciable stability 
and a relatively long mean life. The width of 
the state is connected with the mean life by the 
indeterminacy relation so that a long life corre. 
sponds to a sharp level. The compound state is 
formed with a maximum intensity when the 
incident particle has a suitable energy and with 
smaller intensity when the energy of the bom- 
barding particle is varied through the resonance 
width. The compound state disintegrates in 
several energetically possible ways producing 
either scattered incident particles or disintegra- 
tion fragments. 

This picture is good enough for a qualitative 
understanding of some reactions. It is not entirely 
satisfactory for quantitative applications. The 
compound state used in the above description is 
not defined by the description itself and it is 
understood in different ways by different physi- 
cists. Some think of it in strict analogy to pre 
wave-mechanical pictures as being like an e- 
cited state of the electronic system of an atom 
which can emit photons of several different 
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wave-lengths by making optical transitions. This 
picture is correct except for its pre-wave- 
mechanical features. Without wave mechanics 
one can describe qualitatively the decay of the 
atomic level by introducing Einstein’s emission 
probabilities and saying that the atomic system 
jumps to one or another of the lower levels with 
the simultaneous emission of a photon. So far 
the picture is correct. It is then said that the 
level width is increased because it is not a real 
stationary level and that the emitted lines will be 
fuzzy. For a justification of the latter view an 
appeal is made to the indeterminacy principle. 
At this point there is a bad mix-up, the word 
“level” being used in two different ways. In the 
first part of the picture it was used in the sense 
of an atomic level, the atomic system being 
pictured as an idealized mechanical unit with its 
interaction with radiation neglected. This kind 
of “level” is a clear concept as long as it lies in 
the discrete part of the spectrum, as it is sup- 
posed to lie. It is not diffuse but absolutely 
sharp. When, toward the end of the explanation, 
one speaks of the “‘level” as diffuse, it is the 
level of the compound system consisting of 
atom+radiation field that one is thinking about. 
This kind of level is not well defined in the above 
explanation. 

The actual situation in the optical case has 
been cleared up by Weisskopf and Wigner.‘ 
They showed by a perturbation calculation that 
if an atom is put into an excited state j, so that 
it is definitely in that excited state at the time 
t=0, then at later times there is (a) an expo- 
nentially decaying probability of the atom being 
in the excited state j, (b) a growing probability 
of lower atomic states 7, k, / and of photons with 
corresponding energies. The energies of the 
photons emitted in the ji jump are approxi- 
mately equal to the energy difference E; — E;= Ej; 
between the upper and lower atomic states. The 
band of photons emitted in such a jump is found 
to have an intensity distribution of the resonance 


type 
const. 
Te 
*V. Weisskopf and E. Wigner, Zeits. f. Physik 63, 54 


(1930); P. A. M. Dirac, Zeits. f. Physik 44, 594 (1927); 
V. Weisskopf, Ann. d. Physik 9, 23 (1931). 


(1) 


with a peak at the emission center E;;. The 
half-width of the line 2(1;+1T.+TI) was shown 
by them to be the sum of half-widths due to 
separate emissions ji, jk, jl with emission proba- 
bilities ,./h, They have also 
explained that for the compound system atom 
+radiation the energy of the initial state (no 
photons, atom in state j) is not sharp but has a 
width connected with 


(1’) 


by the indeterminacy relation. In their theory 
the ‘‘compound state’”’ is the state of the system 


atom-+radiation field 


in which there are no photons and the atom is in 
the state 7. This compound state does not have a 
definite energy. The energies of the compound 
system out of which the state is composed have 
a probability distribution of the form (1). 

In the optical problem the optical transition 
probabilities are small and it is natural as well 
as useful to define the compound state in the 
above manner in most cases. But even here this 
definition is not the most natural under all 
circumstances. If a group of atomic levels E;, Es, 
etc. have energy differences etc. 
smaller than the energy widths T, Ts, etc. the 
compound state obtained by exciting the atom 
to a single level, say E2, and having no photons 
at ¢=0 does not correspond to a simple physical 
condition. The photons emitted from E, are re- 
absorbed so that the atom is partly in E, and 
partly also in other states E,, E;---. It is here 
also possible to define? the compound state in a 
sensible manner, at least approximately. This 
corresponds to an initial state of the atom with 
a wave function which is a proper linear com- 
bination of the wave functions y¥, corre- 
sponding to E,, Es, E;---. If, at t=0 the atom is 
put into such a state, represented by ci¥i+cowe 
+---, and if initially there are no photons, then 
at later times the probability of the atomic 
state - - decreases exponentially and 
photons appear within a band having a width 
connected with the mean life of the atomic 
state by the indeterminacy relation. It is to be 
noted that the atomic states E,, Es, --- have 


7G. Breit, Rev. Mod. Phys. 5, 91, 117 (1933); G. 
Breit and I. S. Lowen, Phys. Rev. 46, 590 (1934). 
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different energies. One could, if one liked, speak 
of ¥2 as the compound state. This would not be 
an especially useful procedure for physical appli- 
cations, however, because (a) this compound 
state would not decay exponentially; (b) the 
scattering of light is expressed more simply in 
terms of the states c:~itcepet+---. 

It is seen that even in the optical case the 
definition of the compound state requires caution. 
The natural and useful definition is that of using 
‘linear combinations of atomic states which give 
an exponentially decaying probability of atomic 
excitation. This definition is slightly more com- 
plicated than that for the problem of Weisskopf 
and Wigner. It is definite as long as the inter- 
action of the atom with the radiation is suffi- 
ciently weak. If the interaction were strong the 
exponential decay of the atomic excitation would 
not be obtained and the compound state could 
not be defined in the above manner. 

Nuclear resonances differ from the optical 
problem in the following respects: (a) There are 
cases of nuclear resonances in which the peaks 
of excitation curves overlap. The coupling 
through emission which has just been discussed 
can be expected, therefore, to be of more fre- 
quent occurrence. (b) In the optical problem the 
direct change of one type of photon into another 
is of secondary importance. With Dirac’s electron 
equation it even disappears altogether. The 
scattering of photons in this formulation can be 
explained entirely by a succession of emissions 
and absorptions affecting the wave function. In 
the general nuclear problem such a simplification 
cannot be made. Direct transitions between 
states in the continuum are possible. (c) The 
effect of a nuclear resonance is of interest also 
outside the main resonance peak. For slow 
neutrons this effect gives the 1/v law. There is 
evidence for it also in the y-rays emitted by Li’ 
under proton bombardment below the 440-kev 
resonance. (d) The interaction of the atom with 
radiation is always small. The level of the com- 
pound system responsible for the resonance is a 
state in which there are no photons present. 
The interactions between nuclear particles are 
not sufficiently small to make the analogous 
simplification possible. 

It is desirable to be able to use experimental 
yield-energy curves and to draw from them 
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conclusions regarding properties of the nuclear 
system. Data on widths of nuclear levels due tp 
neutron emission should be comparable with 
similar data for widths due to proton emission, 
The comparison can be made only if 
allowance is made for the effect of the Coulomb 
barrier. A quantitative understanding of the 
problem is necessary for this purpose. Similarly 
it is desirable to know whether peaks in yield. 
energy curves indicate the position of levels of 
the compound nucleus with any precision and 
just what relation these energy levels have to 
ordinary nuclear levels. For example, the com. 
pound state (C”)* formed as B"+H! shows 
resonances in gamma-emission at proton energies 
that correspond to a mass difference between 
(C)*—B” which would be expected if ag 
well as B” were stable and if the binding energy 
of (C)* differed from the binding energy of Be 
only on account of the Coulomb energy. A ques. 
tion of direct physical interest is whether the 
hypothesis of a symmetric Hamiltonian for 
nuclear particles makes one expect a coincidence 
of this sort and to what accuracy. Another 
type of question arises in observations of the 
same level under different circumstances. The 
position of the level is not expected to be 
exactly the same and some calculations regarding 
the magnitude of the energy shift are made in 
the present paper. 

For the discussion of such questions it js 
difficult to use too general and formal treatments 
such as have been attempted by Bethe and 
Placzek* and by Kapur and Peierls.* The point 
of view of these authors is that one must repre- 
sent the collision cross section of a reaction as the 
square of the absolute value of a sum of dis 
persion-like terms. Even if this were accom- 
plished the result would not necessarily be useful 
for the above types of questions because it is 
impossible to compute all the terms in such 
formulas starting with expressions for nuclear 


8’ P. L. Kapur and R. Peierls, Proc. Roy. Soc. 166, 277 
(1938). For the one-body problem at low energies this 
theory gives levels which correspond to d§/dr=0 at the 
nuclear radius. There are no perceptible resonances at 
such energies in the deuteron problem. This is due to the 
fact that the KP levels depend on the total energy and 
the caution, which must be used in attributing to their 
levels ordinary resonance properties, may be em 
These KP levels show —— resonance to change of 
depth of the potential well but not to change of total 
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forces. It is practically hopeless to analyze data 
in terms of resonance levels having variable 
positions depending on the energy of incident 
particles, as well as on the assumed nuclear 
radius, as becomes necessary in the formulas of 
Kapur and Peierls. 

The unsatisfactory state of the theory of 
nuclear resonance makes it desirable to consider 
some typical examples. Only partial answers to 
some of the above questions will be arrived at 


below. 
OnE-Bopy PROBLEM 


(A) Method of complex eigenvalues ‘ 


For the one-body problem a single and reason- 
ably sharp level in the continuum has simple 
properties. The absolute value of the wave 
function, normalized to unit amplitude out- 
side the nucleus, becomes a maximum inside 
the nucleus for the resonance energy. An inci- 
dent plane wave of particles asymptotic to 
exp (ip cos @) and ‘scattered by a central field is 


represented by 
v=> 14(2L+1)Pr(cos 0)Fx1(p)/p. (1) 


The functions §}z(p) are normalized to unit 
amplitude at p= and, therefore, a sharp 
maximum of §z inside the nucleus corresponds 
also to a maximum of y. The asymptotic form 
of Fx is 

sin 


For a given L a resonance can be expected in 
the neighborhood of a complex eigenvalue of 
the energy. The possibility of using complex 
eigenvalues of the energy has been first intro- 
duced into quantum theory by Gamow.® A con- 
nection with resonances has been described by 
Breit and Yost.'° The latter connection has 


(t>0); ¥=< 
(t<r/v), 


(t<0); 


(t+r/v9>0). 


24 (Ea) exp | now) | 
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recently been pointed out again by Siegert." 
A discussion of the one-body resonances making 
no use of complex eigenvalues has been given by 
Wigner and Breit.? Kapur and Peierls* have given 
a discussion of one-body resonances also using 
complex eigenvalues but along different lines. 

By the method of complex eigenvalues one 
considers the asymptotic form of the solution of 


@V2/dr?+ (2u/h*?)(E— V)¥e=0, (2) 
which is 
We~Ae*r+ Bet", (2.1) 
The normalization of Vz is such that 
is kept constant for all energies E. If 
(dVz/dr),.o=0 
then one uses for C the first nonvanishing 
derivative of Vz. The function A(£) is supposed 
to have a root Ey—iIy so that 
A(Eo—iIl) =0. (2.2) 
For small values of E—E)+iIy one has 
A(E) 4-0, 
B(E) 


It is assumed that for real E, A*=B. The con- 
nection with Gamow’s picture of an exponentially 
decaying state is obtained by forming a wave 
packet out of the solutions of (2) with real values 
of the energy. One has then as in Eqs. (20) of BY 
the following wave packet 


To exp (—iEt/h) 
Ve -dE=W. (2.4) 


This wave packet outside the nucleus i.e., in the 
region where (2.1) applies is given by 


(2.3) 


(2.5) 


(2.6) 


*G. Gamow, Zeits. f. Physik 51, 204 (1928); H. Casimir, Physica 1, 193 (1934). 
“G, Breit and F. L. Yost, Phys. Rev. 48, 203 (1935). Referred to as BY; G. Breit, Phys. Rev. 40, 127 (1932). 


"A.J. F. Siegert, Phys. Rev. 56, 750 (1939). 
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Here vp is the classical value of the velocity of 
the particle outside the nucleus and is obtained 


from 
Eo = 


According to Eq. (2.5) the wave packet flows 
out of the nucleus for ¢>0 as a diverging wave. 
The place at which ¥ becomes zero moves with 
the constant velocity vo. For a fixed r the 
amplitude of WY decreases exponentially at the 
rate exp (—Iot/h). The sharp break in W at 
r=vol is not strictly correct. It is brought in by 
the approximation in the calculation of BY in 
the neglect of terms in (E—£o)? in the expansion 
of k in powers of E—Ep. For t<0 one obtains 
only incoming waves as indicated by —ipvor/h 
in Eq. (2.6) and for a fixed r the function in- 
creases exponentially with ¢ up to t=rvo. After 
this time ¥=0 at this r. The wave packet (2.4) 
represents particles flowing out of the nucleus for 
t>0 and flowing into it for <0. 
The mean life of the nucleus is r given by 


1/r=21o/h=4nT o/h, (To in energy units) (2.7) 


1/r=4nTo. (Io in frequency units) (2.8) 


In energy units Tp is the coefficient of —7 in the 
complex eigenvalue. This is convenient and this 
quantity will not be denoted here by y/2 as is 
often done. The value of I) can be computed in 
some cases by finding the root of Eq. (2.2). 
This is not necessary, however, because one has 


(Eq. (21”) of BY] 


where §£o is the real solution of Eq. (2) for 
E=E, normalized to have unit amplitude at «. 
The integral in Eq. (3) is taken through the 
nucleus and through a part of the barrier outside 
the nucleus along the rapidly decreasing portion 
of §. The qualitative meaning of Eq. (3) is clear 
when it is noted that the numerator is the 
velocity of the escaping particle while the de- 
nominator is the probability of the escaping 
particle being inside the nucleus. The correctness 
of the numerical factor 1 in the formula is not 
immediately obvious since the real function 


(3) 
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represents a standing wave and corresponds to a 
density (per unit thickness of a spherical 

of } rather than 1. The correctness of the facto, 
1 is shown in BY. It may also be seen from the 
fact that inside the nucleus there is little differ. 


ence between §£0() and the solution of Eq, (2) 
for E=E,)—iT» normalized so as to make 


2|A(Eo)| =1 (3.1) 


in Eq. (2.5). If we assume for the moment that 
this is the case, the coefficient 2A*(E») which 
multiplies the exp { } in Eq. (2.5) has absolute 
value one. The rate at which particles flow oy 
for this wave packet is vo and Eq. (3) follows. The 
fact that Jo"|}|*dr is practically the same fo, 
the normalization (3.1) of the complex eigen. 
value solution as for and 
follows from the assumed sharpness of resonance. 
This in fact makes /("|§|*dr dependent prag. 
tically only on C and not on E as long as E js 
not changed by amounts greater in absolute 
value than Io; and the real function Fap(p) 
normalized so as to be a sine function of unit 
amplitude at p= satisfies Eq. (3.1) so that 
WEo= £0 makes the coefficient of the exp in 
(2.5) have absolute value 1 in agreement with 
Eq. (3). 

The correctness of the coefficient 2A*(E,) in 
(2.5) can be seen without integration if one 
observes that this coefficient must be B(E,—iI\) 


and that 
B(Eo—il 9) = A*(Eo—ilo+2il 9) =2A*(E)). 


Equation (3) is definite only if the upper 
limit can be defined without ambiguity. For 
cases where this cannot be done a more precise 
formulation has been given by BW. This makes 
no use of the complex eigenvalue point of view 
and allows a better definition of the upper limit. 
It follows from Eqs. (2.3) that for the wave (1" 
the probability of finding the particle inside the 
nucleus is proportional to 


1/[(E—Eo)?+T 0]. 


In fact outside the nucleus the normalization of 
Eq. (3.1) holds. Hence according to (2.3), |C} 
is proportional to (3.2) and ||? inside th 
nucleus is, therefore, also proportional to (3.2) 
A nuclear reaction which is proportional to the 
probability of finding the incident particle inside 


(3.2) 
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the nucleus is also proportional to this expression 
and is thus represented by a resonance formula. 
According to Eqs. (2.3) the phase of § changes 
by ~™ as E changes from Eo—nI to Eo +nTo 
for n appreciably >1. One goes through a reso- 
nance for scattering in the vicinity of Eo. The 
explicit connection with the direct calculation of 
shifts is seen through Eqs. (9), (9’), (10’), 

(13), (17) of BY. The discussion of BY and of 
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discussion of Siegert brings in a distinction 
(mentioned in reference 13) between resonance 
and potential scattering which is rather artificial 
and does not correspond to the usual one. It 
may be permissible, therefore, to consider the 
connections again and to extend them so as to 
include the combined effect of several complex 
eigenvalues for the energy. It may be possible 
to represent the function A(EZ) defined by Eq. 


BW at this point is perhaps too brief and the (2.1) as 


= +E ( + ; ) (3.3) 
A(E) a(E) A'(E,—i0',) \E—En+ily 


and to have 1/a(E) a sufficiently smooth function of the energy to make its more detailed considera- 
tion unnecessary.” The chance of finding the particle in the nucleus, according to Eqs. (1”) and (2.1) 


is proportional to 
|1/A(E) |? 


and shows, therefore, the interference of N terms of the resonance type. For special functions A(Z) 
it is possible to have a(Z)=const. and N= ~. The analytic character of A(E) is not always easy to 
ascertain and it appears preferable, therefore, not to assume that a(Z) is a constant. The representa- 
tion of 1/A(E) by the sum on the right amounts then to splitting 1/A(Z) into a smoothly varying 
part 1/a(EZ) and a part represented by >, which is designed to take out the bumps of 1/A(E). The 
partial scattering cross section corresponding to angular momentum L is 


oL= (2L+1)(A?/7) sin? Kt. (3.4) 
The phase shift Kz, is obtained from (2.1) which yields 
exp {i(K,—Lx/2)} =iB/|B| =iA*/|A| 


so that 
(2L+1)A A* Ay (2L+1)A? |1/A+(—)4/A*|? 
on 2 ) = 3.6 


Equation (3.6) shows that the partial scattering cross section oy is a fraction in which both the 
numerator and denominator are sums of the type of Eq. (3.3). This form of the scattering cross 
section is different from the formulas frequently given in which ¢, is represented as the square of 
the absolute value of a sum of terms of the type 


‘ (3.7) 


an 
| E-E,'+:T,’ 
Even if the form (3.6) should be reducible to (3.7) there is no reason to expect E,—iI’, to be equal 
to E,’—iI’,’. The representation of partial cross section by “dispersion theoretic formulas” is seen 
to be somewhat arbitrary in the case of interfering levels. The interference is seen to take place in 
two sums occurring in the numerator and the denominator of a fraction separately as in Eq. (3.6) 
rather than in a single sum such as is written in the form (3.7). 


® Professor G. Y. Rainich has called the author's attention to the fact that for entire functions a may be set =const. 
and the summation is then an infinite series. 
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In the vicinity of a single sharp resonance level one may write 
1 1 
A(E) B(E) 


and consider 1/8(£) a smooth function of E. Since the whole scattering is due to the potential V no 


reasonable distinction between 1/8(Z) and 1/A(E)—1/8(£) as referring to ‘‘potential’’ and “reso. 
nance”’ scattering can be made in this problem.” Both terms are in this case due to potential scatter. 


ing. For large B(Z) one has 


(3.8) 


A? —@n(E— Eq) }? 


9(2L+1)— 


where 

1/A’(E,—iT =an—tb, 
and a similar formula for odd L. Thus even neglecting 1/8(£) in (3.8) one obtains a shift of the 
maximum of o, with respect to the maximum of 


1/((E—E,)?+T,?]. 
If we neglect the variation of A with E the maximum of ¢, is at an energy 


and is displaced with respect to E=£, by an amount of the order I',. There is thus no exact agree. 
ment between the maximum of the scattering cross section and the maximum of the probability of 
the bombarding particle being in the nucleus. The two maxima are displaced relatively to each 
other by an amount of the order of magnitude of the width of the level. A displacement of the same 
order of magnitude can be expected between the maximum of the y-ray yield curve and the maximum 
of the scattering cross section. Additional shifts are produced by 1/8(E£) which has been neglected in 
Eq. (3.9). In the case of Coulomb scattering there are additional shifts which are due to the inter. 
ference of the Coulomb wave with resonance scattering. These have not been considered above. 


(B) The method of regular and irregular functions 

The method of complex éigenvalues has some advantages in introducing from the beginning the 
quantities E,—iI, which appear also in the final formulas. It is not a convenient method, however, 
for the computation of special problems. The values E,,—1I’, are troublesome to find and an analysis 
of the series in Eq. (3.3) by Cauchy’s theorem is complicated. The possibility of a continuous dis- 
tribution of poles in 1/A(E£) may also enter. For physical applications one needs 1/A(E£) only for 
real E. For this reason the method of complex eigenvalues has been partially related to the straight- 
forward use of wave functions for real E by BY to the extent of Eq. (3) of the present paper and by 
BW in a more thorough manner. The frequency shifts, discussed at the end of the preceding section 
have not been covered in these papers and will now be considered. 

The value of § , outside the potential well (r>a) is obtained from the requirement of the con- 
tinuity of §} z and its derivative. This requirement gives 


1 1 
—=—(1 F,G16—iF =—(G,+1F_) (Gx'+iF 


Fr'=dF./dp, 
18 Such a distinction is made in the one-body problem by Siegert, reference 11. 
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The second of the three forms for 1/§z follows from the invariance of §1/(Gi+iF 1) —§1(Gi'+iF 1’) =1 
as is seen for large p. The functions F;(p), Gi(p) are the regular and irregular solutions of the wave 

equation 

@F, L(L+1) 
dp? 


2u 
=0; k?=2uE/h? 


defined by their asymptotic forms sin (p—Lx/2), cos (ep—L/2). The only potential barrier con- 
sidered for these functions is that due to the term —L(L+1)/r* in the wave equation. The nuclear 
potential is supposed to be zero for r>a. The maximum of §, for a high potential barrier corre- 
sponds approximately to 


The first and second forms of the above equation are equivalent on account of 
(4.2) 
and both give for the resonance of Gyeate/A?* 
(4.3) 


This relation means that at resonance the wave function outside the nucleus is a multiple of the 
irregular function G,. For the resonance conditions given by Eqs. (4.1), (4.3) one has 


tan K,;=, sin? K,;=1. (4.4) 
These conditions correspond exactly to maxima of 
o1,/A? (=maximum) (4.5) 


and for sharp resonance give also maxima of the scattering cross section. Equations (4.1), (4.3) have 
been used by BY and BW both in the above sense and as an approximate way of finding the reso- 
nance energy of the maximum of | ,|* inside the nucleus. For future reference there is given below 
also the formula for Kz, 


1—F 


The value of §1’/}x at r=a is determined by the wave equation and the boundary condition at 
r=0. Equations (4.6) and (4) may be used, therefore, for the calculation of tan Ky, and 1. The 
invariance of 1’ 1, (Wronskians) for changes in p together with Eq. (4.2) 
shows the correctness of the second form of Eq. (4) and the first form of Eq. (4.6) practically without 
calculation. One has from the asymptotic form for 


Kr+Gzsin Kz ]exp (iKz), (r>a) (4.7) 


which can also be verified: by means of Eqs. (4.6) and (4). 
For high barriers (i.e., high values of L(L+1)/r?) one has 


Gi(ka)> F (ka). (5) 


Equation (4.7) shows then that the maximum of | (ka) |? should fall approximately at cos K,=0, 
sin K,=1. From Eq. (4) one obtains the condition for a maximum or minimum of ||? 


Resonance to scattering maximum of o,./A? occurs for F,Gz,é=1. At the energy for which this 
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happens one has 


= 


aE (5.2) 


and the maximum of the wave function inside the nucleus is not reached at the same energy as the 


maximum of o,./A?. An estimate of the difference in energy for the two maxima will now be made. 
From the second of the three forms for 1/9 given in Eq. (4) one obtains 


One has the relations 


OF, k? OF k* 
) f F,2dr, — =— f (5.4) 
F.0r EF;? 0 Es z* 0 
and the resonance width is determined by 
—_ = -ar 
ROE Gg 


with AE determined so that at a distance AE from the resonance energy |x|? has one-half its maxi. 
mum value. In the derivation of Eq. (5.5) it has been assumed that the resonance is so sharp that 
the maximum of | § |? occurs at the same energy as the maximum of ¢,./A*. It was found by means 


of Eqs. (5.4) and 


which is suggested for sharp resonance by the last of the three expressions for §, in Eq. (4). At 
E=E,+AE the scattering cross section is also approximately } of its maximum value as is seen from 
the fact that (5.5’) makes tan K,=1. For large |G,/F | the scattering cross section and |¥z(ka)| 
are nearly proportional to each other. This is seen from Eq. (4.7) which gives 
[1+(F1/G1) cot Ki sin? (5.6) 
The condition which must be satisfied for the proportionality of the scattering cross section in units 
A? and the probability of finding the particle in the nucleus is 
| cot Kr const. (5.7) 
Close to the resonance peak one may neglect (F,/Gz) cot Kz. The ratio is then G,? which varies 
with the energy. For high barriers G;? is likely to decrease with the energy. The scattering cross 
section is relatively higher, therefore, on the higher energy side of the peak around its main part. 


There is thus no exact correspondence between the maximum chance of the bombarding particle 
penetrating the nucleus and the scattering cross section and the correspondence becomes poor when 


|tan 
on account of the first ee in Eq. (5.6). By means of Eqs. (5.4) one obtains 
F,! 
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In these formulas E, is the energy for which the scattering resonance condition (4.1) is satisfied. 
The derivatives 0/0E are supposed to be evaluated for E= Ey. The same convention will be followed 
in the following approximate calculations of the energy shift. In the vicinity of E=E one has 


according to Eq. (5.3) 


The second derivative of [F1(§x'/§}.—F1'/F1)} is neglected. Minimizing the above expression, 
one obtains the condition 


p+ = 


G;? 
] 
kaE 


Similarly the corresponding approximate condition for the minimum of |Gz/§z|? 
is 


F 
E'-E,= 


pak Guar 


Here AE is introduced by means of Eq. (5.5). For high barriers | F,/Gz|<«1 and | E’—E,| AE. 
The shift of the peak in such a case is small in comparison with the resonance width and the barrier 
penetration factors enter in the same power in (E’—Eo)/AE as in AE/E. According to Eq. (5.8) 
E-E, is of the same order of magnitude as E’—Eo. The term (E»/F1)0(1/Gz)/dE contains the 
penetration factors in the same power (zero) as the other two terms in the numerator of Eq. (5.8) 
and there is no special general distinction between them. It is thus seen that the energy shifts indi- 
cated by Eq. (3.9) are the less important the smaller | Fz/Gz|. 
This conclusion can also be seen from Eq. (5.6). For E=Eo+AE the value of G,? is 


~ (Gxt) / f 


The latter estimate is very rough. The point is One-body problem with strong repulsion 
that the addition to G,? is of zero degree in G2. If the nucleus is surrounded by a strong re- 
The fractional change of the factor Gz? at  pulsive field acting in addition to the centrifugal 
E=Ey+AE is roughly 1/G,’ and becomes negli-  parrier! then it is more natural to modify the 
gible for high barriers. oe above procedure and to introduce wave functions 

The essential point of the approximation tO in the repulsive field. The repulsive field has to be 
resonance by the method of — lex eigenvalues continued to r=0 and there is some arbitrariness 
and the method of regular and irregular functions involved in this step. The procedure is illustrated 


is the same. Resonance occurs nearly at the . 
‘ " Pa in Fig. 1. The heavy line is the potential energy 
energy for which |1/§z|* has a minimum. At the curve. The light line is the continuation of the 


minimum this quantity can be approximated by 
a parabola and in this approximation formulas of 4 The repulsive field considered here is not the Coulomb 


: field. The calculations are easily modified to include this 
the resonance type are obtained for cross case also. The notation used here is convenient for appli- 


sections, cations to many dimensional problems. 
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Fic. 1. Potential well and potential barrier. Heavy line 
represents the actual potential. The light and dotted lines 
represent auxiliary repulsive potentials. The height of the 
straight horizontal line is the energy. 


repulsive potential to r=0. At r=a the two 
curves part. The region from r=a to r=d is the 
potential energy barrier and the region a<r<b’ 
is the region of negative kinetic energy. The 
auxiliary repulsive potential could be drawn in 
also as the dotted curve and in an infinite number 
of other ways. 

In the auxiliary repulsive potential one obtains 
a regular solution f, for r times the radial 
function. For r>6 this function is a linear combi- 
nation of corresponding functions without any 
potential. One has 


F, cos sin 
sin 6,+Gz_z cos 6;. (6) 


Here 6, is the phase shift of f, with respect to 
F, and the function gz is defined so as to be 
asymptotically the cos of the same argument of 
which f, is the sine. For r<6 the function gz is 
defined by the radial equation and continuity of 
the function and its derivative at r=). 

From Eq. (6) one obtains 


=F (6.1) 


The phase shift of § with respect to fz will be 
denoted by «xz and one has 


For r>06 one has according to Eq. (4.7) and 


(6.2) 


BREIT 


Eq. (6) 
= (fr cos sin xz) exp {i(xz+8,)}. (6.3) 


But from r=a to r=6 the actual potential and 
the auxiliary repulsive potential are the same. 
Therefore Eq. (6.3) holds also for all r>q and 
can be applied at =a where the attractive region 
begins. 

The usefulness of this model is seen in cases jn 
which the region r<a can be pictured as an 


- approximation to the interior of the nucleus 


somewhat as in the case of the Coulomb barrier. 
In such a case one is interested in | ¥,|? in this 
region and one can obtain an idea of these by 
studying ||? at r=a because in 0<r<a the 
shape of 1 does not depend critically on the 
energy. At r=a one has, on the other hand, 
\gx/fx|>>1 for high barriers and one ‘obtains, 
therefore, higher reaction yields for approxi- 
mately those energies for which sin x; = +1. The 
scattering cross section depends, on the other 
hand, on sin? K, and one has thus an energy 
shift between the peaks of the scattering cross 
section curve and the curve for |§}z|? which js 
due to 6;. This shift is of the order of magnitude 
of the width of the resonance curve multiplied by 
5,/m because the half-value width of the reso- 
nance curve corresponds to a change in Ky by 

It is natural to expect a shift between the 
energy of maximum penetration and the energy 
of maximum scattering. Suppose, for instance, 
that the maximum of |§z(kc)|? occurs at the 
same energy (Fig. 1) as the maximum of sin? K,, 
The maximum of |§ (ka) |? will not occur at the 
same energy because the barrier between 
c and a affects the wave function. The ratio 
| §1(ka)/§xr(ke)|* varies by large amounts in 
the vicinity of resonance, provided a and ¢ 
are sufficiently far apart and an energy shift 
results. 

Since the repulsive potential can be continued 
to the left of a in an infinite number of ways the 
method of the auxiliary repulsive potential is not 
unique. Even though the functions fz, gz are 
natural ones to use they are only auxiliary 
quantities, while §, is the actual wave function. 


W Vv Fic. 2. Illustrating the 
pe cusp at the intersection 
—w and 


anal 


The i 


516 

V 

This 
eval 
ener; 
appli 
A 
ther 
sugge 
(6.43' 
Greer 
Cor 
in the 
with b 
It w 
ponent 
numbe 
possibl 
change 


RESONANCES IN NUCLEAR REACTIONS 


1 fx! ] 
6.42 


analogously to Eqs. (4.6), (5.3). For the width of resonance (AZ)’ as measured by | (ka) |* one has 


| ar+| -f [Maximum of x(ka)|*]. (6.43) 


This formula is analogous to Eq. (5.5). The function § 1 and the other functions occurring in it are 
evaluated for the energy which gives maximum |§z|? inside the repulsive barrier, however, and this 
energy is not the same as the energy for which sin? K; is a maximum. Equation (6.43) should not be 
applied, therefore, to an estimate of the width of the scattering resonance in the general case. 

A superficial comparison of Eq. (6.43) with Eq. (5.5) appears to indicate that (AE)’=AE because 
the right-hand sides do not depend on the choice of the upper limit and Eq. (6.43) could equally well 


itten 

(AE)’ 0 kdE g.or r=b 


The integral in this formula has the same form as in the corresponding modification of (5.5) viz. 


E b G,2E@ 
| 
AE 0 kdE rah 


suggesting that AE=(AE)’. This is not the case because the energy which must be used in Eq. 
(6.43’) is for resonance at a while in Eq. (5.5’) the energy is for resonance at b. 


Two-DIMENSIONAL PROBLEMS 
Green’s function for separable two-dimensional differential equations 
Consider the differential equation 


») =| —+— - »)=0 (7) 
dx? 
in the domain 
—a<x<a, 
(7.1) 
with boundary conditions 
¥(a, y)=¥(—a, y) =0. (7.2) 


It will be required also that for y=-++ © and fixed x the wave function should either vanish ex- 
ponentially or else be asymptotic to a sum of a finite number of sine functions of multiples of y with a 
number of fixed but arbitrary phase constants. In order that the latter boundary condition be a 
possible one, it will be supposed that V(y) approaches a constant value faster than 1/y. This constant 
value can be made = 0 without loss of generality by changing «? so as to absorb the difference and this 
change will be supposed to be made. 

The differential equation 


(d?/dx?+d, — U(x) Jun(x) =0 (7.3) 
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with the boundary conditions (7.2) defines in —a<x <a an orthonormal set with eigenvalues a: The 
function U(x) will be supposed to have either no singularities or if it has singularities these wi] be 
supposed to be sufficiently innocent so as to allow the existence of the set \1<A2<A3--- witha mini. 
mum \, ¥ — ©. A particular solution of Eq. (7) is obtained by letting 


v(x, y) = uUn(x) Ve(y) (7.4) 
(7.5) 


[d?/dy*+k?— V(y) ]¥i(y) =0. (7.6) 


The dependence of & on n is not explicitly indicated here so as not to complicate the notation. There 
are two linearly independent solutions of Eq. (7.6). These will be 


we(y), 


respectively. The solution v;(y) will be chosen so as to be acceptable at y= + © and w(y) so as to be 
acceptable at y= — ~. Thus if at y=>+ © the solutions of Eq. (7.6) are of the exponential type, 
v(y) will be taken as the exponentially decreasing branch and similarly w,(y) will be arranged to be 
exponentially increasing with increasing y. This condition is schematically indicated in Fig. 2. If, on 
the other hand, k? >0 then one may use any pair of linearly independent solutions of (7.6) for v and w, 
These solutions will be restricted later so as to give either diverging or standing waves. In all cases 


the functions v(y), we(y) will be normalized so as to have 


ve! (y)we(y) = 1, 
v'(y) =dve(y)/dy. (19 


If Eq. (7.7) holds for one y then it holds also for all y as a consequence of Eq. (7.6) according to the 
theorem regarding the invariance of the Wronskian. A kernel will now be defined as follows 


K(x, & 1) = (y>n), 
K(x, & 2) = (y<n), (8) 


where it is understood that k varies with m in accordance with Eq. (7.6). It will now be shown that 


with 


and with 


f f K(x, &, »)e(é, n)dédn= p(x, y). (81) 


This equation will be first verified for the special form 
0) =un(E)f(n). (8.2) 


One has 
+00 +a” 
£ f f K(x, & = Ltun(x) f Mf Yoly, 


where 
M=d?/dy?+x?—d,— V(y) (83) 


(84) 


and 


Yi(y, 2) =v(y)we(n) (y>n), 
Yi(y, (y<n). 


+00 


f(n)dn=f(y). 


It remains to show that 
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This is the case according to the construction of Green's function for one-dimensional problems as 
explained in Courant-Hilbert.'® The special case in which »(y),w, (y) are the same is excluded for the 
present. The verification of Eq. (8.5) is indicated next for the sake of completeness. 


da? a? 


0) f(n)dn=— ve(y)wa(m) f(n)dn+— f ve(n)we(y)f(n)dn 
dy’ dy dy 


y d fr” 
f (y)we(n) f(n)dn+— f ve(n)we' (¥)f(n)dn (y)we(y) we’ 
dy dy 


+ f f 


Substituting this into Eq. (8.5), taking into account (7.7) and Mu(y) = Mw,(y) =0 one verifies Eq. 
(8.5). In order that the calculation in Eq. (8.6) be valid it is necessary that the integrals occurring in 
it converge. In applications p(x, y) will vanish except in a finite range of values of y so that no 
difficulty will be met in this connection. Having verified Eq. (8.1) for the special form (8.2) one sees 


also that Eq. (8.1) holds for 
p(x, ¥) = (8.7) 


If the sum is finite the above verification applies directly. If the sum in Eq. (8.7) is an infinite series an 
investigation of convergence is required. This has not been carried out in the general case. No 
difficulty is expected in this connection in view of the fact that in special cases Eq. (8.1) can be 
verified by other methods. 
According to Eq. (8) two different expressions for K are used depending on whether y>» or y <7. 
These approach the same value }om(y)we(y)un(x)un(£) as the line y= is approached in the two 
regions and K(x, y; £, ») is continuous on that line. The derivative normal to the line appears to be 


discontinuous being given by 
[aK (x, & 0)/Ay LOK (x, & —ve(n)we'(m) (8.8) 
On account of Eq. (7.7) the right side of the above formula becomes 
= 5(x—£). 
The discontinuity in the derivative occurs, therefore, only at the point 


(x, 9) (é, 


_ If one wishes to work with complex eigenfunctions the product u,(x)u,() in Eq. (8) can be replaced 
by 
Un(X)Un*(E) 
and the above argument remains essentially unchanged. It should be noted that Eq. (8.1) has been 
verified only for such forms of p which can be expanded as in Eq. (8.7). This is a restriction which 
must be watched especially at the boundaries of x since the boundary conditions for y have no direct 
significance for p. 


General plan for applications of Green’s function 


In proper units one can apply the wave equation (7) to the discussion of a quantum- 
mechanical system consisting of two particles with the following interpretation of quantities: 
«=total energy, \,=energy of x particle in state n, k?=energy of y particle, and u,(x) =wave 
function of x particle in state n. 


*R. Courant and D. Hilbert, Methoden der mathematischen Physik (Springer, Berlin, 2nd edition, 1931), Vol. I, p. 302. 
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In Fig. 3 a schematic picture of the two- 
dimensional condition is sketched. The wave 
function is confined to the strip of width 2a on 
the boundaries of which it is supposed to vanish. 
The wave equation (7) holds outside the shaded 
region A. Inside this region another and in 
general a nonseparable wave equation is sup- 
' posed to hold. The potential energy function 
V(y) is schematically represented as a barrier 
sticking out of the plane of the paper. In order 
not to complicate the figure U(x) was taken as 
zero and V(y) was made to vanish in the shaded 
region. Neither of these simplifications matters in 
the applicability of the method. 

Outside the shaded region the wave function 
will be represented as 


ve [KG 


X y) (9) 


and it will be supposed that 


p(x, y)=0 (outside A); Lyo=0. (9.1) 


In consequence of Eq. (8.1) one has 
L£y=0, (outside A). (9.2) 


The function yo will be used to represent the 
incident wave and ¥y—yo will give its modifi- 
cation due to the interaction potential in the 
region A. The reflection of the waves is thus 
pictured as due to a source density p placed 
within the interaction region. In the analogous 
electrostatic problem the source density is re- 
placed by a charge density. Since for k*?>0 the 
functions v;%, w, are not uniquely determined by 
the boundary conditions and by Eq. (7.7) one can 
arrange for y— yo to have the desired asymptotic 
forms (outgoing, incoming, standing waves) at 

Some qualitative information can be obtained 
from the form of Eq. (9) without special as- 
sumptions concerning the location of p(x, y). It 
should be pointed out, however, that in a number 
of cases it is possible to condense p(x, y) to a very 
narrow strip on the boundary of A. This is a well- 
known procedure in electrostatics and it corre- 
sponds to representing the field due to a col- 
lection of charges by an equivalent surface charge 


density over a surface enclosing the charges, The 
possibility of doing so depends on the existence of 
solutions of elliptic differential equations in 
closed regions with prescribed boundary values 
and may be seen as follows. A will be bordered 
bya thin strip of uniform thickness just outside of 
it and shown shaded in Fig. 4. The original prob- 
lem will be modified by requiring that the wave 
equation (7) holds everywhere except in the 
shaded area of the strip. The values of y on the 
boundary of A (inner boundary of strip) will be 
taken to be those given by Eq. (9). These values 
determine a new y inside A which will be called 
¢. ¢ is a solution of Eq. (7) and is continuoys 
with y just outside the shaded border when the 
border is thin. The normal derivatives dy/an, 
d¢/dn are not continuous in going across the 
border as in electrostatics. Leaving the thickness 
of the strip finite but small one can join ¢ to y by 
smooth curves across the strip and define a 
function 


V=y outside A, 
W=¢ inside A, (9.3) 
W=joining function in strip. 


Now form 
a(x, y) = £Y. (9.4) 


From the definition of ¥ it follows that o(x, y) 
has a value different from zero only in the shaded 


Barrier V 


Fic. 3. Illustrating the 
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model. The wave function 
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cal strip of width 2a. The 
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strip. Consider 
x(t, =Vot f f K(x, no(, (9.5) 


This function satisfies 


£x =0 outside strip, 
inside strip. (9.6) 


Comparing (9.6) with (9.4) one sees that x and ¥ 
satisfy the same inhomogeneous equation. Also 
y) and K determine the asymptotic behavior of 
x and y at y=+ and one can see that this 
behavior is the same for y and x. In fact the 
choice of u(y) in (8.4) determines completely 
whether for y=+ the mth state of the x 
particle is excited together with y particle moving 
towards y= +2 or y=—o. It is usually as- 
sumed that the specification of boundary con- 
ditions in the above sense determines the wave 
function. Adopting this view, it follows that 


x=¥, 
x=y outside strip. (9.7) 


9) + f K(x, 93 &, f R(x, & adds, 


dYo(x, OK (x, y; &, aK (x, &, 


Ons Ns. 
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According to the above one expects to be able to 
replace a surface distribution of sources by a 
distribution of sources on a contour in agreement 
with physical expectation. One has for an 
infinitely narrow strip 


f Ke, n)o(, n)ds, 


ds?=dt?+dn". 08) 


The integration is carried out along the outer 
boundary of A and o in (9.8) is the previous ¢ 
integrated through the strip thickness."® 

The value of ¥ on the boundary of A and the 
differential equation inside A determine y inside 
A and hence also dy/dn at the boundary. The 
connection between the values of ¥ and those of 
dy/dn on the boundary is linear, therefore, and 
the problem reduces itself to the solution of linear 
equations with the density function ¢ of Eq. (9.8) 
on the one-dimensional boundary of the region A. 

This connection can be shown more explicitly 
by introducing Green’s function for the interior 
of A for the differential equation obeyed by y 
inside this region. In terms of it one can write 


= f R(x, & nds. (10) 


The continuity of y and its normal derivative on 
the boundary gives 


(10.1) 


(10.2) 


the outward drawn normal being denoted by n, just outside of the boundary and the inward drawn 
normal by n_ just inside, as in Fig. 4. The above equations involve values of ¢ and é¢ over the contour. 
To emphasize this write s for (x, y), s’ for (¢, ») and 


a(s, s’)=K(x, y; & 


a(s, s’) = R(x, y; & 0), 


b(s, s’) =8K(x, y; n)/dn4, 5(s, s’) (x, y; 0) /dn_, 


Sols) =Yo(x, 
One has then 


fo(s)+ a(s, s’)a(s’)ds’ = f G(s, s’)e(s’)ds’, 


go(s) + f b(s, s’)o(s!)ds! = f B(s, 


* The possibility of singularities in o is not excluded. 


= OYo(x, y)/dn,. 
(10.3) 


(10.4) 
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The construction of K and K thus reduces the problem to the solution of two simultaneous integral 
equations in a one-dimensional domain in two unknowns g, é. 
The linear connection between y and dy/dn_ on the boundary is expressed in the above by 


y= f G(s, s’)@(s’)ds’, dy /dn_= f s’)a(s’)ds’. 


It can also be expressed in other ways. Thus, for example, for the interior of a circle, inside of which 
the potential is axially symmetric, it is convenient to use polar coordinates and a Fourier expansion 
according to the azimuthal angle. The values of ¥ on the boundary determine then the values of 
the radial functions by integrals involving ¥(s) linearly and d¥/dn_ is a sum of such integrals, The 


formulas are 
Qn 
Ws) == ds=add 
0 


ay/an_=X [Ra'(a)/Ra(a) f 
0 


Substituting into the last of these equations by means of ¥(s) =fo(s) + Sa(s, s’)o(s’)ds’ and dy/an_ 
=go(s) + S(s, s’)o(s’)ds’ one has a single integral equation in o(s). The radial function belonging 
to e*”* is denoted here by Ry. 


Degeneracy due to excitation 

The spectrum of eigenvalues of Eq. (7.6) may be entirely a continuum. It is then always possible 
to find both » and w, for every « and . If the spectrum of (7.6) has also a discrete part then it may 
happen that for certain values of « and =7 it is not possible to find both v,* and w+. This happens 
whenever k? is an eigenvalue of Eq. (7.6) and in this case the energy of the two-particle system jn 
the incident state yo is just equal to the sum of the energies of the x particle in the state 7i and the y 
particle in the state k. The same function »,+(y) is (to within a constant multiple) the only acceptable 
one both for y= + © and y= — ~, and it is impossible to construct a satisfactory Yi+(y, 7) by means 
of Eq. (8.4) on account of the failure of Eq. (7.7). This degeneracy of the kernel can be avoided bya 
slight shift of the total energy and no difficulty arises as long as the energy balance for simultaneous 
capture of the two particles is satisfied approximately rather than exactly. As this condition of 
resonance is approached the product of the normalizing factors in 1°, w,* must be increased in order 
to satisfy v,°’w,*—v,-w,*’=1 and Y;i*(y, 7). This does not mean, however, that one obtains an 
infinite probability of the formation of the state because a(s, s’) in Eq. (10.3) also becomes large. 

For k~k both »% and w are to within a constant factor ~»- and it will be shown next that 


Yi(y, 0) 


This relation can be obtained as follows. From the differential equations satisfied by v, we, v4 one 
obtains 


The asymptotic forms of the logarithmic derivatives of the functions are 
(+ ©)/ve(+ ©) /we(— ©) =(—h?)!, 
0) 20) = — ©) = (— 
and the functions vanish exponentially as follows 


©) =w,(— ©) ©) ©) =0. 
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From Eq. (11) one obtains on setting a=y, b= + © a formula for (y)/ve(y) — (y) Similarly 
forming an expression analogous to (11) with w:(y), u(y) one obtains a formula for wx’ (y)/we(y) 
—v,'(y)/v*(y). Subtracting the second formula from the first, one obtains 


we(y) 
~ In the limit of k=& the right side of this equation can be made equal to unity if one sets 

= W, (11.1) 


and normalizes 


f =1. (11.2) 


The two equations (11.1) can never be satisfied for all y. For fixed y, however, one can make |k—&| 
sufficiently small so that the shape of v;.(n), we(n) is practically the same as that of v*(n) for |n| <|y|. 
For fixed y and sufficiently small |k—&| one has then the limiting form 


Yi(y, 2) (11.3) 


A special problem will now be considered. It will be required that the wave function be zero over 
the area of a small circle surrounding the point (£, 7). The wave equation (7) will be supposed to 
hold outside the small circle. The interaction between the two particles is thus entirely inside the 
circle where it is infinitely great. This example will be considered in more detail later but at present, 
by means of Eq. (11.3) one can see the effect of resonance due to excitation. Close to resonance the 
values of K(x, y; &, ») over the circle are’’ 


Here > c’ represents the effect of other states in K and suffix C means: on circle. Denoting the incident 
wave by yo(x, y) one has the approximate solution 


This formula is approximate only. Its main inaccuracy lies in the assumption that }¢’ in the de- 
nominator has the same value for all points on the circle. This is the case to a first approximation, 
the logarithmic term in (x— £)?+(y—n)? being predominant in }\¢’ for small circles. The effect of a 
single term in the summation in Eq. (12.1) may include in it a wave extending to y= «. The above 
formula shows that for k~k these waves become small. An individual term behaves as 


Deo! 


and shows a maximum for k?=k*?—1u,*?(£)v,**(n)/Soc’. The latter equation cannot be satisfied accu- 
rately because >’ is only approximately real. The resonance has, therefore, a finite width. At the 
maximum the above calculation breaks down because then the parts of t,*(x)tuq*(£) Y:(y, 9) which 


have been neglected become important. For k=k, however, the above scattered wave is zero. This, 
of course, is also only an approximation. For k=& one has from (12.1) 


y) =Polx, ¥) —WolE, 2) J. (12.3) 


"It is natural that une and vx occur symmetrically in this case. 


1) ten (12.2) 
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This expression is an exact solution of the wave equation and it vanishes at (£, 7). It is, therefore, 
also small on the little circle. It is understandable that the scattering due to the interaction within 
the circle is small at this energy.'® The incident wave is in a sense “‘sucked in”’ to produce excitation, 


Extension to several dimensions 
An extension of Green’s me given by Eq. (8) may be made by generalizing the variable x tp 
a set of variables: x1, x2, x;. The functions u,(x1, x2, are supposed to satisfy 


— 


0x? Ox,” 


The \, and u, are supposed to form a complete and discrete set of eigenvalues and of eigenfunctions 
for this equation. The wave pam: outside the region of configuration space in which the particles 


interact is 
Ly (x, +++, =0 (13.1) 


and 
3? 3? 
4 — — U +++, x.) — (13.2) 


Defining k? as before and subjecting v;(y), w.(y) to Eqs. (7.6), (7.7) one defines the kernel 


=> We(Y)Un(X1, Xs) Un(E1, é,) (y<n). 


One verifies in the same way as before that 


In applying the above generalization one may think of x1, x2, -- + x, as describing the internal con- 
dition of the bombarded nucleus and of the coordinate y describing the incident particle. Green's 
function K is designed so as to represent coherently scattered waves and scattered waves with excita- 


tion of the bombarded nucleus by means of the source density p(x1, -+-, Xs, ¥)- 
The restriction of the operator £ to the form of Eq. (13.2) is not essential. The method still works 


for 
where N is the wave equation operator (multiple of Hamiltonian) for the bombarded nucleus 
Equation (13) should then be replaced by 
(N—An)Un(X1, Xs) =0. 
The operator N may contain exchange forces. The separability of N into parts corresponding to 


X1, ***, %, is not assumed and is not necessary for the working of the method. The degeneracy due to 
excitation works out very similarly in the many-dimensional and the two-dimensional models. 


Modification for central fields 
A schematic discussion of radial functions for particles colliding with a nucleus can be carried out 


along the above lines. The differential equation 


LY(ni, 1) = +x?— r)=0 (14) 
or or? 


18 Minima in scattering can be produced also in other ways as has been brought out by M. R. MacPhail, Phys. Rev 
57, 669 (1940) from the general point of view of interference. 
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will be considered subject to boundary conditions 
¥(0, r)= (ri, 0) =0. (14.1) 


The distances of particles 1, 2 from the center of gravity of the system are here denoted by ri, 7. The 
potential energies Vi, V are supposed to include the centrifugal force terms L(Z+1)/r* which enter 
if the angular momentum L+0. Particle 1 will be supposed to be confined to the nucleus and, as in 
the preceding section, the scheme may be generalized to include the case of a composite system 
instead of the single particle 1. 

Equation (14) corresponds to the two particles being in separate fields of force described by 
Vi(r:), V(r). By an extension of the previous scheme one can discuss the interaction between the 
particles provided the potential energy is taken to be a function of r; and r only. Actually the potential 
energy is also a function of the polar angles of the two particles and the present discussion is incom- 
plete in this respect. One may hope, however, that this idealization spoils only the quantitative side 
of the problem but does not affect radically the qualitative conclusions regarding effects of resonance. 


The eigenvalues A, of 
e Vi(r1) Jua(ri) =0 (14.2) 


are supposed to form a discrete set and the eigenfunctions u, are normalized by 


f | =1. 
0 


The number & is defined by 
(14.3) 
and functions f;, gz are defined by 


[d*/dp?+1—k-*V(p/k) =0, 
[d*/dp*+1—k-*V(p/k) =0. 


The function f; will be regular at r=0 and will be normalized so as to be asymptotic to a sine function 
of unit amplitude at large p. The function g, will be asymptotic to a cosine function and g,+if; will 
contain p only as e**. One has the relation 


=1, (14.5) 


the differentiations being performed with respect to p. If the functions f;, gy are not oscillatory 
(k? <0) then f; will still be understood to be regular at r= 0 and g; will be made to vanish exponentially 
at r= ©. In this case the normalizing constants in f;, g will be adjusted so as to satisfy Eq. (14.5). 
The exceptional case analogous to that already discussed under ‘‘Degeneracy due to excitation’’ may 
also arise and can be treated analogously to Eqs. (11)--- (12.3). 

The function g,(p) is analogous to »%(y) inasmuch as it is always admissible at + ©. The smaller of 
the two arguments of the kernel cannot occur in g; since g; is not admissible at r=0. One is thus led to 


—Ks(ri, 73 =D  (r>1’) 
(r<r’). 


The above kernel has the suffix S so as to indicate that its asymptotic behavior at r= © is that of a 
standing wave. The — sign present in Eq. (14.6) is due to the fact that substitutions vg, wf made 
in (7.7) would give —1 on the right side of Eq. (14.5). Practically the same calculation as that carried 
out in Eqs. (8.5), (8.6) gives 


of 1’) p(ri’, r’)dri'dr’ = p(n, 7). (14.7) 


(14.4) 


(14.6) 
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Another possible kernel is 


This kernel contains for large r only diverging waves and the index D refers to “diverging,” By 
adding to the incident plane wave solutions of the type 


f f Ko(n, 7; ri’, r’)p(ni’, r')dr;'dr’, 


one is sure that at a large distance the asymptotic form of the wave function is the incident plane 
wave plus a diverging wave. 


QUALITATIVE APPLICATIONS 


Large repulsive interaction in small region 


The interaction between the two particles will be supposed to be confined to a small circle ag jp 
Eqs. (12), ---(12.3). On the circumference of the circle the wave function is supposed to vanish g 
as to correspond to an infinite repulsion within the circle. The circle will be referred to as an “‘obsta. 
cle” because its interior is a region into which the wave cannot penetrate.’® 

In the absence of the obstacle the wave function will be taken to be of the form 


7) = Uno(1)feo(kor) exp ] (15) 
and for any real k the asymptotic form at r= ~ will be 
Si (kr) ~sin 
The functions f,(kr), gx(kr) are analogous to fz, gz used in the discussion of the one-body problem ip 
Eas. (6), - - (6.43). The function yYo(r1, 7) replaces the term 
tUno(?1) Fx (kor) o(r1, 7) (15.1) 


in the expansion of the plane wave 
Uno(r1) exp (tkoz) = Uno(71) +1) Pi(cos 6) Fi (kor)/ (kor). (15.2) 


The substitution of (15) by means of (15.1) into (15.2) gives a modified plane wave consisting of 
the plane wave uno(r1) exp (tkoz) plus a scattered wave. The initial state of the particle 1 is mp and 
the particle 2 is in the state exp (ikoz) and the substitution indicated by Eq. (15.1) gives, there. 
fore, the solution of the scattering problem, taking into account the central field potential V(r) 
but neglecting the interaction between particles 1, 2. To the function yo(r:, 7) one may add 
JSS Ko(ni, 7; 11’, r’)dri‘dr’ with p(r,’, r’) =0 outside the obstacle and the result may still be 
substituted in place of uno(r1) Fx(Ror) in formula (15.2) without adding anything to the wave function 
except diverging and exponentially decaying waves. If the obstacle is sufficiently small to make 
¥o(r1, 7) practically constant on the circumference, then one may use as an addition to yo a constant 
multiple times Kp(r1, 7; 710, 70) where (710, 79) is the center of the circle. In this approximation one has 


to replace 


tUno(?10) feo(Roro) exp [48(ko) (ri, 73 rio, 70) 


(15.3) 
(Kp(riota cos 8, ro+a sin 8; m0 


Vo(71, 7) —>Uno(?1) fro(Ror) exp [45(ko) ]— 


19 In the present model the interaction is large whenever one of the particles is on the sphere 7; = 719 and the other on the 
sphere r = 79. It does not correspond to a simple law of force on a two-particle model. The requirement of a vanishing wave 
function in the triangular region r; <r would puiage be more natural since it would correspond to the incident particle 
never penetrating closer to the center than particle 1. This type of interaction is of the type considered in Eqs. (16 

-+ (16.4). No such model can represent the actual case perfectly and extreme cases are, therefore, worked out. 
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Here a is the radius of the circle. The value of Kp on the circumference is averaged over the angle @ 
between the radius and the r; axis. For a small circle the value of (Kp)« depends logarithmically on a 
as for the electrostatic potential and its dependence on the energy is not especially critical. The 
second term in Eq. (15.3) contains a part in Uno(?1)geo(Ror) +4fx0(Ror)_] which combines itself with the 
divergent wave part of the first term to give the coherently scattered wave. The contribution due to 
this wave to the plane wave in Eq. (15.2) is 


Vs,coh = 


0 D/™ 


Distinction of potential scattering.—Here the specific nuclear scattering adds itself to the potential 
scattering represented by sin 6(ko). It should be noted that in this problem a definite distinction 
between pure potential scattering and potential scattering + nuclear scattering can be made. For 
pure potential scattering one has 


sin 5(ko)P1(cos 0)tno(r1) exp {tL Ror 


Vs, coh, pot > 


For a known y,con one obtains a coefficient sin 6(o) exp [48(Ro) ] and the point is that an arbitrary 
number cannot be represented as e® sin 6 with a real 5. This complex number has the special property 
of having a square of its absolute value equal to its imaginary part. 
For incoherent scattering with excitation m the above model gives 


L+1 P 1 (cos 6) 
un(ri) exp {iLp+4(ko) }}teno(i0) fe(Rro) feo(Roro). (15.5) 
hor k(Kp) 


2 
Vs, incoh = + 


Equation (15.4) gives the scattered wave for which the nucleus is left in the original state mo. Equation 
(15.5) gives the scattered wave for which the nucleus has been transferred into another (e.g. an 
excited) state. 

The denominator (Kp) contains a number of the f; and the barrier penetrabilities enter in compen- 
sating products fig, since all quantities in (Kp) refer to the inside of the nucleus. The dependence of 
(Kp) on energy is, therefore, much less marked than that in the numerator. One sees that primarily 
it is 

Un(rio)fe(kro)/k, 
that matters for the amplitude of the scattered wave in the present example and in this case |f;|? is 
a measure of the probability of escape to the state n. 


Repulsive interaction in large region 


Some of the qualitative features of the above example hold also in more general circumstances. 
Several small repulsive regions inside the nucleus give a dependence on energy of the scattered wave 
similar to that of Eqs. (15.4), (15.5). Sources of Kp put inside each obstacle have then strengths 
determined by (Kp)~! for each obstacle and the intensity of the incident wave + wave scattered by 
other sources. If the obstacles are sufficiently small and there are not too many of them the wave 
scattered by other sources can be neglected and the energy dependence of scattering or excitation is 
the same as before. 

If, however, the scattered wave inside the nucleus is comparable with yo(r:, 7) then the energy 
dependence of y,,coh, Ws, incoh Outside the nucleus is also changed. Equation (15.3) may be used for a 
crude estimate of some of the effects. The repulsion will be supposed to take place in a circular region 
having dimensions comparable with the wave-length inside the nucleus and it will be supposed high 
in this region. As a crude approximation the wave yWo(r1, 7) will be canceled approximately by a source 
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of Kp at the center so that Eq. (15.3) still will hold. Only in this case Kp does not depend Primarily 


on the radius of the circle and the denominator may have a decided energy dependence. 
If one of the terms in Eq. (14.8) for Kp is appreciably larger than the others inside the Nucleys 


and if for it ge:(k’ro)>fx(R’ro) then one may set 


—(Kp)~ tn: (rio) fe (R10) (R’r0)/R’. (16) 


Here x?=),+k” so that the state n’ goes with the wave number &’. If, besides, n’= mo, k’=ky if 
the incident state is the one that is important inside the nucleus then Eq. (15.4) becomes 


f 


2L+1 
Vs, coh ~———Uno (71) P i(cos 6) sin 5(Ro) 
k £x0(Roro) 


and from Eq. (15.5) one obtains 
2L+1 


Un(P10) Si(Rro) 


exp exp [i(kor+5(ko)} (16) 


P1(cos exp {i[p+6(ko)+4(k) ]}. (16.2) 


Vs, incoh“~ — 


Equation (16.1) corresponds to the scattering cross section 
Feoh~ (Ao?/m)(2L+1) | sin 5(Ro) _] exp ]|* 


and (16.2) to 


The last formula can be compared with the 
result of the considerations of Kapur and Peierls 
as used by Konopinski and Bethe” and those by 
Weisskopf and Ewing.”! Konopinski and Bethe 
apply the theory to Li’+H'—He‘+He‘* and 
they use 


« E4/G/, 


where G; according to their paper is the irregular 
function of the proton for the Coulomb field. 
Equation (16.4) contains the factor A| gio(Roro) 
which shows a similar variation with energy. 
The following points may be noted in a com- 
parison of the two results. 

(a) Equation (16.4) does not give a solution 
of the Li7+H'—He‘+He! problem in a direct 
way. It is obtained only as an approximate 
solution of the excitation problem. One may 
nevertheless expect some similarity between the 
energy. dependence of the excitation and dis- 
integration processes. Since the factors A|fi(kro) |* 
X|un(rio)|? and geo(Roro)|~* are 
characteristic, respectively, of the final and 
initial states one has reasonable grounds for ex- 


o~ AT; x (16.5) 


20E. J. Konopinski and H. A. Bethe, Phys. Rev. 54, 


130 (1938). 
21V. F, Weisskopf and D. H. Ewing, Phys. Rev. 57, 


472 (1940). 


P10) 


Fineoh ~ (AAo/1)(2L+1) | | | 


(16,3) 
(164) 


pecting the dependence Ao| uno(710) | ~*| geo(Rors)|* 
to hold also for disintegrations in the approxima. 
tion in which Eq. (16.4) has been obtained, 

(b) Since the formalism used is essentially the 
same for the case of several coordinates r;, --+,1, 
replacing the single coordinate r; the above cap 
be applied to the case of a nucleus excited by 
bombardment of the particle whose coordinate 
has been denoted by 7. If the nucleus disin. 
tegrates into two products (with a small probe. 
bility) from the state m a three-body disintegn.- 
tion is described. 

(c) The approximation used to obtain Eq, 
(16.4) is very crude. Essentially it neglects the 
effect of all states except uno(71)geo(Ror) /Ro which 
are produced by the sources ,of Green’s function 
inside the nucleus. The complete expression 
according to Eq. (14.8) contains other terms 
besides and these contribute to (Kp). In addition 
the approximation of replacing the many sourcs 
of diverging waves inside the interaction regia 
by one at the center neglects interference effect 
which contribute to an additional energy de 
pendence. 

(d) The function gio(or) which enters into the 
present discussion is not the irregular function 
for the Coulomb field but the irregular function 
for the solution of the radial wave equation wit 
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Fic. 5. The two-dimensional resonance model. 


the Coulomb field modified by the average 
nuclear potential taken into account by the 
difference between V(r) and the Coulomb poten- 
tial. (The Coulomb field is here thought of as 
breaking off to zero at a reasonable distance 
outside the nucleus. If this is not done the above 
formulas must be changed to take into account 
the effects characteristic of the 1/r potential.) 


Schematic representation of two-body resonance 


The interaction potential will be supposed to 
be large and attractive in a small circular region 
surrounding the point (rio, 79). Outside this 
region there is a narrow ring [Fig. 5] within 
which the potential is repulsive and sufficiently 
high to serve as a barrier for waves emerging 
from (rio, 70). The two-dimensional region formed 
by the circle and the ring is somewhat analogous 
to the Gamow, Gurney-Condon model of radio- 
active decay. The depth of the potential well and 
the height of the barrier will be adjusted so that 
for a suitable energy the wave function is large 
inside the circle and small just outside the ring. 
For this energy the circular resonance region 
will be found to act as a strong source of two- 
dimensional waves which will give rise to strong 
scattering. The wave function within the circular 
region is in this model a schematic representation 
of the compound nucleus which is formed by the 
coupling of the two particles to each other. In 
this model the compound nucleus is screened 
against external disintegration by the repulsive 
region in the ring and in addition also by the 
barrier V(r). 

The resonating properties of the circle-ring 
system will be first defined by coupling it to an 
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infinite plane with constant x. In plane polar 
coordinates the axially symmetric solution 
satisfies 


(17) 
which is the same as 
d?(Riy) 
X(Riy)=0. (17.1) 


Here V(R) is the potential within the circle-ring 
region and V(R)=0 for R>Rp». Outside the 
circle-ring system 


¥=Jo(«R) cos g—No(«R) sin g (17.2) 
is a solution of Eq. (17). One has 


Jo(xR) cos (kxR—7/4) 
(xR>1) (17.3) 


No(kR) ~(2/r«R)! sin (xR—72/4), 
so that by Eq. (17.2) 
cos (kR>>1). (17.4) 


Equations (17.4), (17.2) show that ¢ is caused 
by V(R) in the same way as the phase shift K is 
produced in the single-body problem. It will be 
assumed, therefore, that close to the resonance 


tan p= AE/(E—E)), (17.5) 


where AE is a positive constant. This assumption 
is reasonable in view of Eq. (17.1) which is like 
the radial equation for § in the central field 
problem. Eq. (17.2) does not apply for R< Ro. 
At R=R, it joins smoothly to the internal 
solution. One has 


(17.6) 


( dy ) _So'(«Ro) cos No'(«Ro) sin 
VdR/ Jo(xRo) cos y—No(xRo) sin 


and this formula does not depend on whether 
the ring is coupled to an infinite plane or to the 
actual potential provided one may assume y to 
depend only on R for R< Ro. In fact this axially 
symmetric part of y is determined by Eq. (17) 
and the requirement of regularity at R=0. The 
requirement of axial symmetry is good only if 


x= 
For small x 


Jo(x) =1—2°/4+-->; 
tNo(x)/2=Jo(x) In 
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The formula 
Jo(x).No' (x) — No(x) Jo’ (x) =2/xx 
when used to eliminate No’ in Eq. (17.6) gives 
(= 2/mRo 
¥dR/ Ro 2 cot g—(2/z) In (yx/ 2) 


Since x<1 this expression will be approximated 
by 


( ) —2AE/Ro 
xp E—Ey—(2AE/z) In (yx/2) 


where the approximation (17.5) has been used. 

In the presence of Vi(r:), V(r) and with the 
boundary conditions y(0,r)=y(r:,0)=0 the 
approximation 
¥(r1, 7) = Poli, 7) 

+ayol(rio, o)Kp(ri, 7; 710, 70) (18) 

will be used. Here yo is given by Eq. (15) and 
it is assumed as in Eqs. (15), --- (15.5) that the 
density p in R< Rp can be replaced by a point 
source. This assumption is a poor one close to 
R=R, but the axially symmetric part of y can 
be represented by a point source. One has 


dy /dR ao(r10, ro)(dK p/dR) Ro 


( dy (dK p/dR) Ro 
VdR/ Ro 


and 
(18.1) 


1/a+(Kp) Ro 


The scattered wave at a large distance is 
2L+1 
kor 


sin 8(ke) exp [ikor-+i8(bo) ] 


G. BREIT 


Sxo(Ror) exp [46(Ro) sin 5(ko) exp {iL Ror ]} + Fro(Ror) 
gu(kr) +if,(kr) exp {i[kr+4(k) ]}. 


1 
—attno(? 10) fro(Roro) (Rk? wo) x(Rro) exp [ikr +15(k) +15(ko) 


The first term in the braces represents scattering due to V(r). The remainder is due to the resonance 
interaction. Neglecting the former one obtains the scattering cross section due to resonance scattering 


k 
0 0 


The variation of Yo over the circle-ring 

is neglected and in (18.1) the average values of 
d¥/dR over R= Ry are used. These values will 
be denoted as (y), (dy/ydR) and one obtains 


(dy/ydR) 
— /ydR)+(dKp/dR) 


cot g— No 1 


One has approximately 

4iK. D™ J 0 
and hence 4(dKp/dR) ~KNo so that 
4 


(No) —cot 


4\E 


On account of the potential field V(r) the 
quantity tno(71)Fx1(Ror) in (15.2) should be re. 
placed by yYo(ri, 7) as in Eq. (15.1). On account 
of the resonance interaction in R< Rp one must 
make the substitution 


yoy 


with y given by (18), (18.3). These substitutions 
are made in the right side of (15.3) and 
Uno(r1) exp [tkoz] is subtracted from the result, 
The difference is the scattered wave. Its value 
for large r is obtained by means of the asymptotic 


expressions 


(18.2) 


(18.4) 


(18.5) 
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where ys is the part of y, due to excitation into the state » and emission of the wave gi+if:. The 
factor k/ko is introduced into (18.5) because the number of emitted particles is proportional to their 
velocity. Substituting into (18.5) one finds 


Ao? 1 1 
Ores = (2L+ x0*(Ro?10) (Rk? > (18.6) 


Introducing a by means of Eq. (18.3) one has 


QL+1) 
| 


(18.7) 


1 


The expression (14.8) for Kp is now used to obtain (Kp). The imaginary part comes only from 
(#>0) and consists of a finite number of terms. The real part is an infinite series and has a loga- 
rithmic singularity like that of No/4. Denoting the real part by Re and the coefficient of the imaginary 
part by Im one has from Eq. (14.8) 


(18.9) 
which when introduced into (18.7) gives 


(2L-+1)AcT 


(19) 


This is the ordinary formula for resonance competition. It has been obtained with the following 
restricting assumptions. 

(a) The two-dimensional region within which the interaction takes place is small. This assumption 
does not mean that within this region the interaction is small. The result of this assumption is to 
simplify the result by making I, depend only on u,?, f,? at r1=110, r=1o. (b) Resonance was assumed 
to be sharp. 

According to Eq. (18.9) the damping constants I’, depend on the constant 4AE which is charac- 
teristic of the resonance region. Besides, there is present also the dimensionless factor u,?(r10)fi2(Rro) /k 
which is characteristic of the binding of particle 1 in the residual nucleus, the energy of the escaping 
particle and the potential barrier through which this particle must escape. The barrier enters im- 
plicitly through f,?(Rro). The presence of this factor could be expected by analogy with the result of 
Wigner and the writer by a method analogous to the Weisskopf-Wigner theory of the absorption and 
emission of light. The damping constant is then proportional to the square of the matrix element of 
the interaction energy divided by the average energy interval, | M,|*/Av, in the notation of the above 
reference. This quantity varies with the energy as f;*/k in exact agreement with the result obtained 
on the present model. It is claimed, on the other hand, from the considerations of Kapur and Peierls 
that ' « £4/G,*. The difference between the results is in a sense quantitative rather than qualitative 
and it is perhaps unnecessary to emphasize it. For slow neutrons both results give I « E! and for 
small r one has F,G, «<p « E* so that F,2/k and E!/G, vary with energy in the same way. For low 
energies f,;? and F,? also vary similarly. On the other hand f,2/k, E/G,? do not vary with energy in 
exactly the same way and it should be noted that f;,2(ro) is here taken at r=7 i.e., at the mean r of 
the compound state rather than at the arbitrary nuclear radius of Kapur-Peierls. 
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SPECIAL Two-DIMENSIONAL GREEN’S FUNCTIONS 


(A) Infinite plane 
The differential equation 


is supposed to hold through an infinite plane. Green’s function is 
Kp= (1/41) Ho ([(x— (20) 
where Ho” is the Hankel cylindrical function of the first kind. The method of Eqs. (7) will now be 
applied to the computation of Kp and the result will be shown to agree with Eq. (20). The point of 


doing so is to make sure that the method works also for a continuous spectrum of X,. 
Quantizing in the interval —L/2<x<L/2 by means of the periodic boundary condition one has 


tt, (x) = 
t=2nn/L; n=0, +1, 


and 


Un(x)Un*() exp [ir(x—&)]. (20.1) 


Equation (7.6) gives 


ve(y) = wily) = 


with 


k=(1-—7°)! for 1>77 
k=i|(1—7?)!| for 1<7? (20.2) 


and according to the normalizing condition (7.7) one has 
kN (20.3) 


Substituting (20.1), (20.2), (20.3) into Eq. (8) and remembering that the interval between successive 


t is 2r/L one has for large L 
1 

K (x, & gt! 


According to Eq. (20) one has then 


exp (y>). (204) 


(1—7?)! 


=- exp y| ]Jdr+— 

e (i- mid; 
This expression for H is not obviously symmetric in x, y. It shows the importance of states with 
negative kinetic energy for motion in the y direction. The /{% part of the formula consists of such 
states. It should be noted that these states of negative kinetic energy may contribute to the asymptotic 
value of Hy for x. Setting x=0 in Eq. (20.5) one obtains 


2 1 
Ho (y)=— sinh di+— >0), 20.6) 
(y) a) f e (y>0) ( 


#0 


cos Tx 


exp [—(7?—1)?|y| Jdr. (20.5) 


which is obtainable from Eq. (2) in Watson’s Theory of Bessel Functions, p. 178.” Setting y=0 one has 


2 ¢' cos rx 2 ® cos Tx 
f d (20.7 


Ho (x) =— 


2G. N. Watson, Theory of Bessel Functions (Cambridge University Press, 1922). 
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which is verified by means of 


i 2 
Jo(x) =— f cos (x sin 6)d@, No(x)=—— f cos (x cosh f)dt, Ho =Jo+iNo. 
0 


For large x, Jo(x), No(x) have comparable values and Eq. (20.7) shows that No arises in this case 
from waves with a negative kinetic energy of the y motion. The most important values of r are those 


slightly >1 which correspond to only slightly negative kinetic energies of the y motion. In Eq. (20.6), 
on the other hand, the negative kinetic energies contribute the first integral and one has 


2 1 
f e~tul sinh tdp< f tdi = —_., 
0 0 ly| 


This is negligible in comparison with Ho (y) for large |y|. Along the y axis, therefore, the negative 
kinetic energies of the y motion do not contribute to the dominating term of the outgoing wave. 


(B) Quadrant x>0, y>0 

The same differential equation as in case (a) is being solved subject to boundary conditions (0, y) 
= (x, 0) =0. Green’s function can be obtained again by the general method and this time by the 
modification for central fields, as in Eq. (14.8). The substitution r;—x, r—y is made, the functions u, 
are quantized in the finite interval 0<x<L and are (2/L)! sin rx while the interval between successive 
ris Ar=2/L. One obtains 


sin [(1—7*)!y] exp ]dr (21) 


2 f sin rx sin ré 
0 


(1-7)! 
with the understanding that for 
(21.1) 


Rearranging expression (21) one expresses it in terms of integrals of the form (20.5) and obtains 


1 
§)?+(y¥— n)*}!) — Hy n)* }') 
— Hy (((x— §)?+(y +0)? + Ho (21.2) 


which can be also obtained from (20) by placing negative images of the source at (—&, »), (, —n) 
and a positive image at (—é, —7). 

It will be noted that in both of the above examples Kp is asymptotically an outgoing wave at a 
large distance even though this is not obvious in the forms (20.6), (21). 


(C) Straight channel 
The region considered is —-a<x<a, 0<y< @ and ¥(x, y) is subject to the boundary conditions 
v(a, y) ¥(—a, y) =0 


and for y= it is required that Kp be an outgoing wave. These are the boundary conditions required 
for Eq. (14.8) for use with central fields. The substitution 


(ri, 7) =(x-+a, y) 
is made. The wave equation in the strip is taken to be 


(0°/dx?+0?/dy*+ =0. 
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One has 
Un(x) sin [wn(x+a)/2a]; n=1, 2,3,--- 
(22) 
(#>0)— [cs cos sin xn: 
2a|k| A 2a 2a 
For x0 this expression becomes 
1texp {(x/2a)[—(y—n) +8)]} 1—exp {(w/2a)[—(y +n) 
1—exp {(/2a)[—(y—n) +i(x—8)]} 1+exp 
(y>m). (221) 
The expression under the logarithm is a function of x+7y and, therefore, the right side is a solution 


of Laplace’s equation as it should be. The boundary conditions are verified directly in Eq. (22, 1), 
Also one sees that Kp and its derivatives are continuous except at x+iy=£-+7n. At the latter point 


Kp behaves as 


1 
—Kp=—R.P. In 
2r 


2a 


Constant terms have been dropped. Except for constant terms Kp is seen to be the electrostatic 
potential due to an electric charge of linear density (— 1/47) on an infinitely long filament perpen. 
dicular to the plane (x, y) and crossing it at the point (&, 7). This is in agreement with Poisson's 
equation in electrostatics which is in two dimensions 


(=+ —p(€, 0) In 1)? = 9), 


which is the same as Eq. (14.7) for Kp in the limiting case of Eq. (22.2). 

If x >0 but <zx?/4a? one may still apply Eq. (22.1) in the rectangle 0 <y < 2a as an approximation, 
Since in the form of Eq. (22) the series converges slowly one sees here the importance of taking into 
account many states of excitation of the residual nucleus. 

For any « the series (22) for Kp behaves for high like the series (22.1) and as long as | y—»| #0 
this series remains convergent when differentiated with respect to y. The convergence is due to the 
presence of the factor exp [—n|y—7| /2a]. It will now be verified that 0Kp/dy is continuous on the 


line For a small ¢>0 one has for x* <2?/4a? 


1 
(OK p/dY) (8K OY) L —] COS £)—(—)" cos sinh 


The second term in square brackets gives an absolutely converging series which vanishes with ¢ and 
may, therefore, be disregarded. The first term also gives an absolutely converging series if ¢>0. 
For e=0 the series does not converge. As -—- the sum can be seen to —0. In fact ¢ can be made 9 


small that erN/2a<1 for N>2ax/x. The xi is split into E+E and the series is compared with its 
N 

value for the same e and x=0. The difference of the > is negligible on account of the assumption 
1 


about ¢ and the difference of > is negligible on account of the assumption about N. On the other 
N 
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hand for x=0 the series can be summed and gives 0. It is thus seen that 0Kp/dy is continuous on the 
line y=. The point (x, y) =(&, ») is, of course, an exception. 

The main object of the above example is to illustrate the importance of excited states in the 
vicinity of the region in which interaction takes place. This is seen in the limiting form of Eq. (22.1) 
in which the outgoing waves combine to give the electrostatic solution. If there are two obstacles, 
instead of one, and if their distance apart is smaller than the wave-length the wave scattered by 
(&:, m) gives a contribution to y at (£2, 72) which contains many states of excitation. 


(D) Transformation of Hankel function 


The relations of Eq. (20) suggest formula (20.6) for the Hankel function Ho. This will now be 
verified by means of contour integration. According to Courant-Hilbert.* 


1 
Hy (r) = f (23) 


where the path of integration is from —i~ to 0, then to —7 and finally to —++i~ as indicated in 
Fig. 6. The differential equation for Ho“ is satisfied also for other paths of integration provided 


[cos ¢ sin (23.1) 


where {=a and {= are the limits of integration. The path will be deformed into L;’ which starts at 
and goes over fo, to with [This path is not the path 
denoted by L;’ on p. 408 of Courant-Hilbert I.] The condition (23.1) is satisfied by Li’ because 
cos £)>1. The integrand of (23) vanishes at the ends of the path for both LZ; and L,’ and no poles of 
the integrand have been crossed in the deformation. One has, therefore, 


1 
Hy(r) = —- f sin (23.2) 
Ly’ 


Setting {= and introducing the abbreviations 


y=r cos fo, sin (23.3) 
one obtains 


2 1 
Hy (r) =— f cos (x cosh n)e~¥ "dn+— f sin 609, (23.4) 
Tt 


T 


Substitution of r=cosh 7 in the first integral changes it into the second integral of Eq. (20.5). The 
second integral in (23.4) is transformed as follows 


1 1 2 
- exp (ir sin 0)d@=— f exp (iy sin @—ix cos @)d@=— f cos (x sin 
0 0 


The second integral in (23.4) is thus the first integral in (20.5) so that the latter representation is 
equal to (23). It will be noted that the restrictions 


0<in<x/2, x>0, 


are essential for the deformation of L; into L;’ and that these conditions are also essential for the 
applicability of the construction of Ho as 4iKp. 

The somewhat unexpected features of (20.5) are its dissymmetry between x and y and the fact 
that it is not obvious from its form that it represents outgoing waves in any direction but that of 


* Reference 15, Vol. I, p. 407. 
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increasing y. The way in which (20.5) gives an outgoing wave for y=0 has been discussed in connes, 
tion with Eq. (20.7). It may be pointed out that 


2 sinrx 2 ¢' sin rx 2 sinrx 


is also an outgoing wave for x >0. One has in fact 
2 sin rx 2\! 
f —dr= —iJ (x) ~ -+(—) cos (x-") 
mid, (7r?—1)! mx 4 


and 


(1—7?)! 


Combining these relations one has 


2 f® sin rx 2\3 
T (1 TX 


This formula is another example of an outgoing wave which is compounded of stationary waye, 
It and Eq. (21) show that the general method of Eq. (14.8) when applied to a continuum can give jp 
special cases outgoing waves in all directions within the quadrant x>0, y>0. This does not mean, 


however, that the same is true in the general case. 


FINAL REMARKS parison with the difference between the actul 
The calculations of the last section show that and inchtent sine functions ineide the nucleus 
the behavior of Kp close to the singular point is = many-Gimencional Gagram. The importance 
in fact such as has been used in the calculation of the — ediate state is = this eee 
of scattering by a small obstacle in the (rn, r) pa aie the = ory of interacting Particles 
plane and the related calculations on the inter- model treated 
action in a large region as well as with resonance. above gives the expected anewer. It Galle 
The example of the straight channel shows how 
at the singularity the logarithmic dependence on a 
(x—£)?+(y—7)*? is obtained. In fact for small 
values of this quantity the contributions of high 


1 


terms in the series of Eq. (22) approach the A oA 

limit for x=0 and give the electrostatic term. 

It is clear, therefore, that close to (¢, 7) one needs 

a larger number of excitation states to represent cE Pe <E 

the wave function. Ry 
-T 0 


In the two-dimensional model the singularity 
is logarithmic. In m dimensions the dependence 
of K close to the singularity is of type 
const. [>> (x;— &;)* The effect of a source 
of a diverging wave is the more localized in the 
m space the higher the number of dimensions 
and the region in which the higher states are ’ 
especially important becomes more confined as m L LL 
increases. The scattered wave function outside 
the nucleus may be expected to be small in com- 


> 


Fic. 6. Contours of integration for Ho. 
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include the effect of competitive emission of 
several particles but only the effect of competi- 
tion between various types of emission of the 
same particle. For the discussion of more general 
competitive processes it is obviously more suit- 
able to bring in the intermediate state in a more 
pronounced way. It is possible to do so by a 
generalization of the method of complex eigen- 
values to many-dimensional problems. This will 
be treated in another publication. 

The occurrence of the regular and irregular 
solutions of the radial equation (f; and gx) in the 
formulas for the scattered waves is not special 
to the two-dimensional model as is clear from 
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the many-dimensional extension in Eqs. (13.1) 
+++ (13.4). The examples show that in special 
cases the factor f;,” occurs in the formulas. If the 
interaction is distributed within the nucleus, 
however, then the general form of Green's func- 
tion shows that the g; occur in the denominators. 
In the general case the g, for all excited states 
will occur and very special assumptions about 
the relative importance of excited states are 
necessary to leave the simple factor [gio(Roro) |? 
in Eq. (16.4). Quantitative agreement with 
experiment obtained by application of formulas 
with one or another simple choice of such factors 
does not appear therefore to be very significant. 
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A study of the M-series lines of tungsten and of the K-absorption edges of magnesium and 
aluminum has been made by means of a double crystal x-ray spectrometer. The extension of the 
double crystal technique to these long x-ray wave-lengths has been made possible by thorough 
evacuation of the spectrograph and a windowless x-ray tube. The crystals were the aquamarine 
form of beryl in which the 1010 planes have a grating constant of 8.06A. The detector was a 
Geiger-Miiller counter with an aluminum window of 2.5X10~* cm thickness. The width and 
relative intensities of the prominent M-series tungsten lines were measured. Accompanying the 
diagram lines are some 33 satellite lines. In the study of the K-absorption edges of magnesium 
and aluminum, structure not heretofore resolved photographically has been recorded. A quali- 
tative comparison of the intensity of the continuous radiation from tungsten to that from 
aluminum at 10A indicates that the ratio is less than one, whereas theory predicts about six. 


INTRODUCTION 


RECISION measurements with the double 

crystal x-ray spectrometer have shown the 
effectiveness of the x-ray method in checking 
experimentally the theoretical calculations on the 
electronic conduction bands of solids. Ruled 
grating methods for the longer wave-lengths have 
given valuable information on the conduction 


'W. W. Beeman and H. Friedman, Phys. Rev. 56, 392 
(1939); H, Friedman and J. A. Bearden, Phys. Rev. 57, 
1085A (1940); T. M. Snyder and J. A. Bearden, Phys. 
Rev. 57, 1085A (1940); W. W. Beeman and J. A. Bearden, 


Phys. Rev. 57, 1085A (1940). 


bands of the light elements.? Photographic meas- 
urements in the intermediate range have not 
yielded results of definitive accuracy because of 
low resolving power and inaccuracies of intensity 
measurements. This region, however, is one of 
considerable theoretical interest, because here 
the present theory of metals makes sufficiently 
accurate predictions that their comparison with 
experiment may be expected to be significant.* 


* See, for example, Reports on Progress in Physics, Vol. 5 
(Physical Society, London, 1939). 

3M. F. ex hy H. M. Krutter, Phys. Rev. 51, 
761 (1937); J. C. Slater, Phys. Rev. 45, 794 (1934). 
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Fic. 1. Detailed construction of the spectrometer. The 
locations of the first and second crystals are shown. The 
walls J of the x-ray tube and the target T are of solid 
copper, kept cool 4 water. The target, which rotates 
about the first + ges may be removed from the instru- 
ment through the connection EZ. A baffle connected 
electrically to the high voltage prevents electrons reflected 
from the target from passing into the vacuum chamber. 
A ent arm indicates the setting of the first 
crystal on the circle A graduated in degrees. The position 
of the second crystal is read by a pointer on the circle B 
which also serves to locate the counter C. The counter 
rotates about the second crystal by the slip-ring pivot H. 
The arm D leads to the lever and micrometer head for 
slow motion of the second crystal. 


A vacuum double crystal spectrometer capable 
of making accurate measurements in a part of 
this intermediate region from 5A to 15A has been 
constructed and used to measure the width and 
relative intensities of the M-series lines of 
tungsten and the K-absorption edges of mag- 
nesium and aluminum. 


EXPERIMENTAL 


The double crystal spectrometer was assembled 
on a circular steel bed-plate as shown in Fig. 1. 
A steel bell jar permitted a vacuum of 10-* mm 
to be obtained. The x-ray tube was placed inside 
the vacuum chamber, eliminating the window 
which would otherwise be necessary. 

The crystal mountings were set in fixed posi- 
tions on the bed-plate and rotated on pivot 
bearings.‘ The second crystal could be rotated 


‘J. A. Bearden, Phys. Rev. 56, 1023 (1939). 


from outside the chamber by a micrometer head 
the motion being introduced through the bed. 
plate with the aid of a Sylphon bellows and lever 
A large dial containing 200 divisions was a ' 
to the micrometer head, making it possible tp 
read the position of the second crystal within two 
seconds of arc, or 0.07 x.u. at 7A. 

The detector was a Geiger-Miiller counter 
containing air at a pressure of 7 cm of m 
Due to the extremely thin window (aluminum 
foil 2.5X10—~ cm thick), a means was Provided 
for lowering and raising the pressure in the 
counter simultaneously with that in the chamber. 
The Neher-Harper circuit’ with a scale-of-eight 
thyratron circuit® and a Cenco counter were used 
for recording the counter pulses. 

The crystals used were the aquamarine form 
of beryl (BesAl2:SicOis), cut from a single speci. 
men along the direction of the 1010 planes which 
have an x-ray grating constant of 8.06A.’ The 
parallel rocking curves were recorded for three 
different wave-lengths and the full widths at 
half-maximum in seconds and x.u. are given in 
Table I. 

The percent reflection* was measured at 8559 
x.u. with the result P(0)=3 percent. This low 
percent reflection is due to the low atomic 
population of the 1010 planes.’® 

The physical resolving power of the crystals” 
from the measured (1, —1) widths is 1900 at 
7000 x.u. and 1700 at 8559 x.u. 


RESULTS 


M-emission lines of tungsten 


In recording the M-series lines, the target used 
was a thin disk of tungsten, spot-welded to 


TABLE I. (1,— 1) widths of crystals. 


X.U. w (SEC.) w (x.v.) 


135 4.8 
159 5.15 
199 6.5 


== 


5H. V. Neher and W. H. Pickering, Phys. Rev. 53, 
316 (1938). 

* J. Giarratana, Rev. Sci. Inst. 8, 390 (1937). 

7 A, Stenstrom, Ann. d. roe 57, 347 (1918). 

A. Compton and S. K. Allison, in Theory 
and Experiment (D. Van Nostrand, 1935), p. 722. 

*W. L. Bragg and J. West, Proc. Roy. Soc. Alll, 691 


(1926). 
10S. K. Allison, Phys. Rev. 38, 203 (1931). 
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nickel, which was in turn soldered to the base 
of the x-ray tube. The curves were all recorded 
with a current of 10 ma and at a potential of 
4800 volts; each curve was repeated a sufficient 
number of times to reduce the probable error 
to a satisfactory value. A sample of the data is 
given in Fig. 2 which shows the tungsten M6 line 
and the accompanying satellites. The other lines 
were recorded with comparable accuracy. The 
widths of the observed curves could be read to 
about 2 percent. In the absence of the data 
necessary for the development of a correction 
formula" for these crystals, the true width 
values have been obtained by decreasing the 
observed width by an amount equal to the full 
width at half-maximum of the (1, —1) curve for 
the corresponding wave-length. Table II gives a 
complete list of all the lines and satellites studied 
in the present work. 

The width of a diagram line varies from 0.31 
volt for the a2 line to 12.5 volts for the £; line.” 
This tremendous variation in width probably 
arises mostly from unresolved lines. For example, 
the line Mn which arises from the transition 
MyyOu,m in which the Oy,m levels are so close 
as to be indistinguishable, is doubtless an un- 
resolved doublet. 

Certainly the abnormal width of the M¢, 
cannot be explained on the basis of the sum of 
the widths of the levels from which it comes, 
since both these levels are involved in the transi- 
tions of Mé and Ma, which have the widths 
0.84v and 1.68v, respectively. 

A check on the relative intensities of Ma; and 
MB can be made on the assumption that they fit 
the equation” 


=a/(1—(x/d)?), 


in which a is the maximum intensity, } is the 
half-width at half-maximum, and x is the 
abscissa coordinate in appropriate units. The 
area under this curve is rab, from which it is a 
simple matter to calculate the relative intensities, 
defined for an instrument of high resolving power 
as the ratio of the areas under the curves. 


“L. G, Parratt, Rev. Sci. Inst. 6, 387 (1935). 
= ui ae and T. Magnusson, Zeits. f. Physik 88, 
# A. Hoyt, Phys. Rev. 40, 477 (1932). 
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Fic. 2. Intensity curve of Mé, showing two satellites. 
It is possible that 5” is more complex than is apparent, 
as suggested by the large width. The accuracy of the data 
in this curve is typical of all those reported in this paper. 


Experimentally this gives 
Ma:/MB=20.8/14; 


whereas the theoretical ratio is 20 : 14. 

The satellites were obtained by drawing in 
arbitrarily the background due to the parent line. 
Special methods were used in some cases, for 
example Mf, was drawn in as logical a way as 
possible, taking into account the fact that theory 
predicts the relative intensities of the three 
lines Mg,, Mé and M¢; to be in the ratio9 : 5 : 1. 

The satellite contours were resolved into satel- 
lite lines as shown by the dashed lines in Figs. 3, 
4, 5 and 6. The individual satellites. are of 
approximately classical shape and so chosen that 
they add to give the recorded curve. While the 
error may in some cases be very large, this is 
nevertheless the most satisfactory and accurate 
method available. The riumber of satellites 
appearing on some of the curves is surprisingly 
large. Three of the Ma, and three of the MB 
satellites have been photographed by Hirsch ;“ 
all others appear here for the first time. 


The K-absorption edges of Mg and Al 


In view of the necessity of a good source of 
continuous radiation when studying absorption 
phenomena, a test was made to compare the 


“F.R. Hirsch, Phys. Rev. 38, 914 (1931). 
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TABLE II. Wave-lengths of peaks and full width at half-maximum of emission lines of the M series of W(74) 
WavE-LENGTH SEPARATION CorREC | 
LINE TRANSITION X.U. X.U. SEc. x.U. — Ons, 3 
as My Nyt 6978 33.8 1.2 0.31 676 
My (6969) 0 186.0 6.6 1.68 6967 
a)! Sat. 6962 7 
a? Sat. 6960 9 50.7 1.8 0.46 
as Sat. 6956 13 50.7 1.8 0.46 6956 H 
a! Sat. 6947 22 152.0 5.4 1.37 6947 H 
a> Sat. 6939 30 101.0 3.6 0.92 
a8 Sat. 6931 38 169.0 6.0 1.53 6933 H 
B My Ny1 (6743) 0 220.0 7.85 2.14 6941 H 
B' Sat. 6735 6 67.6 2.42 0.62 6734 H 
B* Sat. 6729 12 101.0 3.64 1.00 
BS Sat. 6723 18 67.6 2.42 0.62 
Bt Sat. 6717 24 50.7 1.83 0.49 6719 H 
ps Sat. 6712 29 60.8 2.18 0.59 
Bs Sat. 6705 36 152.0 5.42 1.42 6707 H = 
Nun (8559) 0 4.97 0.84 
5 Sat. 8544 15 
e Sat. 8503 56 840 28.1 2.75 
o My Nin (8943) 0 2470 80.5 12.5 
ti Sat. 8867 76 335 10.9 1.69 z 
Sat. 8865 78 
Sat. 8858 85 
if Sat. 8851 92 44 2.32 0.36 zZ 
we Sat. 8842 101 
8 Sat. 8837 106 
7 Sat. 8828 115 173 5.6 0.87 
m8 Sat. 8812 131 258 8.4 1.3 
o° Sat. 8802 141 105 3.42 0.53 
ed Sat. 8795 148 54 1.76 0.27 
of Sat. 8783 160 
oi? Sat. 8775 168 460 15.0 2.32 
$2 My Nu 8972 0 410 13.3 2.3 8977 L 
$2! Sat. 8953 24 274 8.8 1.37 
2? Sat. 8934 43 54 1.76 1.27 
Sat. 8928 49 
$e Sat. 8920 57 105 3.42 0.53 669 
Sat. 8913 oF 
$28 Sat. 8909 68 
$27 Sat. 8905 72 
n Myv0n..101 6794 304 10.9 3.01 6794 L ine 
L: Lindberg; H: Hirsch; ( ) in col. 3 from Lindberg, used as reference points. the col 
alumin 
: intensity of the general radiation from a heavy W in this wave-length region. The efficiency total ¢ 
absorpt 


relation predicts, at this potential, e=0.029 per- 
cent for W, and 0.005 percent for Al. Further 
investigation in this wave-length region is needed. 

The present study of the K-absorption edges 
of Mg (12) and Al (13) was made using the con- 
tinuous radiation from a tungsten target, the 


element with that from a light one. The elements 
tested were W (74) and Al (13), and radiation of 
10A was used, since it was far removed from the 
characteristic radiation of either element. Under 
identical operating conditions (4500 volts, 8 ma) 
the counting rate with the Al target, over a 


INTENSITY 


twenty-minute period, showed a 30 percent 
increase over that for the tungsten target, the 
recorded intensity being due to both scattered 
and reflected radiation. About 50 percent of the 
measured intensity was due to scattering, which 
should be proportional to the intensity of the 
radiation incident on the first crystal. Since a 30 
percent increase is not likely to be due to the 
reflected 10A radiation alone, there must be an 
increase in x-ray production for Al over that for 


intensity of which was constant over the ranges 
used, which is not the case for an aluminum 
target. Although the aluminum in the beryl 
crystals might be expected to affect the re 
flectivity of the crystals in the region of the 
absorption edge, the percentage of the aluminum 
is low and probably affects the reflectivity in- 
appreciably. For the aluminum, the counter 
window served as the absorber with a thickness 
of 2.54X10-* cm. In the case of magnesium, at 
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absorbing foil 3X 10~* cm thick was hung before 
the counter window. The transmission curve for 
aluminum was repeated eight times, giving a 
total of 2500 counts per point on the high 
absorption side and 5000 on the low side. The 
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probable error is indicated on some of the points 
on the curve (Fig. 8). The magnesium curve was 
repeated twelve times, giving a probable error 
slightly less than that for aluminum. The curves 
have been plotted with mass absorption coeffi- 
cient against wave-length to eliminate the effect 
of foil thickness.'* 

By application of the quantum-mechanical 
radiation theory of Weisskopf and Wigner to 
absorption phenomena, Richtmyer, Barnes and 
Ramberg"* arrived at an expression for the shape 
of an x-ray absorption limit in the case where 
the absorption is due to the ejection of electrons 
from inner orbits to the continuum of unoccupied 
Fermi-Sommerfeld levels. The expression is in 
the form of a simple arctangent curve in which 
the horizontal separation of the two points, where 
the curve reaches } and ? of its asymptotic 
ordinate for high frequencies, is the width of the 
limit. Applying this relation to the observed 
contours of the Mg (Fig. 7) and Al (Fig. 8) curves 
yields the results shown in Table III for the 
width of the K states of the two atoms. 

The Kronig theory applies strictly to the 
structure 100 volts or more from the edge, where 
the energy of the electron is great enough that 
it may be considered practically free in the 
lattice. The present curves for Mg and Al 
indicate that the effect of the crystal structure 
is appreciable to within a short distance of the 
edge. Because of the similarity between cubic 


%F. K. Richtmyer and S. W. Barnes, Phys. Rev. 45, 


754A (1934). 
®F. K. Richtmyer, S. W. Barnes and E. G. Ramberg, 
Phys. Rev. 46, 843 (1934). 
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close-packed lattices and close-packed hexagonal 
lattices as regards the shells of nearest and next 
nearest neighbors, there should be a correspond- 
ence between the maxima and minima of the 
absorption edges in the two types of lattices. 
This has been pointed out by Coster,!” and the 
present curves bear out this fact. The distance to 
nearest neighbors in the Al lattice is 2.87A. On 
reducing the Mg lattice to an ideal close-packed 
hexagonal lattice, the distance to nearest neigh- 
bors is 3.18A. The distances of the maxima and 
minima of the Mg curve from the edge, multi- 
plied by (3.18/2.87),2 may be compared with 
those of the Al curve, the results being given in 
Table IV. The agreement is seen to be within 
experimental error. The first few maxima and 
minima do not check as well as might be ex- 
pected, since within about 15 volts of the edge 
the atomic fields predominate over the crystal 
potentials. The Mg curve matches each of the 
fluctuations in the Al curve and in addition con- 


TABLE III. Width of the K states of Mg(12) and Al(13). 


AU 


tains three maximum-minimum pairs, inserted 
in the spaces corresponding to the B—8, B-C, 
and y—D intervals of the Al curve. The appear. 
ance of these may be explained as due to the 
additional reflection planes (zones of discon. 
tinuity) made possible by the increase in com. 
plexity of the hexagonal crystal over the cubic 
crystal. 

These same two contours were recorded photo- 
graphically by Sandstrom’ and it is interesting 
to note that the first maxima observed there were 


w CORRECTED FOR 

- Epce w Not CorRRECTED RESOLVING POWER 
Mg 1.5 volts 0.5 volt 
Al 1.25 volts 0.35 volt 


17D. Coster, Physica 2, 606 (1935). 


ul 
1320 30 1300 i200 


Fic. 7. The absorption coefficient of Mg (12) in the 
neighborhood of the K-absorption limit. width of 
the K state is indicated by the line directly below th 
arrow. 


18 A. Sandstrom, Dissertation, Uppsala, Sweden, 1935. 
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for Al, a=27v and Mg, a=26v from the edge. 
These correspond to y and ¢, respectively, in the 


present Cases. 
A comparison between the K and L edges in 


crystals of the same type is also possible, for | 


while the selection rules involved in the two cases 
differ, the electron removed by absorption may 
be considered free. The crystal structure is the 
determining factor and the positions should 
agree quite close to the edge; the shape of the 
maxima and minima, depending on the transi- 
tion probabilities, is mostly determined by the 
lower state. 

TABLE IV. Positions of structures of Me 2) 


(reduced), Al(13) and Ag(47) measured in volts 
main edge of each curve, respectively. 
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J. Veldkamp, Physica 2, 32 (1935). 


ALUMINUM 
K EDGE 
795A 


1600 1550 VOLTS 


Fic. 8. The ie coefficient of Al (13) in the 
neighborhood of the K-absorption limit. The width of the 
K state is indicated by the dash below the arrow. 


Silver has a grating constant of 4.08A, or a 
distance to nearest neighbors of 2.89A, almost 
identical with Al, both having the same crystal 
structure. The Lm edge of Ag was recently 
recorded on a double crystal spectrometer” and 
in the final column of Table IV are given the 
distances of the maxima and minima from the 
edge for the Ag curve for all points corresponding 
to e on the Al curve. In spite of the nearness to 
the edge, the agreement is quite remarkable. 

The effect of nuclear (ionic) charge is also 
noticeable. Al has three conduction electrons, Mg 
two. The variation in potential through the 
lattice should be greater for Al, the effects of 
which show up in the curves in the sharpness of 
the fluctuations. 


2°. G, Parratt, Phys. Rev. 54, 99 (1938). 
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PHYSICAL REVIEW 


A Theory on Cathode Luminescence 


Uco Fano 
Washington Biophysical Institute, Bethesda, Maryland 
(Received May 31, 1940) 


A tentative explanation is given of experiments concerning the dependence of the brightness 
of a luminescent material on the voltage of electron bombardment. A model is assumed in which 
surface effects prevail over luminescence in using radiation energy. The possibility of a quanti- 
tative theoretical treatment is largely due to symmetry properties in the diffusion of scattered 


electrons in a solid body. 


1. INTRODUCTION 


B. BROWN! performed accurate meas- 
*urements of the amount of light emitted 

by a luminescent material known as “artificial 
willemite’’ (Zn2SiO,+about percent MnO; 
activator) when bombarded by cathode rays of 
200 to 1300 volts.2 Measurements were made 
varying both voltage and current density. 
Luminescent energy was emitted within a band 
at about 5500A. The important feature of the 
experimental result is expressed by the formula: 


W=5.10-*V?+0.17. (1) 


- W is the luminescent energy per impinging 


electron measured in electron volts, V is the 
voltage of the impinging electrons. This formula 
is valid for small currents and represents a very 
good approximation of the experimental results 
between 200 and 800 volts; above 800 volts W 
becomes slightly larger than indicated in (1). 
We shall call efficiency the fraction of the 
electronic energy converted into light, that is 
W/V. The efficiency is small: at 200 volts its 
value is 0.002; that is, one obtains an average 
of 2/10 of a light quantum per electron. If the 
slight correction arising from the term 0.17 is 
neglected, the efficiency increases proportionally 
with the energy. The term 0.17 will be neglected 
in the general lines of the explanation; its 
meaning will appear only in the final discussion. 
It is very striking to find a steady variation of 
efficiency following a simple formula, since the 
phenomenon itself is quite complex and since 

1T. B. Brown, J. Opt. Soc. Am. 27, 186 (1937). 

2An independent research by W. H. Nottingham (J. 
App. Phys. 8, 762 (1937)) confirmed Brown’s results and 
was extended to a wider interval of voltage, in which, 


however, the quantitative relations are not so impressive 
as in Brown's interval. 


the elementary processes involving slow elec. 
trons are not known to follow any simple 
quantitative law. 

Further information may be obtained from 
experimental results quoted by N. Riehl. This 
information concerns production of luminescence 
by a-particles in ZnS+activator Cu which takes 
place with very high efficiency (order of magni. 
tude 80 percent). The mechanism of the phe. 
nomenon is not supposed to be very different in 
the two cases studied by Brown and Riehl. 
Assuming that each atom of Cu can only emit 
one or two quanta under the influence of the 
same a-particle, Riehl calculates the number of 
Cu atoms which are needed to convert into 
light the energy of a single a-particle. From this 
number he calculates, then, that an a-particle 
is able to excite Cu atoms at a distance of 
4X10-5 cm from its path. 

It is striking to find a large difference in the 
order of magnitude of the efficiency of a-rays 
and cathode rays whose velocities are nearly the 
same. (Because of general arguments of the 
theory of the stopping power, the density of 
production of excited atoms cannot be very 
different if the velocities are nearly equal.) 


2. GENERAL LINES OF EXPLANATION 


We shall call A the effective range of diffusion 
of the energy transferred by charged particles to 


*N. Riehl, Ann. d. Physik 29, 636 (1937). 

‘The information given by a-particle experiments 
precludes any explanation of Brown’s small efficiency 
which would rely on the crowding of excitations along the 
path and on the fact that crowding decreases with in- 
creasing energy. Ultraviolet rays also produce luminescence 
with very high efficiency (about 50 pow however, 
no strong argument may be derived from this since the 
mechanism of luminescence may conceivably be quite 
different in the case of ultraviolet rays. 


544 


4 
| 
| 
| a. 
| dis 
for 
| far 
hay 
| Alf 
per 
cry 
strc 
pen 
pen 
cm. 
the 
muc 
the. 
a-pé 
its 1 
of el 
emis 
the 
po 
are 
inter 
effici 
it inc 
tion 
| 
| excite 
it me 
rule, 
motio 
transt 
symm 
the j 
irregu 
face 
irregul 
becaus 
We 
that e 
very | 
proces: 
surface 
= 


a8 


=. 


Bea ers? 


s 5 


a luminescent crystal. Thus A is the maximum 
distance at which charged particles are able to 
excite luminescent atoms. It is reasonable to 
assume that A is of the same order of magnitude 
for a-particles and slow electrons and not very 
far from the value 4X10-° cm calculated by 
Riehl. 

The energies of a-particles and slow electrons 
having the same velocities are quite different. 
Alpha-particles are not strongly scattered and 
penetrate very deeply (~3X10-* cm) into a 
crystal before being stopped ; slow electrons are 
strongly scattered and readily stopped before 
penetrating deeply: One may guess that the 
penetration of 200-volt electrons is about 10-7 
cm. Hence the penetration of cathode rays into 
the crystal is much smaller than A, whereas it is 
much larger in the case of a-particles. Therefore 
the diffusion of energy takes place in the case of 
a-particles far away from the surface, whereas 
its very origin lies near the surface in the case 
of electrons. 

Parasitic processes, as production of heat or 
emission of secondary electrons, may prevent 
the diffusing energy from reaching luminescent 
atoms. We shall assume that parasitic processes 
are much stronger on the surface than in the 
interior of the crystal. This explains that the 
efficiency of slow electrons is very small and that 
it increases with increasing energy and penetra- 
tion of the electrons. Our assumption cannot 
be proved, since sufficient knowledge of the 
excited states of crystals is not available, but 
it may be easily made plausible. As a general 
rule, the coupling between electronic and nuclear 
motion in solid matter is very effective in 
transforming excitation energy into heat; special 
symmetry and periodicity relations may reduce 
the intensity of this process strongly, but 
irregularities of symmetry or periodicity will 
conversely destroy any metastability. The sur- 
face of a crystal constitutes a very important 
irregularity, not only as a discontinuity, but also 
because the surface itself may not be smooth. 

We shall first make the particular assumption 
that energy reaching the surface by diffusion is 
very likely to be lost for the luminescence 
process. The actual probability of loss on the 
surface will be deduced as an experimental result. 
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‘In §3 we shall consider the diffusion of energy 
Excitation energy may be transferred from an 
atom to the neighboring ones in the lattice; 
wandering excitations may be treated as diffusing 
particles, called excitons. Excitons are scattered 
(elastically or inelastically) by any irregularity 
in the lattice: thermal motion, foreign atoms and 
so on. The energy of the excitons is absorbed in 
one or several steps. We shall assume that the 
efficiency of an exciton is proportional to its 
probability of diffusing to a distance A without 
having previously reached the surface. Elemen- 
tary considerations of diffusion theory show then, 
that the efficiency is proportional to the distance 
of the origin of the excitons from the surface and 
that the proportionality coefficient is about 1/A. 
Hence one obtains agreement with the experi- 
mental efficiency 10-* for 200 volts, taking as 
reasonable figures A~10~* cm and the distance 
of the origin of the excitons from the surface 
~ 10-7 cm. 

One will then investigate (§4) how the 
“average depth’ of production of excitons 
depends upon the initial energy of the electrons. 
It is helpful to apply symmetry considerations 
to the process of multiple scattering and diffusion 
of the electrons. The ‘average depth’”’ may be 
expressed in terms of simpler physical quantities : 
free path of scattering and stopping power of 
electrons, and “effective albedo’’* of the material 
for impinging electrons. 

Free path of scattering and stopping power 
will require some particular discussion, since the 
approximations generally used are rather inac- 
curate within our range of energies. Knowledge 
of the “effective albedo”’ would require a detailed 
solution of the diffusion problem, but one may 
show that its influence on the results is very 
small. One will finally show that the average 
depth of production of excitons is nearly pro- 
portional to the electronic energy (§5). The final 
formula may be fitted quite well to Brown’s 
experimental results (§6). 


* The properties of excitons may vary within a wide 

range, since they are very sensitive to quantitative relations 
of the particular crystal lattice (reference may be made 
toa od by J. Frank and E. Teller, J. Chem. Phys. 6, 
861 (1938)). 
_ *Following the analogy with neutron diffusion, we 
indicate as “effective albedo” a physical quantity con- 
cerning that fraction of the impinging electrons, which 
emerges again from the surface after diffusion. 
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3. ExcitonN DIFFUSION AND SURFACE 
ABSORPTION 


To determine the efficiency 7 of an exciton of 
energy ¢ produced at a small distance x from the 
surface, we assume first that the exciton is lost 
if it reaches the surface. We conclude from 
a-particle experiments that the efficiency is 
practically unity if the surface losses are avoided. 
The efficiency 7 is then equal to the probability 
p(x), that the exciton does not reach the surface 
and we want to calculate this function under 
the assumption: p(x)<1. A second exciton 
produced at a depth 2x has the probability p(x) 
of never coming within a distance x from the 
surface; if, however, it does come within such a 
distance—the probability for this is (1—p(x))— 
it still has the same probability p(x) of escaping 
as the exciton considered first. Hence we find a 
total probability: p(2x)=p(x)+(1—p(x))p(*) 
~2p(x), neglecting p(x)*, and we deduce, that 
n= p(x) is proportional to the depth: 


n~kx if kx<1. (2) 
Moreover, since »~1 for x~A, we have: 
k~1/A, n~x/A, if x<A. (3)? 


Calling w the average amount of energy radiated 
as luminescent light, one has: 


w=ne~xe/A. (4) 


We shall now take into account two factors 
neglected heretofore. First the probability of 
an exciton’s reaching the surface and being 
reflected by it without absorption is the product 
of the probability of its reaching the surface, 
and the reflection coefficient r (which is still 
supposed to be small). This product is nearly 
equal to r, since the first factor is (1—x/A)~1 
within the limits of our approximation. Second, 
we have so far considered the free path of 
diffusion / as infinitely small as compared with x 
(according to the usual assumption of diffusion 
theory). If 10, but still <x, one has to add to 
(3) a constant of the order of //A. Then, putting 


Cas: 
C=1/A+r (S’) 


7 The formula without the condition x/AX1 is of the 


t : 


one has: 
n~x/A+C; w~xe/A+Ce. (5) 


The generalization to the actual Productiog 
of many excitons by the electrons is carried out 
by introducing the amount of energy €(x)dy 
lost by the electrons between the depths x ang 
x+dx. Hence (5) goes over to: 


f 
W=— ——+€ J 


which is the theoretical expression to be com. 
pared with (1). The physical meaning g 
So°xe(x)dx is evidently : the excitation energy time 
the average depth at which the excitation takes 
place. 


4. SYMMETRY PROPERTIES OF ELECTRON 
DIFFUSION 


Electrons entering a solid body undergo 
multiple scattering, which goes over into diffusion 
and may eventually lead to the re-emission of 
the electrons. This phenomenon is rather compli- 
cated by the parallel action of diffusion and 
energy losses; slowing down changes the prob. 
ability of scattering, as the diffusion goes on, 
The following qualitative considerations of 
simple diffusion phenomena will be helpful for 
the calculation of J/o°xe(x)dx. 

We shall consider the diffusion of particles 
with unit mean lifetime produced at a rate of V 
per second (their total number is then JN). The 
particles diffuse in a three-dimensional space 
without boundaries and are produced along the 
plane x=a. Call m(x)dx the number of such 
particles between the planes x and x+dx; we 
want to discuss the average value of the x o 
ordinate of the particles, thatis, (1/N) /_2xn(x)de. 
The symmetry of the system requires n(a+i) 
=n(a—6), and hence: 


f xn(x)dx=aN. (7) 
We shall now distinguish as particles d 


“class 1” those particles which have already 
crossed the plane x=} (b<a) at least once; al 
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(6) 


the remaining particles shall be called of ‘‘class 
2,” so that: The particles of 
class 1 (iin number) are distributed as particles 
produced along the plane x=); hence: 


f (7’) 


Since no particle of class 2 is found to the left of 
x=b, we have: 


aN= f f 


+ f f xn2(x)dx. (8) 
6 


We consider now particles which are produced 
at x=a and which are lost if they reach the 
surface x=0. Their distribution is represented 
by n(x)=m2(x), if one takes b=0; then (8) 
becomes: 


f =aN. (9) 
0 


The integral on the left side of (9) is independent 
of the number of surviving particles Jo°n(x)dx: 
the more particles flow out through the surface 
x=0, the farther the survivals diffuse in the 
interior towards x= 

If one takes into account that the value of 
the free path \ of the particles is not negligible 
one cannot consider the particles of class 1 as 
if they were produced at x=) because they are 
all directed towards x= — when they cross 
this plane. The coherence of their directions is 
lost at the end of their free path; we have then 
to substitute in (7’): 


b—b—X{|cos 8! (10) 


'The mathematical treatment of the same problem is 
this: Call g(x) the distribution of production of the 
particles, r their mean lifetime, D the diffusion coefficient 
of the icles. n(o)=0 with its derivatives is the 
general boundary condition; »(0)=0 is the effect of the 
absorption at x=0. Then: 


— n(x) +9(x)=0 
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X 


Fic. 1. m(x): Plot of the approximate distribution ‘of 
particles originating along the plane x=a. mi(x): Plot of 


the approximate distribution o icles having already 
the plane m2(x): Plot of the approximate 
ene of particles which have not yet crossed the 
plane 


where @ is the angle between the directions of 
the particles at x=} and the x axis, and where 
{|cos 6|)w is the average of the absolute value 
of cos 6. The same method may be applied to 
the case of a nonisotropic source of particles, 
for instance directed towards x = © , substituting : 


a—a+X!cos (10’) 


We consider now the diffusion within a body 
extending from x=0 to x= ©; particles enter 
the boundary surface x=0 and are lost, if they 
flow out again through it. In this case we have 
to introduce the substitutions (10) and (10’) 
into (8) and then to put a=b=0. Calling J the 
current of particles, we have then: 


f xn(x)dx = J inAin( | cos Gin | Yaw 
+ Joutdout( | COS Bout | (1 1) 


The average depth of the particles is thus de- 
termined by the flow of particles through the 
boundary and by their free path, independently 
of the actual details of the diffusion problem. 
The problem of finding /o°xe(x)dx for electrons 
scattered in a solid body is still more complex. 
One may consider electrons having been slowed 
down as classified according to their energy; 
slowing down shifts electrons from one class to 
the next one. One finds that symmetry con- 
siderations still apply to each class so that a 
final result analogous to (11) may be found. 
Anisotropy of scattering has also to be taken 
into account; this must be done here—as in 
many other problems—by dividing the actual 
mean free path by 1—(cos #),, (8 =deflection). 
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We shall now go over to quantitative calcula- 
tions. We introduce a system of polar coordinates 
6, ¢, having x as polar axis, and call F(x, s, @) 
the current of electrons of residual “‘true’’ range 
s * at a depth x with direction 6, per unit residual 
range, per unit solid angle. We call # the angle 
between the directions (6, ¢) and (6’, ¢’), o(#)dQ’ 
the differential cross section for elastic scattering 
from (6, ¢) into (6’, ’) within the solid angle dQ’ 
and N the number of scattering atoms per cc. 
F(x, s, 0) is a solution of the equation: 


OF 
cos @—=—— NF(x, s, af 
ax OSs 
F(x, s, 0’) o(8)dQ’. (12) 
4x 
Putting: 


f cos"6@ F(x, s, @)dQ= F,,(x, s) 


4n 


and taking into account the properties of the 
spherical harmonics, one has: 


OF, OF; 
—=——— | o(8)(1—cos #)dQ’ 


Ox Os 
OF, 


1 
=—-Fr, (13)" 
Os 


where: 
1/A=NSiz0(8)(1 —cos 8)d0’, 


and hence: 


* OF a(x, s’ ds!" 
Fi(x, s) = f - f —| (14) 


*“True” range is the range measured along the actual 
path. We consider s as a suitable measurement of the 


energy of the electrons. 
10 for an analogous equation: W. Bothe, Zeits. f. 


Physik 54, 161 (1929). 

11 Expansion in spherical harmonics may be considered 
as a general method of attack to the integral-differential 
equation (12) and transforms it into a chain of differential 


equations ; 


F(x, s, 0)Pi(cos 


F,(x,s) represents the positive current of 
particles with residual range s across the plane 
So’ E(s)Fi(x, s)ds is the total stream of kinetic 
energy across it towards x= «© where E(s) jg 
the energy corresponding to the residual range 5. 
Hence the negative gradient of fo" E(s) Fi(x, s)ds 
represents the density of energy loss of the 
electrons; 


J 


Integrating by parts, using (13), changing the 
order of integrations and using the condition 
F=0 for x= ~, we obtain the final formula: 


dx= | dsF,(0,s)&(s), 
J xe(x)dx J (0, s)&(s) 


&(s)= J ‘ds'E(s") exp| - J 


This formula closely corresponds to (11), 
since F,(0, s) may be represented as: 


Jin{ | cos Bin | w+ Jout{ | Cos Bout | Dave 


(15) 


&(s) replaces \; its dimensions are: energy loss 
of electrons times depth of their penetration, 
The dependence of &(s) on the elements de- 
termining the penetration, scattering cross 
section and stopping power dE/ds, appears more 
clearly in the asymptotic formulae: 


&(s)~E(s)Xs), if A(s)<s, (16) 
&(s)~E(Sav)s, if d(s)>s. (17) 
A better approximation than (16) is: 


16 


For the second integral in’ (6) one has, obvi- 
ously: 


f f "asE(s) Fx(0, (18) 
0 0 


5. DiscussION 


We have now to collect the evidence available 
at present concerning the dependence of the 
quantity (15) on the initial energy of the 
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electrons. Equation (15) is given as a function 
of the free path A(EZ)=X(s), of the stopping 
power dE/ds (appearing as E(s)) and of F;(0, s). 
The quantity F2(0,s) is an integral over 6 and 
consists of the contributions of the impinging 
electrons (@< 7/2), which is known, and of the 
outgoing ones (02/2). The unknown features 
of the diffusion” affect only the integration for 
62 1/2. 

The contribution to (15) arising from the 
incoming electrons is proportional to the square 
of their voltage V? if any one of the following 
conditions is fulfilled : 

(a) if \(s)s and \(Z)~const. E£, on account 


of (16); 
(b) if A(s)>s and dE/ds~const., on account 


of (17); 

(c) if both \(Z)~const. E and dE/ds~const., 
whatever the average of \ and s may be, since 
then, putting \=aE, s=BE: 


B 
= 
2a+B 2a+B 


The contribution of incoming electrons is 
sufficient to determine the dependence of (15) 
on V completely ; 

(A) if the contribution of the outgoing elec- 
trons is negligible; 

(B) if the contribution of the outgoing elec- 
trons is proportional to the contribution of the 
incoming ones and the proportionality factor is 
independent of the energy. 

Fulfillment of one among the conditions (a), 
(b), (c) and of (A) or (B) is sufficient to show 
that (15) is proportional to V*. Proportionality 
of the second term of (6), that is of (18), to V 
is also a consequence of fulfillment of (A) or (B). 

We shall next discuss how (A) and (B) are 
connected with (a), (b) and (c). Condition (b) 
leads automatically to (A), since then practically 
no electron is scattered back. Condition (a) 
involves an approximation opposite to that of 
(b); however, (a) is consistent with (A), if the 
relative energy loss along a free path is not 
negligible. In this case every free path is shorter 
than the former one; thus an electron is pre- 


(19) 


_ ™The determination of the outgoing part of F(0, s, 0) 
is analogous, but more complex, to that of the “albedo” 
of neutrons impinging on a paraffin block. 


vented in most cases from reaching the surface 
of the crystal soon after its first penetration 
with a large residual energy. (A) is then fulfilled, 
since the contribution of outgoing electrons to 
(15) depends strongly upon their energy (possibly 
upon its square, according to (a)).'* Condition 
(B) is fulfilled if: 

(B’) o(#), as a function of energy, is inversely 
proportional to s™ and the relation between 
energy and range has the form E=as* where a 
and d are constants. The proof of this is obtained 
through a similarity transformation of (11) 
showing that, if (B’) is fulfilled, a change of 
initial energy affects the solution F of the 
diffusion equation only as a change of units of x 
and s. Condition (c) together with the additional 
assumption that the dependence of the angular 
distribution of scattering on the energy may be 
neglected, is sufficient to fulfill (B’), so that (B) 
follows then as a consequence. The additional 
condition is required, since proportionality of \ 
to s is less stringent than the condition o(#) 
=const./s. 

Going over to discussion of the quantitative 
relations involved in (a), (b), (c), we have first 
to evaluate by a wave-mechanical calculation 
the effect of shielding on scattering of slow 
electrons. The scattering cross section in a 
central field is: 


h?2 


= (21-+1)Pi(cos 1) |2, 
l=0 (2 


0) 


(m=electron mass, E=electron energy, 4 
= phase shifts of the wave functions determined 
by field and energy). 


1 
—cos 3)dQ 


Ni? 


1 sin? (6:;— (21) 


% This picture is also supported by the experimental 
evidence, that all but a small fraction of the electronic 
energy is absorbed in a solid body. See the results discussed 
by W. Bothe (Handbuch der Physik, Vol. XXII/2, 1 
(1933)), particularly those of B. F. J. Schonland (ébid., p. 
44) ; very soft re-emitted electrons are not taken into account 
in these experiments. Secondary electrons emitted by the 
surface are considered in our theory as “‘lost excitons.” 

“This condition is fulfilled for instance (neglecting 
logarithmic factors) in the case of high energy electrons, 
for which Rutherford scattering is inversely proportional 
to E’, while the range is directly proportional to it. 


of 
1e x, 
etic 
) is 
re 
s)ds 

the 

the 
tion 
loss 
‘ion. 

de- 
TOSS 
(16) 

able 

the 

the 


550 UGO FANO 


30} 
Zn, SiO, 
20) 
IS 
=~ 
ws 
10 
Os 
5 


0 "700 200 300 400 500 600 700 800 
Volts 
Fic. 2. The effective average cross section for elastic 
scattering of electrons: ¢= is given 


by the square of the wave-length of electrons times 
sin? by x. (From calculations based 
el. 


on Thomas-Fermi 


Constant 6,’s (and therefore constant \/E) are 
the consequence of a potential field 1/r* for all 
distances from the nucleus. The 6, for a Thomas- 
Fermi field have been calculated by Henneberg.” 
For high energy electrons all 6;’s are very small 
and decrease with increasing energy; therefore 
the terms sin? (6;—4::) also decrease with in- 
creasing energy. One finds \ to be proportional 
to E*. At E=0 all 6,’s are multiples of x and \ 
tends to a finite value. In the intermediate 
region (which is the interesting one) 


/ sin? 


—that is E/\—goes through a flat maximum. 
Inasmuch as this function may be considered as 
a constant in the neighborhood of its maximum, 
one finds the relation \~const. E."* Because of 
the characteristic features of a Thomas-Fermi 
field, only the waves corresponding to the low 
values of / !” are strongly distorted, giving rise to 


15 W. Henneberg, Zeits. f. Physik 83, 553 (1933). 

16 The Thomas-Fermi potential field is: Ze*/r-o(r/u) 
=Ze*/r?-re(r/u); the function xg(x) has a very flat 
maximum and varies slightly if one changes the distance 
from the nucleus within a factor 10. Further evidence on 
the inverse linear dependence of scattering from the energy 
is obtained by plotting Lenard’s experimental data of 
a ent absorption coefficients, as given by W. Bothe 
(Handbuch der Physik, Vol. XXI1I/2 (1933), p. 41). 

17 These values correspond essentially to the possible 
values of / for atomic electrons (see for instance: E. Fermi, 
Zeits. f. Physik 48, 73 (1928)). 


large values of their 6;; hence the first 
sin? (6;— 5:1) are of the order of 1 and the sum 
converges rapidly after a certain value of | 
which is characteristic for the atomic number 
but does not depend upon the energy of the 
scattered electron. 

Figure 2 represents the values of : 


lsin? (6,;— 61-1) Jan sin? 5:1) Js; 
sin? (8:— 814) 


and of its single terms as functions of the 
energy, deduced from Henneberg’s results by 
graphical methods. Some approximations in 
these or analogous calculations could hardly be 
avoided. The result is fairly constant, since jt 
varies only by 28 percent within the range of 
E=200 to 800 volts; as a consequence, } jn. 
creases in this interval by a factor ~5.6 as 
compared to 4 in the case of a 1/r* field and 16 
in the case of a 1/r field. Variations of E/) are 
chiefly due to the large amount of oxygen 
contained in the material, since the typical | ir 
region corresponds for oxygen to energies lower 
than 200 volts. Henneberg’s calculation neglects 
the effects of electronic exchange and of polar. 
ization of the internal electronic shells. Polariza- 
tion is not expected to influence our result; 
electronic exchange might be effective in im- 
proving the constancy of \/E appreciably, 
Both effects are difficult to take into account 
quantitatively. 

Through discussion of the stopping power as 
a function of energy, we wish to justify two 
fundamental qualitative assumptions and two 
subsidiary ones: (1) the range s is definitely 
larger than A (this is also known from experi- 
ments) and hence (a) is nearly fulfilled; (2) the 
relative energy loss along a free path [~(h/B) 
x (dE/ds) ] is of the order of 3, hence (A) may 
be fulfilled together with (a); (3) the stopping 
power is fairly constant over our range of 
energies, therefore (c) as well as (a) is nearly 
fulfilled; (4) the ratio of the free path to the 
“true” range and the relative energy loss along 
a free path are fairly constant over a wide 
range of energies, thus also (B) is not very far 
from being fulfilled and one feels—though 
indirectly—encouraged to use formula (19). 
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One can hardly evaluate the stopping power 
in our range of energies, since every theory 
assumes that the particle to be stopped is much 
faster than atomic electrons, while this condition 
is not fulfilled in our case. Considering, however, 
that we need only qualitative results, we shall 
tentatively use Bloch’s formula’* in the form: 


dE h? Eo E 
—=0.385N— — log —; 
ds mE Eo 
(22) 


6.04 
Z-13.54~11.2Z volts. 
\e 


This formula relies on Thomas-Fermi model and 
might be not too inaccurate for slow electrons, 
since it has been applied successfully to a- 
particles having a velocity equal to that of 
900-volt electrons. Comparison with experiments 
is difficult in this region ; experiments on ionizing 
action in gases are available’® and show a more 
irregular dependence on energy than (22) does. 
The zero point of (22) at E=£p is obviously 
incorrect but a rapid improvement in accuracy 
may be expected for increasing energy. For 
E/E>>1 the variation of the log factor is 


3.0 


0 300 600 900 1200 1500 
Volts 
Fic. 3. Stopping power of willemite, calculated from 
Bloch’s formula, assuming density 4. The contributions 
of the single elements are also indicated. 


4 F. Bloch, Zeits. f. Physik 81, 363 (1933). It might be 
questionable, which way one should take into account the 
effects arising from the largest energy transfers to atomic 
electrons. Formula (22) contains elements from Bloch’s 
formula and from a formula by H. A. Bethe (Handbuch der 
Physik, Vol. XXIV/1 (1933), Pp. 521 (56.10)). This is because 
the characteristic of Bloch’s formula concerns large values 
of the impact parameter, while the factor (e/2)# concerns 

impacts. Anyway, a change of this factor would 
probably not affect the results too much. 
287 (18 Ay instance: R. B. Brode, Rev. Mod. Phys. 5, 
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negligible and one has the well-known result 
s~const. E?. A maximum occurs at E= Eye~30Z 
volts, which is very flat. Hence dE/ds is nearly 
constant over a wide region (point (3) above). 
This last fact is also due to the different position 
of the maxima for different elements. Figure 3 
represents the stopping power of Zn,SiO, and 
of its single components. The result is certainly 
very inaccurate below 600 volts since Ep» is 
~ 336 volts for Zn; the stopping power should 


0.05, 


200 


1000 (200 
Volts 


Fic. 4. Correspondence between the “true” residual 
range of electrons and their energy and between “effective” 
free path and energy. The ratio (A/E)/(ds/dE) roughly 
ae the fraction of electronic energy along a 
ree path. 


probably be larger in this region than is indicated 
in the figure. 

The quantities \/E and ds/dE, which are 
directly comparable, and, if constant, should be 
equal to the parameters a and £ of (19), are 
plotted in absolute units in Fig. 4. Points (1), 
(2), and (4) are verified by inspection of this 
figure. The asymptotic behavior of both \/E and 
ds/dE for large energies (1/r field) is proportional 
to the energy and their ratio is then 4Z,,/(Z*) 
=1/5.3; this result gives further evidence for (4). 

Figure 4 shows that \/E may be fairly well 
represented by a straight line between 200 and 
800 volts; then: 


(23) 


where E is expressed in volts. 
The results (1), (2), (3), (4), indicated above 
allow us to use any one of the formulae: 


& E(s)X (: 16’) 
&(s)~ E(s)Xs). (19) 

2a+B 
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Remembering our former results it is easily 
seen that the factors multiplying EX in both 
formulae are equivalent and that their value is 
approximately 0.75. Putting then (16’) or (19) 
into (15), we have: 


f 
X(1+7.3X10-*V) V2 voltscm. (24) 


‘ includes: the correction 0.75, as mentioned 
above, a factor <1, if the incidence is not 
normal, and the possible contribution (factor 
>1) of outgoing electrons to (15); 71 is expected 
to be of the order of 3. We call y2 the fraction 
of incoming energy absorbed within the crystal, 
i.e., not taken away by re-emitted electrons; 
‘2 is constant if (A) or (B) is fulfilled and has 
to be very near to 1. Then (6) becomes: 


V? 
+y2CV. (25) 


6. CONCLUSION 


‘The results of the quantitative discussion 
seem to justify the following picture of the 
luminescent phenomenon. The average penetra- 
tion of slow electrons in a crystal is essentially 
given by their effective free path of scattering \ 
(defined in formula (13)). Because of character- 
istics of the atomic field as represented by the 
Thomas-Fermi model, \ is nearly proportional 
to the energy of slow electrons. A better approxi- 
mation is given in Fig. 4 and formula (23). 

Slow electrons undergo multiple scattering 
and diffusion within the crystal, but it is not 
necessary to calculate this process in detail. 
One may calculate directly /¢*xe(x)dx, which is 
the product of the energy lost by the electrons 
within the crystal and the average depth of 
energy loss. A careful discussion is necessary to 
show, that /o*xe(x)dx may be considered as 
proportional to the voltage of the electrons and 
to their penetration \. The efficiency of pro- 
duction of luminescence is proportional to the 
average depth of production of excitons (i.e., of 
energy loss by electrons) as stated in §2. Thus 
one finds finally, that the efficiency may be 
considered as a linear function of the penetration 
\ of the electrons. Formula (25) is then obtained, 
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Fic. 5. Luminescent energy emitted by willemite wher 
bombarded by cathode rays, as a function of electrop 
energy, expressed both as brightness and as energy output, 
The theoretical curve (26) agrees with the experimental 

wo parameters have been adjusted to fit the points at 
200 volts and 800 volts. . 


which gives the amount of luminescent energy, 
that is, the efficiency times the incident energy, 

v1 and y2 may be considered constant and 
slightly smaller than 1. A is the range of diffusion 
of excitation energy within the crystal. C js a 
physical constant defined by formula (5’) and 
still to be determined. The third term of (25) 
corresponds to the constant term in the linear 
representation of the efficiency as a function of 
the penetration of the electrons (see formula (5)), 

Comparison of (25) with (1), neglect of the 
second and third term of (25) and the second 
term of (1), gives: 


Aw~1.1y:10-*~0.7 X 10-4 cm. 


This comparison corresponds to a rough theory 
aiming to show that the efficiency is proportional 
to the energy of impinging electrons. The 
agreement is improved, if one tries to compare 
the complete expressions (25) and (1) by 
suitable determination of the unknown constants 
A and C: 


aV?+7.3X10-4a V?+5V—5 XK 10-* V2+0.17. 
Numerical evaluation shows that these two 


analytical expressions may be made neatly 
coincident over a wide range of values of V by 
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suitable choice of a and b. If we take, for instance: 


X 10-8, 


(26) 
b=Cy2=0.133, 


in order to fit the points at 200 and 800 volts 
exactly, that is: 


A~1.6X10-*cm, C~0.14, (26’) 
then the discrepancy between (25) and (1) never 
becomes larger than 0.03 volt in absolute 
value, or 5 percent in relative value, between 
200 and 800 volts (Fig. 5). This agreement is as 
good as one could expect. Above 800 volts (25) 
correctly represents the faster increase of W 
(see §1). Equation (25) is certainly not very 
reliable below 200 volts, but it is important 
that the theory indicates W=0 for V=0O since 
the extrapolation of (1) to V=0 is obviously 
incorrect. 

The interpretation of the experimental results 
now becomes slightly different from the simple 
one outlined in §2. The theory does not consider 
that the experimental points of the plot W(V?) 
lie on a straight line as represented by (1). 
Instead it yields an S-shaped curve which fits 
the experimental points as well as a straight line 
between 200 and 800 volts and still better than 
a straight line above this interval.?° The result 


*No important deviation of \ from (23) is likely to 
take place above 800 volts, since the slope of \/E already 
corresponds to its asymptotic value given by the theory 
of Rutherford scattering. 
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is essentially due to the approximate constancy 
of \/E, but the apparent linear dependence of 
W on V? is due to the opposite action of two 
different factors. The slight increasing trend of 
\/E tends to bend the curve W(V*) upwards, 
whereas the constant term C in the efficiency 
tends to bend it downwards. The accuracy of 
the agreement with the experimental results 
appears not to be due to accidental relations 
among the numerical constants, since it is only 
slightly affected by rather large changes of the 
numerical values (26’). 
From the formulae (5’), (26’): 


C=1/A+r~0.14 


we may obtain further information. The free 
path of diffusion of excitons 7 should not be 
larger than A, according to the assumptions of 
our theory ; now, since \/A0.14, one concludes, 
that: 

1/AKC~0.14. 


Hence we obtain an evaluation of the probability 
that an exciton impinging on the surface is not 
absorbed by it: 


r~C~0.14. 


I am very much indebted to Dr. E. Teller, 
who suggested this work to me and supported 
me with encouragement and advice in many 
discussions. Dr. T. B. Brown kindly furnished 
me with additional information and data on 
his experiments. 
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ar 
Cross sections for capture of thermal energy neutrons were measured for several elements. m 
The method consisted in comparing, by a BF; ionization chamber, the density of thermal co 
neutrons in water and in a solution of the element investigated. In order to avoid corrections for la 
geometry, the effect of the substance was compared with that of boron solutions of different 
concentrations, and consequently the measurements directly give the ratio of the cross sections 
of the elements investigated to the cross section of boron. Data are given for N, F, Na, P, Cl, K, 
V, Mn, Co, As, Se, Br, Cd, I, Ba. 
de 
. an 
OST of the experiments on scattering and captured, either by hydrogen or by boron un 
absorption of slow neutrons give values of nuclei. The same method would be unpractical, we 
the total (scattering+capture) cross section. when applied to most other substances, since jt eit 
When the cross section is of the order of 10-* cm? would require exceedingly large volumes of cop. ter 
or more, it is well known to be almost wholly due _ centrated solutions. This can, however, be fro 
to capture. Our present knowledge about the avoided if the cross section of the element to be net 
character of the cross sections lying between investigated is measured in terms of the capture caf 
10-* and 10-* cm? is much less definite. Measure- ‘cross section of boron. An amount of solution of 10° 
ments of the scattering cross section have been any volume and shape may be taken, provided ‘ 
performed,! but in most cases they are not ac- one measures the concentration cg of boron which con 
curate enough to give a significant value of the produces the same density decrease of the anc 
capture cross section by subtraction from the thermal neutrons as the element investigated eler 
total cross section measured by neutron beam’ under the concentration c. Then the relation con 
experiments. holds: den 
whe 
EXPERIMENTAL PROCEDURE 
To perform a direct measurement of the where o and op, respectively, represent the all 
capture cross section, the simplest method which ©@Pture cross sections of the element in question leve 
suggests itself is the one already employed by 2"4 of boron. ; - 
Frisch, von Halban and Koch.? It consists in _ he neutron source consisted of 250 mg Ra = 
measuring, by a BF; ionization chamber, the Be BF; 
decrease in density of the thermal neutrons 
boron disintegrations were recorded by means of the 
inside a certain volume of water when a given : : N 
a linear amplifier and a thyratron scale of four. 


percentage of an shecrbing clement io added. The solution was contained in a cylindrical glass fact 
ont beaker, 14 cm in diameter and filled to a height pare 
q oron ane of about 20 cm, the ionization chamber being conc 
of the hed to be immersed in the center of the solution. The first 
taken large enough so that practically all  outron source, shielded by 2 cm of lead, was Penc 
neutrons emitted by the source were slowed placed outside of the solution, on the axis of the tion: 


q d down to thermal velocities, and eventually cylinder. The number of counts, under such con- (abo 
} 1M. Goldhaber and G. H. Bri 4. See ditions, was already over one thousand per of th 
= (1938) Mitchell and EJ, Murphy, minute. The effect of neutrons of energies higher 
| ; ood M. D. Whitaker, Phys. Rev. 49, 870 (1936); 50,133 than thermal in the present experiment & bein 

I negligible and was not corrected for. The ioniz- 

; 20. R. Frisch, H. v. Halban and J. Koch, Proc. Danske : ol 

i Vidensk. Selskab 15, ‘fo (1938). J: tion chamber and the glass cylinder were rigidly 
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connected together, and the investigated solu- 
tions could be quickly drained through a stopcock 
at the bottom of the container and replaced with 
water. Many alternate measurements with water 
and with the solution were taken for each ele- 
ment. Each reading comprised about 20,000 
counts and the discrepancies observed were no 
larger than the statistical fluctuations. 


RESULTS AND DISCUSSION 


Since large cross sections have already been 
determined by beam absorption measurements, 
and cross sections smaller than 10-* cm? are 
unpractical to measure by the present method, 
we confined our attention to elements which, 
either from the values of the total and the scat- 
tering cross sections given in the literature,’ or 
from the intensity of activation under slow 
neutron bombardment, were expected to possess 
capture cross sections lying between 10-* and 
cm?. 

The compounds employed had to answer the 
conditions of being sufficiently soluble in water 
and of not containing other strongly absorbing 
elements besides the one to be investigated. The 
concentration was chosen such as to reduce the 
density of the neutrons by at least ten percent, 
when allowed by the solubility of the compound. 

Most of the elements investigated (probably 
all except Cd) do not possess neutron resonance 
levels for thermal energies, and therefore the 
cross section in this region will follow the 1/v 
law. The ratio of the cross sections to the cross 
section of boron therefore possesses a definite 
meaning, being independent of the velocity of 
the neutrons. 

No corrections were required for geometrical 
factors, since the effect of the solution was com- 
pared with the effect of a boron solution of known 
concentration under the same conditions. It was 
first determined how the neutron density de- 
pended on the concentration of boric acid solu- 
tions. Within the limits of experimental error 
(about one percent) it was found that the ratio 
of the intensity J) with water to the intensity J 
with a boron solution of concentration cp 
was a linear function of cg, the empirical relation 
being 

X10*cp/cu, 


where with cq we have indicated the hydrogen 
concentration. In a solution of infinite volume, 
the coefficient k should be equal to the ratio of 
the capture cross sections of boron and hydrogen. 
The smaller measured value evidently depends 
on the fact that about half of the thermal neu- 
trons diffuse out of the solution before being 
captured, either by hydrogen or by boron nuclei. 

In most cases, the amount of substance which 
had to be dissolved in water in order to absorb 
a measurable fraction of the neutrons was large 
enough to alter to an appreciable extent the 
hydrogen concentration. Since a change in 
hydrogen concentration will, in. general, affect 
the density of the thermal neutrons, measure- 
ments were undertaken to determine the mag- 
nitude of the effect. It would not be correct to 
compare water with hydrocarbons or other 
hydrogenous substances, since a difference in 
the chemical bond of hydrogen might produce 
effects of the same order of magnitude. The pro- 
cedure employed was to dissolve in water a 
substance which would decrease the hydrogen 
concentration without appreciably absorbing 
the neutrons. The experiment was performed by 
using either a 30 percent solution of HO, 
(hydrogen density, 5.5 percent lower than in 
water) or a saturated solution of oxalic acid 
(hydrogen concentration, 3.7 percent lower than 
in water). The density of neutrons was found to 
be lower than in water by 4.3 percent in the 
first case and by 3.1 percent in the second case. 
As an average, we assumed that a decrease of 
one percent in hydrogen concentration produced 
a decrease of 0.8 percent in neutron density. The 
correction has been applied to the values reported 
in Table I. 

Part of the decrease in the number of thermal 
neutrons when an absorbing substance is dis- 
solved in water will be due to capture of the 
neutrons through resonance levels before they 
have been slowed down to thermal energies, and 
a correction may be required therefor. In order 
to estimate the importance of the effect, the 
following experiment was performed. A rhodium 
plate 5 by 5 cm, shielded by Cd, was substituted 
for the BF; chamber, and irradiated. Under such 
conditions, it is well known that rhodium is prac- 
tically activated only by neutrons of about one 
volt energy. The activity of rhodium was 
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measured (by means of a pressure ionization 
chamber connected toan Edelmann electrometer) 
when the container was filled either with water 
or with the solutions that had been investigated. 
For most elements, the characteristic resonance 
energies lie higher than the resonance energy of 
rhodium, this applying in particular to As, Br, 
and I. Then, the observed decrease in the 
number of characteristic rhodium neutrons will 
determine the effect of resonance absorption by 
the dissolved element in decreasing the number 
of neutrons which are slowed down to thermal 
velocities. Measurements performed with Cl, Co, 
As, Br, I, indicated that the decrease in activity 
of the rhodium was certainly smaller than 2 
percent, and probably within one percent. Some 
of the above-mentioned elements (As, Br, and I) 
would be likely to produce a more considerable 
effect than others, since it is known* that their 
activation is due, for an unusually high per- 
centage, to resonance neutrons. The effect being 
hardly measurable, even in these cases, we 
decided not to correct our results for resonance 
capture. 

A certain disagreement exists in the litera- 
ture?* about the cross section of boron for 
thermal neutrons, the values given lying between 
10-* and 700X10-**. The value will also, 
of course, depend on the average effective energy 
assumed for the thermal neutrons. We assumed 
op =600X10-*. Any change in the cross section 
of boron would evidently affect all our data by 
the same factor. A measurement of the cross 
section of Cd performed as a check on the method 
gave the correct value. 

Table I includes all of our results. Column 3 
indicates the number of hydrogen atoms per 
atom of the element in the solution. Column 4 
gives the ratio J)/I of the number of neutrons 
with water and with the solution, and the sta- 
tistical error. Column 5 gives the same ratio, 


*H. H. Goldsmith and F. Rasetti, Phys. Rev. 50, 328 


(1936). 

‘J. R. Denning, G. B. m, G. A. Fink and D. P. 
Mitchell, Phys. Rev. 48, 268 (1935); E. Amaldi and E. 
Fermi, Phys. Rev. 50, 899 (1936); G. A. Fink, Phys. Rev. 
50, 738 (1936); J. G. Hoffman and M. S. Livingston, 
Phys. Rev. 52, 1228 (1937). 
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TABLE I. Capture cross sections for thermal neutrons, 


ELE- Compo 
5B |H:BO; 1,000 | 2.02 +0.02 | 2.02 | 
7N |NH.NO; 24.3] 1.124+0.006] 1.086| 4. 
9F |NaF 105 | 0.998+0.006 0.990| <9'5 
11 Na |NaF 105 | 0.998 40.006] 0.990 | <9’5 
15P |HsPO, 43.5] 1.088+0.007/ 1.040] 
17Cl 145 | 1.32 40.01 | 1.32 | 97 
19K |KF 82 | 1.067+0.005/ 1.067] 
23V_ |VOSO, 190 | 1.094+0.006] 1.061| 68 
25 Mn|MnSO, 161 | 1.117+0.005] 1.099| 94 
27 Co |CoSO, 470 | 1.12 +0.01 |1.12 | 33 
33 As |As:05 75 | 1.187+0.006| 1.147} “65 
34Se |H2SeO; 75 | 1.283+0.007] 1.261] 14.5 
35 Br |NaBr 104 | 1.14 40.01 |} 1.116] 7 
48 Cd |Cd(NOs)2 | 10,000 | 1.49 +0.01 | 1.49 |2909 
531 |Nal 115 | 1.133+0.006) 1.101] 68 
56 Ba |Ba(C2H302)2} 200 | 1.001+0.007| 0.985| <j 


corrected for the change in hydrogen density as 
above explained. The last column gives the 
capture cross sections in 10~** cm*. The value of 
o was calculated by means of the formula 


Cu Io 
1). 
c I 


When the element investigated possesses only 
one isotope, no doubt exists about the process 
responsible for the absorption. Some other cases 
may warrant a few remarks. 

Cl—The weakness of the 37-minute activity 
of Cl** indicates that the high absorption is 
mainly due to capture by Cl*, with formation of 
long-lived Cl**, 

K—By a similar argument, most of the ab- 
sorption appears to be due to capture by K® 
with formation of K*°. 

Se—No strong activation corresponds to the 
large capture cross section, which must therefore 
correspond to formation of stable isotopes. 

Br—Both capture by Br” and Br* will be 
important, since all activities are strong. 

As it can be seen from the statistical errors 
reported in column 4 of the Table I, the values 
of the cross sections may be considered as 
accurate within 10 percent or better for the 
higher values, and within about 25 percent for 


the smallest. 
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Self-Consistent Fields and Diamagnetic Susceptibility for Magnesium III* 


ns, 
W. Jacoue Yostt 
i os Brown University, Providence, Rhode Island 
oO (Received June 24, 1940) 
1.2 
<0.5 Self-consistent fields, both with and without exchange, have been calculated for doubly ionized 
= magnesium in the normal state. All of the equations have been solved numerically with the 
n° aid of an electrical calculating machine. The Z» functions have been tabulated along with the 
3.2 corresponding P functions in each case. The former functions do not differ by more than 
6.8 +0.003, at any point, in the two last stages of the work. The diamagnetic susceptibility has 
By been calculated, yielding the value: —3.74X10~, 
6.5 
* INTRODUCTION methods have been employed, together with a 
N the theoretical investigation of the metallic comparative analysis Gwe 
state of a particular metal, it is desirable to METHOD OF SOLUTION 
ity as have the radial wave functions for the free ion. Both Hartree! and Fock? have discussed the 
s the The self-consistent field calculations yield the details of the derivation of the self-consistent 
lue of best set of such functions available by present field equations with the use of the variation 
methods. In order to make available a suitable principle. We shall refer to the equations in 
set of radial functions for the double ion of which the exchange terms occur as the Fock 
magnesium, the self-consistent fields for the equations, those in which they are neglected as 
1s°2s?2p° configuration have been calculated, the Hartree equations. Basically, the difference 
both with and without the inclusion of the lies in the inclusion of the symmetry demands of 
only ch the Pauli exclusion principle in the Fock equa- 
rocess le off tions. In addition, the Hartree equations result 
Cases relatively JA from the assumption of spherical charge distribu- 
of investigating the question raised by Hartree' tion, whereas higher order Coulombic terms are 
tivity concerning the relative merits of two methods included in the Fock equations. With these 
on is suggested for solving the so-called Fock equa- points in mind, the resulting equations for the 
ion of tions. The present paper represents the results radial functions for any ml shell can be written, 
of the calculation for Mg III in which both in the Hartree case, as: 
e ab- 
y P" (nl|r)+(V(nl|r) — (nl) JP(nl|r) =0; (1) 
o the while the corresponding Fock equation is: 
refore P" (nl |r) Ve(nl|r) — ]P(nl|r) +T(nl|r) =0. (2) 
‘Il be In these equations, and throughout the remainder of this work, atomic units are used. In the case 
of Mg III, there are three equations of each type to be solved. If two sets of functions, defined by: 
errors 
values Z,(nl, n'l’ |r) = f P(nl| p)P(n'l’| p) p*dp, 
ad as (3) 
r the as 


the degree of Doctor of Ph 
+ Now at Syracuse University, Syracuse, New York. 


*V. Fock, Zeits. f. Physik 62, 795 (1930). 
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1D. R. Hartree and W. Hartree, Proc. Roy. Soc. A150, 9 (1935). 
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are introduced, the various V(r) and 7(r) functions can be written: 
V(1s) = Vg(1s) =24/r— Wo(1s 1s) —2Wo(2s 2s) —6Wo(2p 25), 
V(2s) = Ve(2s) =24/r—2Wo(1s 1s) — Wo(2s 2s) —6Wo(2p 25), 
V(2p) = 24/r—2Wo(1s 1s) —2Wo(2s 2s) —S5Wo(2p 2p), 


Ve(2p) = V(2p) +2r°W2(2p 29), 


(4) 


T(1s) =[Wo(1s 2s) — 2s) ]P(2s) +3rW,(1s 2p)P(2p), 
T (2s) =[Wo(1s 2s) — e(1s 2s) ]P(1s) +3rW,(2s 2p) P(2p), 
T(2p) =rWi(1s 2p) P(1s) +rWi(2s 2p)P(2s). 


The W, functions are related to the Y; functions, used by Hartree, by the relation: 


2 \ Y,(nl, n’l’ |r) 


W,(nl, n'l’ |r) -( 


(5) 


as can be verified by a partial integration of the defining Eq. (3). 


The method usually adopted for solving these 
equations is one of successive approximations, 
in which the solutions of the Hartree equations 
are used to form the first approximation of the 
Fock solutions. In this method, it is customary 
to continue the work until the maximum change 
in each Zy function in going from one approxi- 
mation to the next becomes smaller than some 
arbitrarily prescribed amount. 

An alternate method has been used by 
Torrance’ in the case of the Hartree equations 
for carbon. The principle of this method is to 
solve the equations as a simultaneous set of 
integro-differential equations, subject to a group 
of boundary conditions. This method will be 
referred to as the method of direct solution. 
Hartree! has suggested the use of such a method 
for the Fock equations and has discussed the 
advantages and disadvantages, in principle, of 
such a method. In the present work, both 
methods have been used in an attempt to 
determine their relative merits as judged by the 
actual computational work involved. The nu- 
merical methods used were the same as those 
discussed by Mooney.‘ These are the most 
valuable whenever the computational work is 
all carried out with the aid of an electric calcu- 
lating machine. On the basis of the present 
investigation, the method of successive approxi- 
mations (s.a.) is definitely superior to that of 
direct solution (d.s.) for both the Hartree and 
the Fock equations. The reasons for this opinion 


*C. C. Torrance, Phys. Rev. 46, 388 (1934). 
4R. L. Mooney, Phys. Rev. 55, 557 (1939). 


are: (1) The d.s. prohibits the inward integration 
of the y-equation, which is of great help in the 
Hartree equations; (2) in d.s., it is necessary to 
estimate each W, function at the origin. This 
introduces several additional parameters which 
are difficult to estimate; (3) the total number 
of trials necessary to obtain the correct ¢ 
parameters is less in the s.a. than in the ds, 
method ; (4) in the d.s., each W; function must 
be calculated for all of the solutions attempted, 
whether they satisfy the boundary conditions 
or not. This last objection is probably the most 
serious since there are seven W, functions in the 
case of Mg III and these had to be calculated 
for sixty-two solutions which were of no further 
value because they did not satisfy the boundary 
conditions. 

It should be emphasized that the boundary 
conditions present a far greater problem in the 
Fock equations than in the Hartree equations, 


0.2 04 06 O8 10 20 3.0 40 50 


Fic. 1. The difference, AP, between the Hartree and 
Fock solutions. 7, represents the position of the maximum 
of the respective function. 
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SELF-CONSISTENT FIELDS 


TABLE I. Radial wave functions. 


EXCHANGE 
P(is) P(2s) P(2p) 


WitTHouT EXCHANGE 
P(is) P(2s) P(2p) 


eessss 


This is due to the fact that the latter are homo- 
geneous, whereas the Fock equations are non- 
homogeneous. Consequently, two aids generally 
used in the Hartree solutions are not possible 
with the Fock solutions. These are the use of 
inward integration of the so-called 7-equation, 
and the normalization of solutions by multi- 
plication of unnormalized solutions with the 
appropriate constants. 


RESULTS OF THE SOLUTION 


In Table I, the three normalized P functions 
are tabulated for both the Hartree and the Fock 


Fic. 2. Total radial charge densi 
from the various 


and the contribution 
shells. 


equations. Three digits to the right of the 
decimal are given, although six were retained 
throughout the computational work. Also, only 
half of the actual tabulations are given.* At the 
end of this table, the various e-parameters are 
listed. The curves in Fig. 1 show a marked 
resemblance to similar graphs given by Hartree® 
for argon-like configurations. This suggests the 
use of such graphs in order to obtain approximate 


TABLE I. (Continued) 


WitHout EXCHANGE 
P(is) P(2s) P(2p) 


—0.791 741 
—0.743 
— 0.697 
— 0.652 


— 0.567 
—0.491 
— 0.423 
— 0.363 
—0.310 
— 0.264 
—0,224 


—0.161 
—0.114 
— 0.080 
— 0.056 


— 0.030 
—0.012 
— 0.005 
— 0.002 
— 0.001 


With EXCHANGE 
P(is) P(2s) P(2p) 


001 —0.756 .728 
001 —-0.705 . 
— 0.656 
— 0.609 


—0.522 
— 0.445 
— 0.377. 
—0.318 
— 0.267 
—0.224 
— 0.187 


—0.129 
— 0.089 
— 0.060 
— 0.041 


—0.023 
—0.010 
— 0.003 
— 0.001 
— 0.000 


003 
‘001 
5.210 


99.692 7.426 99.469 8.944 


* More complete tables, which are essential for certain 
ee can be obtained upon application to the 
author. 

*D. R. Hartree and W. Hartree, Proc. Roy. Soc. A166, 
450 (1938). 
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0.379 0.093 0.001 | 0.378 0.095 0.001 
00150 | 1.008  .246 .010/1.006 011 
0.0200 | 1.266 .309 .018 | 1.264 317 .019 
90250 |1491 362 .027/1488 .372 .028 \ 
0.0300 | 1.686  .038 | 1.683 419 ‘ 
0.0350 | 1.853 .050/1.850 458 .052 A 
9.0400 | 1.996 .063/1.993  .066 4007 ~S 
0.0450 | 2.116  .503 .078| 2.113 515.081 
0.0500 | 2.216 2.213 534.097 ~_ 
(5) 2.363 .546 .127/2.360 .132 ad 04 O8 r—» 12 16 
2451 .552 .163|2447 .170 
2493 .520 .243/2489 .252 
2465 488 .284|2461 487 .295 
— 2413 447 .326/2.410 442 .339 
ation 2343 .399 .369/2340 .382 
2261 .346 .411/2.258 333 .426 
n the 2168 .453|2.166 271 
ry to 2069 .227 .495/2.068  .206 .511 
This 1.967 .535/1.966 
which 1.758 .614} 1.758 
1.552 —0.098 .687| 1.553 —0.138 .707 
imber 1.358 —0.228 .755| 1.360 —0.273 .776 
ct 1.179 —0.351 .817| 1.183 —0.402 
1.017 —0.468 .873| 1.022 —0.522 
e ds. 873 —0.576 .923| .878 —0.633 .945 
must 0.30 | .746 —0.676 .967| .751 —0.735 .990 
pted 032 | —0.766 1.005} .641 —0.826 1.028 
040 | 322 —1.033 1.104] .330 —1.090 1.127 
044 | .227 —1.116 1.126] .234 -1.169 1.149 
in the 048 | .158 —1.171 1.135] .166 —1.217 1.157 = a 
052 | 110 —1.201 1.131] —1.241 1.152 001 
ulated 036 | 076 —1.212 1.119] (083 1.245 1.138 ‘Ol 
urther 0.60 053 —1.206 1.099} .058 —1.231 1.116 ‘001 
0.64 | 036 1.188 1.072] .040 —1.205 1.087 
ndary 0.68 | .025 —1.160 1.042] 1028 —1.169 1.054 
0.72 | 017 —1.124 1.008] .019 —1.126 1.017 566 541 
0.76 | 012 —1.083 .972| .012 —1.078 505 
ndary 030 | [008 —1.037 ‘007 —1.026 937 = an 
in the 0.84 | 005 —0.989 .896] .005 —0.972 .895 = 
ti 0.88 | .004 —0.940 .857| .003 -—0.917 350 318 
0.92 | .003 —0.890 .818| .002 —0.863 .811 309 396 
0.96 | .002 —0.840 .779| .001 —0.809 .769 "261 
1 .208 179 
2 159 132 
2 -120 097 
2 091 071 
| 2 054 041 
3 026 018 
3 012 008 
3 .003 
4 001 
4 001 
5.990 
_ 
ree and 


£23 


= 
= = = 


Fock solutions whenever the corresponding 
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Hartree solutions are available. The various Z» 
functions can be calculated from the proper P 
functions. These functions represent the total 
charge within a sphere of radius r, due to a 
single electron. The Z functions do not vary by 
more than +0.003 at any point in the two last 
stages of the work for either the Hartree or the 


solutions. 


The radial charge distribution due to a 
particular ml shell is given by 


p(nb|r) (6) 


in electrons per unit r interval. These functions 
for the various shells in Mg III are plotted in 
Fig. 2. In order to illustrate more clearly the 
effect of the exchange terms, the change in the 
total radial charge distribution has been plotted 
in Fig. 3. The general result is a shrinking of 
the ion when the exchange terms are introduced. 
This gives an increase in the charge density near 


where 7 is expressed in multiples of ao, the 
atomic unit of length. Table II gives the result 
of the calculations for Mg III based on the 
solutions of the Fock equations. The table also 
includes results based upon experimental work. 
The earlier result quoted by Joos is believed to 
be in error because it cannot be correlated with 


TABLE II. Calculated susceptibilities. 


—x 
Contribution from 1s shell 0.04 
Contribution from 2s shell 0.90 
Contribution from 2) shell 2.80 
Total calculated value 3.74 
Result given by Joos! 


Result given by Pauling? 


2(2/+1)r2(nl) =10-*X 0.790 2(21+1) f P?(nl|r)r°dr, (7 
nl n 0 


1G. Feet Zeits. f. Physik 19, 347 (1923). 
?L. Pauling, Proc. Roy. Soc. All4, 204 (1927). 


0.4 


Fic. 3. The difference, Ap, in the total radial 


density as calculated from the Hartree and the Foct 
solutions plotted as a function of electron-nuclear gj, 


tance, 7. 


its principal maximum, with a corresponding 
decrease at large r values. 


DIAMAGNETIC SUSCEPTIBILITY 
The molar diamagnetic susceptibility of an 
atomic or ionic species in the free state is give 
by 


other data. However, both of the last two valus 
are estimates based on experimental data fy 
various salts. Pauling gives several reasons fy 
believing the value he quotes is more reliabk 
The present calculation seems to substantia 


Pauling’s value. 
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Thermal Conductivity of Metals 
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As indicated in a previous paper, elementary classical considerations suggested a relation 
between thermal conductivity k, atomic heat C, density p, and absolute temperature T of the 
form k/Cp=K/T+X’. This is now shown to be in agreement with recent theories on the metallic 
state. The first term on the right relates to conduction by the crystal lattice, the second term, 
K’, refers to conduction by electrons. It is found experimentally that K’ is the conductivity in 
the molten condition. The present paper shows these relations for lead, tin, and zinc. The data 
reported for lead are new. A modification of the experimental method previously used for zinc 
was required for the work on lead (a much poorer conductor) and is here explained. 


INTRODUCTION 


N a previous paper’ an attempt was made to 
I explain the thermal conductivity of metals 
and its change with temperature on the basis of 
a transfer of heat partly by elastic waves in the 
crystal lattice and their rapid absorption and 
partly by atomic collisions as in a gas. The 
former introduced a temperature coefficient, the 
latter did not, the temperature coefficient being 
associated with the velocity of the elastic waves. 
Elementary classical considerations led to an 


equation of the type 
k/pC=K/T+K’ (1) 


where k is the thermal conductivity; C, the 
atomic heat; p, the density; JT, the absolute tem- 
perature, and K and K’ constants. The theory 
suggested that the intercept should be the value 
of k/pC for the liquid state. The present paper 
offers data on lead, tin and zinc confirming this 
suggestion and interprets the equation in the 
light of recent theory. 


THEORY 


A. H. Wilson? from the standpoint of modern 
theory shows that thermal conductivity should 
require two terms, one due to the lattice, which 
as given by classical physics is k,=4pClu, and 
one due to electrons, which on the basis of 
modern theory is 

nrk?T 


m 


'C. C. Bidwell, Phys. Rev. 32, 311-314 (1928). 
* A. H. Wilson, Semi-Conductors and Metals (Macmillan, 
New York, 1939). 


l is the mean free path in the lattice; u, the 
velocity of sound in the lattice; m, the number 
of free electrons per unit volume; k, the molec- 
ular gas constant; 7, the temperature; m, the 
mass of the electron ; and r, the “‘relaxation time,” 
the average time between collisions. 

We may write k=k,+k, (although these are 
not independent, since the electrons afford a 
mechanism for scattering the lattice waves). 


Hence 
nrk®?T 


k= (2) 


m 

According to Wilson, k, becomes constant at 
high temperatures, while k, reaches a maximum 
and then decreases as 7-'. The experimental 
data on lead, tin and zinc here reported do not 
yield straight lines when k is plotted against 
1/T; but straight lines are obtained if k/pC is 
plotted against 1/7. The theoretical equation 
in the form 

nrk?T 


(3) 
mpC 


must then be equivalent to 


k 1 
—=K—+K’. 
pC 


(1) 


According to this r « 7-! and lu« T-. The Inter- 
national Critical Tables give data on the velocity 
of sound for platinum and copper at 20° C, 100°C 
and 200°C. For these data u« 7-'. According to 
Wilson, / depends on the grain size and is there- 
fore only slightly affected by the temperature. 
Thus with p and C approximately compensating 
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in their temperature effects, Eqs. (3) and (1) are 


in agreement. 
The intercept term must include the atomic 


as well as electronic conductivity for the liquid 
metal. If the atomic conductivity is that of an 
insulating liquid, it will be only one or two 
percent of the electron conductivity. Such con- 
ductivity increases with temperature and may 
thus account for the slight positive slope experi- 
mentally observed. 


The experimental evidence 


The data on zinc previously reported* were 
obtained by a new experimental method and 
agree closely with values found by Konno,‘ the 
only other data published for molten zinc. The 
data here presented on tin are Konno’s, no work 


*C. C. Bidwell, Phys. Rev. 56, 594 (1939). 
*S. Konno, Phil. Mag. 40, 542 (1920). 


Fic. 1. Experimental data for lead at various temperatures. 


on tin having been done by the author. For 
lead, however, Konno’s values for the molten 
state are far too low to give value of k/Cp in 
agreement with the intercept of the equation 
for the solid metal. Hence new observations 
were made on lead with the newly developed 
experimental method. These observations are 
shown graphically in Fig. 1. Figure 2 shows the 
author’s values for the thermal conductivity of 
lead together with values obtained by Konno! 
and Van Dusen and Shelton.’ Konno’s values 
for the molten state are about 30 percent lower 
than the newer values here reported. Figures 3 
and 4 show the course of the thermal conduc- 
tivity curves for zinc and tin. For all these plots 
specific heats and densities were obtained from 
the published tables and are averages of what 


5M. S. Van Dusen and S. M. Shelton, J. Research Nat. 
Bur. Stand. 12, 429 (1934). 


0.0007 


Fic, : 


to be the best values. With the author’s 
yalues for molten lead, the plots for tin, lead 
zinc are seen to be in agreement with the 
ted law, the intercepts giving the value of 

k/Cp for the molten state. 


MopIFICATION OF THE EXPERIMENTAL. METHOD 
IN ITS APPLICATION TO LEAD 


The experimental procedure as reported in a 

vious paper on zinc (reference 3) consisted, 
in brief, in establishing a flow of heat vertically 
down a specimen packed in Silocel, measuring 
the temperature at various levels (from which 
can be computed the gradient d7/dx and the 
change in gradient d*7/dx*) and measuring the 
lateral radial drop AT across the Silocel at the 
various levels. 

The procedure as reported in the work on 
zinc was improved by centering the graphite 
cylinder containing the test metal in an iron 
cylinder of about 12 cm inside diameter on 
which was wound a Nichrome heater ribbon. 
The space between the graphite cylinder (5 cm 
diameter) and the iron cylinder was filled with 
Silocel as was also a space between the iron 
cylinder and an outer asbestos cylinder of about 
20 cm diameter. The Nichrome windings on the 
iron cylinder were in parallel with similar 
winding on the graphite. By adjustment of 
currents in these two windings, the radial 
temperature drop across the Silocel could be 
reduced to small values or zero. The addition 
of the intermediate iron cylinder with its heater 
coil is found to be a great improvement. 


* 
LEAD 
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Fic. 2. Thermal conductivity of lead. & in cal. per sec. 
per cm per degree C. 
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Fic, 4. Thermal conductivity of tin. k in cal. per sec. 
per cm per degree C. 


A further modification of the method of 
handling the data in the case of such a relatively 
poor conductor as lead was found necessary. ‘In 
the method used, the lateral, radial heat loss 
between two levels is equated to the difference 
between the heat entering at the upper level and 
leaving at the lower level. The equation is 


r2 
= f k,2nrldT/dr. (5) 
"1 
This becomes 


2.3 log (r2/11). 


since (d7/dx,—dT/dx:)/l=@T/dx*. k, is the 
thermal conductivity of the test metal; &,, that 


(6) 
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of the Silocel; A», the cross section of the test 
metal; AJ, the radial temperature drop across 
the Silocel at the level in question ; 72, the outer 
radius of the Silocel (inside the iron cylinder) ; 
its inner radius.® 

In the case of molten metals it becomes 
necessary to add to the left side a term to cover 
the heat flow down the walls of the graphite- 
containing cylinder. The equation then becomes 


2rk, AT 


2.3 log (r2/r1) 
—hgrAgr(d?T/dx*) gr. (7) 


kA» (d?T /dx*) 


This may be written 
AT 
kzAm2.3 log (r2/r1) 
RerA gr(d?T /dx*) gr 


It is found for lead that values of d?7/dx* 
plotted against corresponding values of AT give 
straight lines in all cases. The intercept term is 
small and the variation of (d?7/dx*),, from top 
to bottom of the cylinder appears to be negligible. 
This is reasonable, at least when the contained 
metal is a comparatively poor conductor like 
lead, where the graphite has about five times the 
conductivity of the solid lead and about ten 
times that of molten lead. The comparatively 
high conductivity of the graphite means a 
correspondingly small temperature gradient in 
the graphite, changing at an approximately 
constant rate. In the case of zinc (see reference 3), 
the conductivity of the metal was about the 
same as that of the graphite and, for zinc, 
d@*?T/dx*? was found to be constant over a con- 
siderable range. This is in agreement with the 
above explanation of the constancy of (d?7'/dx*) gr. 


(d*T/dx?) »= 


(8) 


For temperatures — 130°, —35°, +35°, +97°, +250°, 
r,=2.6 cm; f2, 4.5 cm; Am = 13.83. For 392°, 620°, 1=2.5 
cm; fz, 5.65 cm; Am =13.83. Values of k, for various 
temperatures are shown in the previous paper. The lead 
——- was prepared by slow freezing from the bottom. 

is was effected by lowering the molten lot slowly 


through a furnace at the rate of about two cm per hour. 
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In the case of zinc it was necessary merely to 
add the cross section of the graphite to that g 
the zinc and use Eq. (7), AT and the slope of 
the d7/dx line both being constant. For lead 
however, d?7'/dx* and AT are not constant but, 
when @?7/dx* is plotted against AT, a strai t 
line results as indicated by Eq. (8). From the 
slope of this line the thermal conductivity g, jg 
obtained. No correction for the graphite walls jg 
required since the quantities having to do with 
the heat flow down the walls all appear in the 
intercept. 


PRECAUTIONS 


In the use of this method it is of the utmoy 
importance that stable unchanging temperature 
conditions be attained before starting the heg 
flow down the specimen. Four §” quartz tubes 
were inserted at the outer radius of the Siloce 
(just inside the iron cylinder) and two such tubes 
in the molten specimen (or frozen in the solid 
specimen). Only three junctions were used, the 
Chromel wires welded at a common point and 
the copper leads connected to binding posts 
One of the junctions was kept at the bottom of 
one of the quartz tubes in the metal, the second 
raised or lowered to the various levels in the 
other quartz tubes and temperature differentials 
between these junctions read on a galvanometer, 
The third junction was kept in ice and differ. 
entials between it and either of the others read 
on a potentiometer indicator gave actual 
temperatures. 

It is usually necessary to wait two or mor 
hours before unchanging conditions are attained. 
Readings up and down the specimen and acros 
the Silocel at all levels are taken and then the 
downward heat flow started by applying heat in 
the d.c. heater at the top of the specimen. 
Final readings may be taken after another hou 
and corrections made for the initial readings 
The fact that one gets smooth curves and 
straight lines as required by the equations is 10 
guarantee of the correctness of the result. Only 
from absolutely stable conditions can sats 
factory results be obtained. 
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Properties of Rochelle Salt. III’ 


Hans MUELLER 
Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received June 28, 1940) 


The interaction theory of Rochelle salt is extended by supplementing the expression for the 
free energy of the deformed and polarized crystal by a term proportional to the fourth power 
of the polarization in the a direction. This second approximation explains the dielectric and 
elastic saturation phenomena above the upper Curie point and the dielectric, elastic, piezo- 
electric and pyroelectric properties in the ferroelectric temperature range. All experimental 
results are compatible with the hypothesis that the clamped crystal has only one transition 
point at about 5° C, where its dielectric constant has a very high maximum. The two Curie 
points of the free crystal are due to the interaction between the elastic deformation and the 
electric polarization. Instead of sudden changes of the heat capacity the new theory predicts 
slight changes of the slope of the specific heat curve at the Curie points. The dispersion of the 
dielectric constant is explained without assuming a relaxation time. 


HE results of paper II suggest the following 

explanation of the anomalous properties of 
Rochelle salt : If the crystal could be clamped in 
such a way that all elastic deformations are 
suppressed it would show only a single transition 
at about 5°C. At this temperature the dielectric 
constant in the a direction of the clamped crystal 
has a sharp and very high maximum, i.e., the 
reciprocal susceptibility x: vanishes (see Fig. 
3, II). This transition is not a Curie point, be- 
cause below this temperature the clamped crystal 
is not spontaneously polarized, its dielectric 
behavior remains normal and x; increases again 
to reach a normal value at low temperatures. The 
transition appears to be similar to those ob- 
served in other crystals, e.g. HBr, HI.2 Most 
probably it is due to a modification of proton 
bonds or to hindered rotation of polar groups. 
The transition of the clamped Rochelle salt 
crystal does not involve changes of the crystal 
symmetry and it leaves the crystal structure 
essentially unaltered. It produces a change of the 
optical properties as is indicated by the fact 
that the temperature gradient of the bire- 
fringence ,—mn, is almost three times smaller 
below — 20°C than it is above 25°C (see Figs. 18 
and 19, I). 

'This is the third publication on this subject by the 
writer. The previous articles, Phys. Rev. 47, 175 (1935) 
and Phys. Rev. 57, 829 (1940), will be referred to as 
paper I and II, respectively. The numerals I and II are 
used also in the references to figures and equations given 
in these previous pa 


*C. P. Smyth and C. S. Hitchcock, J. Am. Chem. Soc. 
55, 1830 (1933) ;G. Damkohler, Ann. d. Physik 31, 76 (1938). 


In the experiments on freely deformable 
crystals the existence of this transition is masqued 
by a series of other effects. As was shown in 
paper II these effects result from the interaction 
between the elastic deformation and the electric 
polarization. They can be accounted for by the 
laws of classical crystal physics. Because of this 
interaction the free crystal has two Curie points, 
one above and one below the transition tempera- 
ture. Between the Curie points the free crystal 
is spontaneously polarized and deformed. The 
spontaneous deformation changes the symmetry 
of the crystal from rhombic hemiedric to mono- 
clinic hemimorphic with the a axis as polar 
axis.* Although the interaction theory makes no 
use of the concept of an inner field F=E+fP, 
it furnishes an equation, 1/«:=x1(1—fie/x1C4s) 
(from Eqs. (4), II), for the susceptibility «x; of 
the free crystal above the Curie point, which is 
analogous to that, 1/x,=(1—0/T)fT/6 (Eq. (9), 
I), given by the inner field theory. A comparison 
of these relations shows, that the interaction 
gives rise to an “apparent Lorentz factor” 
f=fie/cas=0.045, while the Curie temperature 
is 9@=f17/x:C44. The Curie points of the free 
crystal are the temperatures 7=6 at which the 
determinant Dy, = X1C44 = (T 
vanishes. c4=11.6X10" and fi,=7.2X10* are 
the true elastic and piezoelectric constants‘ for 
shearing strains y,. 


*H. von R. Jaffe, Phys. Rev. 51, 43 (1937). 

‘ The present investigation shows that the values of fix 
as derived in paper II from Norgorden’s or Mikhailov’s 
measurements are too small. 
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Fic. 1. Theoretical stress-strain relation for shear in the 
(b, c) plane of Rochelle salt above the Curie point. 


The interaction theory has served to explain 
the relations between the dielectric, piezoelectric 
and elastic properties of Rochelle salt above the 
upper Curie point. We will show here that it 
can account also for the ferroelectric anomalies 
at temperatures between the Curie points. To 
this purpose it is necessary to supplement the 
expression (1, II) for the free energy ® of the 
polarized and deformed crystal by terms of 
higher order. The origin of the new terms can be 
explained as follows: At temperatures near the 
transition point, where x; is very small, the 
P.(E,) relation for the clamped crystal must 
show saturation. Since the clamped crystal is 
orthorhombic this saturation effect introduces 
only terms in even powers of P, and we shall 
therefore use as a second approximation 


(1) 


where 4 is the quadratic form (1, II). Terms 
involving higher powers of the strains or of P, 
and P, or mixed terms, e.g. P,*y,’, are not likely 
to be of importance, because it was shown that 
the true elastic and piezoelectric properties of 
Rochelle salt are not unusual. A second set of 
correction terms arises from the change of the 
crystal symmetry at the Curie points. This fact 
would require the introduction of ten new terms 
involving four new elastic constants C14, C24, C34, 
Css, five new piezoelectric constants f11, fi2, f13, fos, 
fas and a new susceptibility x23. New experimental 
evidence, which will be presented in an other 
paper, shows that these correction terms of the 
second kind play an important role in the theory 
of the electro-optical and photoelastic effects, 
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the anomalous thermal expansion and the electyg, 
striction in Rochelle salt. However, these 

are small in comparison with }BP,‘, and 
they are not entirely negligible they play a 
secondary role in the following discussion of the 
dielectric, elastic and piezoelectric effects, 
second approximation, Eq. (1), furnishes the 
fundamental relations 


(2a) 
(2b) 


To test their validity we shall follow the pro. 
cedure adopted in paper I by discussing firg 
the saturation effects which occur above the 
Curie point, where Diy=x:1Cu—fi2>0, and 
where the experiments furnish more reliable 
data because no hysteresis effects occur. Theg 
data will furnish the value of B and it will they 
be possible to test the theory for the ferro. 
electric range. 


DIELECTRIC AND PIEZOELECTRIC SATURATION 
ABOVE THE CURIE POINT 


The laws for the dielectric properties and the 
inverse piezoelectric effect of the free crystal are 
found from Eqs. (2) by setting Y,=0. The direct 
piezoelectric effect is obtained by elimination of 
yz and by assuming £,=0. If one takes into 
account Eqs. ((4), II) these procedures furnish 


E,=P,/m+BP;'. (3a) 
(3b) 
(3c) 


The P.(E.), and P.(Y.) relations must 
therefore show saturation effects when 1/x; o 
1/di, are small, i.e., at temperatures slightly 
above the Curie point. At higher temperatures 
the cubic terms are negligible in comparison 
with the linear terms. This remark explains why 
it would be futile to search for saturation effects 
in the elastic or piezoelectric properties of the iso- 
lated crystal. In an isolated crystal the polariza- 
tion creates a depolarization field E,=—4P, 
and hence the coefficient of P, in Eq. (2b) is 
(4r+ x1) and is about 100 times larger than 
1 / Kl. 

The existence of saturation in the piezoelectric 
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PROPERTIES OF ROCHELLE SALT III 


effects has been reported by Schulwas-Sorokin® 
and Hinz,‘ but sufficient data for the determina- 
tion of B are available only from the dielectric 
measurements. The validity of Eq. (3a) was 
demonstrated by Kurtschatow’ for static, and 
by the writer for dynamic measurements by the 
method of Sawyer and Towers (Fig. 10, I). 
This equation is equivalent with the relation 
E=(1—0/T)F+86fF* (Eq. (7), 1) of the inner 
field theory when the approximation F=fP is 
valid. It leads directly, without resorting to 
approximations as was necessary in paper I, to 
the experimentally verified relation (see Fig. 8, I 
-and Eq. (19), 1) 

E,/X*=M(t-t) + NX, (4) 


where X = and where xe=0P,/0E, 
is the reversible susceptibility. M=1/C(3B)}, 
N=1/3(3B)* and C is the Curie constant. For 
small fields Eq. (4) reduces to 

kg=k,(t—GE,’), G+=g(t—t,), (4’) 


where g=1/C(3B)*. It was shown in paper I 
that the observations verify the relations 
gi=M?/C, N=CM/3, as demanded by the 
theory. The experimental values g= 1.05, V= 800, 
M=13.5, C=178 (near the upper Curie point) 
- furnish B=5.8+0.7X10-*. B takes the place 
of 6f* of the inner field theory. Its value may 
vary somewhat with temperature, but more 
accurate data would be required to establish 
this dependence. 


ELAsTic SATURATION IN THE 
FoILeD CRYSTAL 


Although the true elastic properties of Rochelle 
salt follow Hooke’s law, this law is not satisfied 
for shearing stresses Y, when the crystal is 
covered with a grounded foil. From Eqs. (2) 
and the Curie-Weiss law’ (10, II) for the elastic 
compliance sq one finds for the foiled crystal 
(E,=0) 


Y,= ( Y,+cuy,)[1 +( 


(t—te)/ocas. (5) 


*R. D. Schulwas-Sorokin, Zeits. f. Physik 77, 541 (1932). 

*H. Hinz, Zeits. f. Physik 111, 617 (1938). 

"I. V. Kurtschatow, Seignette electricity (Moscow, 1933), 
Fig. 15. The isothermal static values may differ from the 
adiabatic dynamic values. To be exact, one should differ- 
entiate between the free energy and the enthalpy. The 
existence of an electrocaloric effect indicates that the 
difference is not negligible. 


At the Curie point, =¢,, the relation is 
(S’) 


Figure 1 presents the resulting strain-stress 
curves for various temperatures t=/,+At. The 
values B=5.8X10-8, 10", 


fu=7.2X104, x1 =0.05+(t—#.)/178, 


_¢=66.7 X10-" were used in the calculation. We 


note that large deviations from Hooke’s law 
must be expected if the temperature of the foiled 
crystal is between 24° and 27°C. The elastic 
saturation should become apparent for shearing 
stresses of only a few kg/cm? and it should 
therefore not be difficult to verify this effect by 
ordinary elastic measurements. At higher tem- 
peratures the deviations from Hooke’s law vanish. 

A simpler and more accurate method for 
studying the elastic saturation in Rochelle salt 
is based on the fact that the resonance and 


Lesu. 


Fic. 2. Influence of a steady electric field on the reso- 
nance frequencies of a plate of Rochelle salt at various 
tem tures above the Curie point: (a) 24.1°, (b) 24.2°, 
(c) 24.5°, (d) 25.0°, (e) 26.8°, (f) 30°, (g) 18°. (f) and (g) 
refer to the lower resonance frequency. The dotted curve 
gives the change of an antiresonance frequency at 24.1°C. 
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Fic. 3. Verification of the limiting law for the change of 
the resonance frequencies in small steady fields. 


antiresonance frequencies of all those vibrations 
of a crystal plate, which can be excited by an 
alternating field in the a direction, must depend 
on the compliance coefficient s4,. For any reso- 
nance frequency f this dependence generally 
will be of the form f=(a+8sq)~!, or since 


(6) 


where the constants a and b depend on the mode 
of vibration and therefore on the size and shape 
of the plate and on the other elastic compliance 
coefficients s;., which are to be considered as 
temperature independent. If the plate is de- 
formed by a constant strain y,° the values of a 
and 8 remain practically unchanged, but if the 
electrodes are in close electrical contact with the 
crystal the compliance s4, must be replaced by 
the “reversible compliance” s4.°= —dy./0Y:, 
which is given by the slope of the y.( Y.) curves 
in Fig. 1 and depends on y,°. Since a constant 
strain cannot be produced by mechanical 
means without altering the boundary conditions 
for the vibrations, it is more convenient to 
produce it by the application of a steady electric 


field Eo in the a direction. A simple calculation 
shows that a field Ey changes the compliance 
from S44 to S44”, where 


1/sas® =Cas—fis?/ (xi +3BP,°), 


and where P» is the polarization due to E, in Re 
free crystal. Therefore, according to Eq. (3a) 
where xp is the reversible 
susceptibility in the field Eo. By combining these 
relations it is found that the steady field g, 
changes the resonarice frequency from f to 
and the calculation 
gives 


(7) 
For small fields E, the approximation (41) 


furnishes 
Af /bf? = 


The change of the resonance frequencies of a 
foiled Rochelle salt plate by a steady field has 
been reported in a previous note.* Figure 2 gives 
the measurements on two resonance frequencies 
of a second crystal of approximately the same 
dimensions as the first, 3X5 cm in area and 
0.096 cm thick in the a direction.* Table I gives 
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Fic. 4. Verification of the law omens the change of the 
resonance frequency to the change of the susceptibility of 
Rochelle salt in a steady electric field. The dotted curves 
result from the elastic measurements, the full curves for 
24.1°, 25° and 26.8° are calculated from the dielectric 


data. 
8H. Mueller, Phys. Rev. 57, 842 (1940). 
® The crystal plates were supplied by the Brush Develop 


ment Company, Cleveland, 
trodes attached. The orientation of the edges could not 


be ascertained without damaging the electrodes. 


hio, with elec- 
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I. Temperature Cpcatne of two resonance 
of a Rochelle plate. 


pc 25.0 26.8 X10 


9 291 %3i2 345 40.6, 3.38 10.1 
IRC 49.0 50.0 54.0 59.7 71.25 103 3.57 


the temperature dependence of these frequencies 
when no steady field is acting. These data satisfy 
the law (6), with = 23.0°C. Thus the constants 
a and 6 can be determined in spite of the fact 
that the modes of vibration are not known. 

Figure 3 is a plot of 2Af/bf* versus the square 
of the steady voltage. In accordance with 
Eq. (7!) both frequencies lead to identical 
curves with initial slopes which vary with the 
temperatures in the ratios (1/1.1)* : (1/1.2)*: 
(1/1.5)* : 2* : (3.8)* as indicated by the dotted 
lines. In Fig. 4 the variation of A(1/f*)/d is 
compared with the variation of (1 —«g/x:)/(t—te) 
as deduced from the dielectric data given in 
Fig. 7, 1. This verification of Eq. (7) involves no 
adjustable parameter. The agreement signifies 
that the fundamental equations are valid and 
that the elastic data lead to the same value of B 
as the one derived from the dielectric meas- 
urements. 


THE FERROELECTRIC PROPERTIES 
OF ROCHELLE SALT 


At temperatures near the transition point x: 
is small, and the determinant Dyu=xiCua—fi? 
is therefore negative between the two Curie 
points. From Eqs. (2) it follows that in this 
temperature range the state of minimum free 
energy of a free and foiled crystal (E,=0; Y,.=0) 
has a spontaneous polarization P,° and a 
spontaneous deformation y,° 
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The spontaneous polarization can be directed 
along the +a or —a axis. Since it is associated 
with a spontaneous shearing deformation, all 
parts of a crystal will usually become polarized 
in the same direction, because, for purely geo- 
metrical reasons, different parts cannot be 
sheared in opposite directions without creating 
internal stresses. A small crystal, when cooled 
slowly, acts therefore as a single Weiss region and 
is pyroelectric, but in a large crystal it is difficult 
to prevent twinning, and frequently the crystals 
crack when they are cooled too rapidly. 

The influence of the spontaneous polarization 
and deformation on the dielectric, elastic and 
piezoelectric properties of the crystal is analo- 
gous to that resulting from the application of a 
steady field above the Curie point. The dielectric 
and piezoelectric constants are decreased, and the 
resonance frequency of a foiled plate increases. 
To arrive at a theory of these effects one takes 
into account that a small stress 5Y, and a small 
field 5E, change the polarization from P,° to 
P,°+6P,, and the strain from y,° to y,°+y,. If 
one neglects all terms containing higher powers 
of 6P, and dy, one finds from Eqs. (2) and (8) 


(9a) 
(9b) 
x1°=x1— /Cas. (9c) 


xi° is the reciprocal of the susceptibility «.;° of 
the spontaneously deformed clamped crystal, 
i.e., of a crystal which has been cooled without 
stresses to below the upper Curie point but which 
is prevented to deform any further if the field 
6E, is applied. The solution of the linear equa- 
tions (9a, b) gives 


—dby.= $44°6 Y,+4,,°6E., (10a) 
(10b) 


TABLE II. Elastic, dielectric and piezoelectric data for Rochelle salt between the Curie points. 


f ke «i? PY (sa —1/cas) X 1012 die X108 fu X10~4 BxX108 
20 1.49 0.51 38 420 12.9 22.2 6.8 7.4 
15 1.62 0.42 20 620 9.1 13.4 8.0 6.5 
10 1.68 0.38 15 720 7.6 10.7 8.5 6.4 
5 1.69, 0.37 15 740 7.4 10.4 8.4 6.1 
0 1.70 0.37 17.5 740 7.4 11.3 7.9 5.2 
-5 1.66; 0.40 23 700 8.3 13.7 7.4 4.4 
—10 1.57 0.47 32 600 10.9 18.6 7.2 4.3 
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Fic. 5. Verification of the relations between the spon- 
taneous polarization, the elastic, piezoelectric and dielectric 
properties of Rochelle salt between the Curie points and 
of their relations to the corresponding properties above 
the upper Curie point. 


where 
—X1/2D 1443/2044, (11a) 
= —f4/2D 4. (11b) 
—C44/2D 14. (11c) 


Saa° is the compliance coefficient of the foiled 
crystal, dy4° is the piezoelectric constant and 
x: is the susceptibility of the free crystal in the 
ferroelectric state. Equations (10) hold only for 
small stresses and weak electric fields. 

If one assumes that in a small. temperature 
range the variation of D,, can be approximated 
by a linear law Dig=(t—t.)€44/C Eqs. (11) can 
be cast into the form of Curie-Weiss laws: 


=0/(t-—t); o=C(fis/cas)?, 
A=Cfhis/cas, 
(P.°)?=h(t.—1); 


Since the ‘Curie constants’ vary with the 
temperature, we shall not use these relations, 
though they give a satisfactory representation 
of the observations near the Curie points.'° 
Instead we follow the procedure adopted in 
paper II, i.e., we calculate the fundamental 


h=1/CB. 


1° Paper I, ra.” and O. Norgorden, Phys. Rev. 49, 


820 (1936) ; 50, (1936). 
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quantities C4, fis, x1 and B from the available 


data. 

Equations (9) and (10) have the same form as 
the corresponding Eqs. ((2), II) and ((3), II) for 
temperatures above the Curie point. Hence we 
can use all the relations of paper II if we syp. 
stitute Sea, dia, «1 and x1 by dix”, «1° and y,0 
respectively. Thus it follows from Eq. ((8), II) 
that the compliance of the isolated crystal js! 


= S44? — 40rd 14"/ 


and the resonance frequency f* for longitudinal 
vibrations at 45° to the b and ¢ axes is 


and has no anomaly at the Curie points. The 
resonance frequency of the foiled crystal jg 
and it is easy to 
prove that Mason’s” relation is still valid. The 
electromechanical coupling coefficient is 


(13) 


and from the equations for f, f‘ and (11a) it 
follows that 


(14) 
According to (11b) and (11c) 
dy4°/ (15) 


This ratio must have the same value as d44/x,, 
because Eq. ((16’), II) is valid also below the 
Curie point and Mikhailov'* has found that 
5f/5D is very nearly constant at all tempera- 
tures. Eqs. (11) furnish also 


Xi (16) 


This relation justifies the procedure used in 
Fig. 3, II for the construction of the x:(¢) curve. 
It should be noted that the susceptibility «,;' of 
the spontaneously deformed clamped crystal 
differs from the susceptibility «..=1/x1 of the 
crystal in which all deformations are suppressed. 
While «.; has a single and very high maximum at 
the transition point near 5°C, x.;° has at this 


1 The approximation e¢°=47x;° is justified also in the 
ferroelectric range, because the dielectric constant 
nowhere smaller than 150. 

12 W. P. Mason, Phys. Rev. 55, 775 (1939). 

8G, Mikhailov, Tech. Physik USSR 4, 461 (1937). 
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temperature a minimum Kez has 
two maxima at the Curie points, where xei°= kei. 
Table II gives the experimental data on 
Rochelle salt at temperatures between the 
Curie points. f and k are from Mason’s” work. 
The values for «:° and are Bradford’s" 
measurements. From these data (S4a°—1/C4,), 
dy® and fis were calculated by using Eqs. (14), 
(13) and (15), respectively. The value of B is 
obtained from Eqs. (8) and (11c), which give 


B=1/2x,°P.™. (17) 


Within the rather large experimental errors 
fu and B have the same values as were found 
from the data above the Curie point. Both con- 
stants seem to decrease somewhat with tem- 


perature. 
The crucial test of the theory is carried out in 


Fig. 5. According to Eqs. (8), (11b), (11c) and 


(14) the quantities (P. 1 / 
and 1/x;° are all proportional to Dy, and should 
therefore have the same temperature variation, 
provided that Cas, fis and B are constant. Further- 
more, the last three curves should be continua- 
tions of the curves for 1/(s44—1/Cas), 1/dig and 
1/x;, as given in paper II for temperatures above 


the Curie point. To illustrate the agreement the © 


scales for the last three curves were adjusted so 
that the straight lines above the Curie point have 
the same slope. This procedure leads indeed to 
the result that the three curves become very 
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nearly identical in the temperature range be- 
tween the Curie points. They also coincide with 
the P,” curve, for which the scale was chosen 
so as to give the same minimum at 5°C. It 
would not be difficult to interpolate the results 
on both sides of the upper Curie point by smooth 
curves, and to ascribe the small deviations to 
experimental errors, which, due tc the gap 
effects, are particularily large near the Curie 
points. However, the elastic data show quite 
clearly that the curves have a small disconti- 
nuity at the Curie point. In a later paper we can 
show that such anomalies are to be expected. 
They are the result of the neglected correction 
terms of the second kind and they can be ex- 
plained as due to small changes of the values of 
the constants c44, fi4 and B at the Curie points. 
It is therefore not possible, on the basis of the 
present theory, to predict the ferroelectric 
properties near the Curie point from the measure- 
ments above this temperature, as was attempted 
with moderate success in paper I. However, 
since the changes of ca, fis and B are small, 
our theory is adequate in predicting the mini- 
mum or maximum values of the various quanti- 
ties at the transition point. At this temperature 
x1 is very small, and the assumption x:=0 gives 
(Dis) min= —fie. By using Eqs. (8) and (11) and 
the values of fis, C44 and B as derived from the 
measurements above the Curie point, one finds 
the following approximate values: 


calculated observed 
(Pz) max =fia(CasB) = 860 740 e.s.u. 
(4X 10-2 = 1.85 X 105 1.71 10° 
(«1°) min =Ca4/2f ig? =11.4 12.5* 
(d14°) min = 1/2f =7X10-* 10.4 10-* 


With the exception of (y.°)mex, Which has been 
measured for only one crystal, the deviations 


“E. B. Bradford, B.S. Thesis, Massachusetts Institute 
of Technology, 1934. P,° is obtained from the height of 
the hysteresis loops. Bradford’s values of P,° agree with 
those given by J. Habliitzel, Helv. Phys. Acta 12, 489 
(1939), but they are larger than those resulting from the 
measurements in paper I. Bradford's values 
or the at 1060 cycles/sec., are slightly 
higher than Habliitzel’s, but they are lower than Mason’s 


ues. The data in Table II are in e.s.u. 


from the observed values indicate that at the 
transition point x1>0. The experimental evi- 
dence, resulting from the three independent 
investigations of Bradford, Habliitzel and Mason, 
supports therefore the conclusion that the 


146 The elastic data are more accurate than the dielectric 
| gy Figure 5 is based on the original data of 
. Mason. 


1* J. Habliitzel, Helv. Phys. Acta 12, 489 (1939). 
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clamped crystal has no Curie point. It possesses 
a transition point where the dielectric constant 
reaches an extremely high maximum. 

The spontaneous polarization and deforma- 
tion lower the free energy of the crystal by 


(18) 


The maximum energy change is about 4500 
erg/cm?=0.017 cal./mol. Since, near the Curie 
point, (P.°)?=h(t.—?t), where h is about 5X10‘, 
the theory does not demand a sudden change of 
the specific heat at the Curie points, but only a 
change of the slope of the specific heat curve. 
This change of slope is only about 3X10~‘ 
cal./mol/°C * and is much too small to be de- 
tected by the usual methods. For the crystals of 
H:KPO, and H2KAsO,, for which our theory 
should also be applicable, a large and sudden 
change of the specific heat has been observed,'’ 
but this is not necessarily in contradiction with 
the theory, because one might expect such a 
change to occur at the transition temperature, 
i.e., at a temperature somewhat lower than the 
Curie point. The experiments leave some doubt 
whether the temperature, where the specific heat 
changes, is identical with that of the Curie point. 
In Rochelle salt no anomaly has been found, 
either at the Curie points or at the transition 
temperature. 

For large electric fields and large shearing 
stresses Y, Eqs. (2) predict hysteresis curves 
for the electric polarization of the free crystal, 
for the deformation y, of a foiled plate and for the 
direct and inverse piezoelectric effect. The dielec- 
tric hysteresis has been studied by numerous 
investigators. Bloomenthal'* and Hinz* have 
reported hysteresis of the inverse piezoelectric 
effect and the existence of elastic hysteresis 
follows from the recent work of the writer.® 
Contrary to the theory, the hysteresis loops have 
not a single Barkhausen jump and the coercive 
field is smaller than the theoretical values. 
These discrepancies are not surprising if one 
realizes that the theoretical Barkhausen jump 


17W. Bantle and P. Scherrer, Nature 143, 980 (1939); 
J. G. Hooley and C. C. Stevenson, Phys. Rev. 56, 121 
(1939); J. Mendelsohn and K. Mendelsohn, Nature 144, 


595 (1939). 
18S, Bloomenthal, Physics 4, 172 (1933). 
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corresponds to a transition from an instable to a 
stable equilibrium state. Structural inhomoge. 
neities, internal stresses or inhomogeneities of 
the electric field may be responsible for the 
transition’s occurring at smaller field intensities, 


DIELECTRIC DISPERSION IN 
ROCHELLE SALT 


Schulwas-Sorokin and Posnov'® and Gog. 
decke,® among others, have reported results 
which are assumed to arise from relaxation 
phenomena in Rochelle salt. The first two authors 
ascribe the relaxation to viscous resistance 
against mechanical deformations, the latter to 
the retardation time of the dipole orientation, 
It should, however, be noted that the measyre. 
ments were carried out in the ferroelectric 
temperature range and that the interpretations 
are based on the erroneous assumption that the 
dielectric and elastic properties of Rochelle 
salt follow the laws which hold for ordinary 
solids. A closer inspection reveals that the ob. 
served phenomena are due to the interaction, 
saturation and hysteresis effects, and that the 
existence of a relaxation time longer than 10+ 
sec. is not justified by any experiment. It seems 
rather obvious that power factor data and volt. 
age decay curves cannot be interpreted in terms 
of a relaxation time when most of the power logs 
is due to hysteresis and when the polarization js 
not proportional to the electric field. The fact 
that at the resonance frequency of a foiled plate 
the impedance minimum is unusually broad in 
Rochelle salt, can be explained without assum- 
ing strong viscous damping. Since an electric 
field changes the elastic compliance sq the 
resonance frequency is modulated by the excit- 
ing a.c. field and the resonance peak is broadened. 
It is sharp only when a large steady field is ap- 
plied across the electrodes. In this case, it follows 
from the data in Fig. 2 that a small variationol 
the field has little effect, and our experiments 
show indeed a considerably sharper resonance 
point when a strong steady field is acting. 
The strongest point in favor of the dipola 
relaxation theory seems to be the fact that at 
high frequencies the dielectric constant in thes 


19R. D. Schulwas-Sorokin and M. V. Posnov, Phys 


Rev. 47, 166 (1935). 
2H. Goedecke, Zeits. f. Physik 94, 574 (1935). 
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direction has considerably smaller values than 
at audiofrequencies. But also this effect can be 
explained by the interaction theory by using the 
discussion in the appendix of paper II. There it 
was shown, for the special case of a plate, 1 cm 
long, in the state of longitudinal vibration in the 
direction at 45° to the } and ¢ axes, that the 
polarization of a volume element at the distance 
y’ from the center of the plate is 


P,=E.[«i— {1—cos (y’w/v)/ 
Xcos (w/ 2v) |. 


The dielectric measurements record the average 
polarization 


and give therefore a susceptibility 


2v 
Ky (1 —— tan w/20) 


For low frequencies tan w/2v=w/2v, and hence 
kw=ki, but for a high frequency, which is not 
near a resonance point of the plate, 2v/w tan 
w/2v=0, whence 


Since = erk®/4a, (Eq. (13), II), where k 
is Mason’s electromechanical coupling coeffi- 
cient, we obtain for the dielectric constant of a 
free crystal at high frequencies 


= €r(1—k?). (19) 


This relation was given by Mason, but his state- 
ment, that ey is the dielectric constant of the 
clamped crystal, is not correct. A high frequency 
field creates a strain distribution in the crystal. 
The average value of these strains is zero, but 
this fact does not imply that ey is identical with 
e.. Above the Curie point ey is larger than e.:, 
but is much smaller than er, the dielectric con- 
stant at low frequencies. Between the Curie 
points Eq. (19) is also valid, but here ey is much 
smaller than ¢,; or er®, however, it is slightly 
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larger than ¢,;°, the dielectric constant of the 
spontaneously deformed clamped crystal. A 
simple calculation leads, for temperatures above 
the Curie point, to the relation 


(xi- 1/x:1) 
2523) +1=1.46. (20) 


For the ferroelectric temperature range an 
analogous relation holds 


(x1°— 1/439) =1.46. (20’) 


From Eqs. (20’) and (16) it follows that the 
susceptibility of Rochelle salt at high frequen- 
cies reaches at the Curie points its largest value 
(ka) max=1.46 Cas/fie. This corresponds to a 
maximum dielectric constant of about 400. At 
the transition temperature near 5°C xy has a 
minimum. Eqs. (20’) and (9c) furnish (xz) min 
= min/(1+}-1.46) which gives, approxi- 
mately, (€7)min=0.6(€r°) min=100. According to 
Eq. (19) this relation implies that k has a 
minimum value of 0.65. Mason measured a 
minimum k=0.60, but since his values of ¢;° are 
rather large, he finds also a large value of 
(€#)min= 160. Direct measurements of the di- 
electric constant at high frequencies have been 
carried out by Mikhailov.'* He checks the theo- 
retical value (€7)min=95, and finds at the upper 
Curie point a maximum of about 250. Also the 
measurements of Busch" and Doborzynski” are 
in fair agreement with the theory. The above 
calculations hold, of course, only if the crystal 
vibrates in the assumed manner. For other modes 
of vibration the dielectric constant at high 
frequencies is different, but in general xq will 
differ only slightly from the dielectric suscep- 
tibility of the clamped crystal (above the Curie 
point) or of the spontaneously deformed clamped 
crystal (between the Curie points). The existence 
of relaxation effects can be accepted only if 
future measurements should lead to still smaller 
values. 
21 G. Busch, Helv. Phys. Acta 6, 315 (1933). 


2 D. Doborzynski, Acad. Polonaise Sci. et lettres, Bull. 
(1938) p. 37. 
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PHYSICAL REVIEW 


The Decay Constant of H® 


R. D. O’Ngeat AND M. GOLDHABER 
Department of Physics, University of Illinois, Urbana, Illinois 
(Received August 2, 1940) 


HE decay constant \ of H? is too small to 

be determined conveniently in the usual 

way from the decay curve.! If, however, we can 

produce a known number WN of Hi? nuclei, and 

observe which number N’ of these nuclei decay 

per unit time, we can obtain \ from the relation 

\=N’/N. We have carried out a determination 

of \ in this manner. H* was produced by exposing 

lithium to slow neutrons, which leads to the 
reaction? Li*-+-n'—He*+ 

(1) In a first experiment lithium metal, sealed 
in a soft glass tube, was exposed inside a large 
paraffin block in a fixed position A to neutrons 
from a 100-mg Ra—a—Be source. The total 
time of irradiation was four months. When the 
irradiation was completed, water was added to 
the Li metal until it was all transformed into 
LiOH whereby hydrogen gas was developed. 
In this way H® can be extracted from the Li 
metal with an efficiency of about 94 percent (see 
below). A known fraction of the hydrogen gas, 
together with some alcohol vapor, was introduced 
into a large Geiger counter, and an activity of 
‘56 counts/min. in excess of a background of 
108 counts/min. was established. This activity 
varied with the pressure in the Geiger counter 
in the expected manner. By ‘“‘scanning’’ the 
counter with a well-collimated beam of y-rays 
we found that it was sensitive over about 90 
percent of its length. If we assume, as is plausible, 
that every H’ nucleus decaying in the sensitive 
portion of the counter is detected, we can calcu- 
late from the above measurements N’ for the 
whole piece of Li metal used. 

N can be put equal to the total number of C 
neutrons absorbed in the Li metal, if we neglect 
two corrections of opposite sign which in this 
experiment should be of the order of a few per- 
cent each: The contribution of neutrons with an 
energy above the Cd absorption limit, and the 
loss of some of the H* nuclei produced near the 


1 L. Alvarez and R. Cornog, Phys. Rev. 58, 197A (1940). 
2 J. Chadwick and M. Goldhaber, Proc. Camb. Phil. Soc. 


31, 612 (1935). 


surface. (Their range in lithium is about 1 /10 
mm.) To obtain the total number of C neutrons 
absorbed in the Li the following experiments were 
carried out. 

_ (a) The ratio of the number of C neutrons 
absorbed in the Li to that absorbed in Mn jp 
position A per unit time was measured with an 
In foil as neutron detector. This foil could be 
wrapped round the samples. 

(b) Mn metal powder was exposed in position 
A, until the 2.6-hr. period produced by sloy 
neutron capture was saturated. A similar exper. 
ment was carried out with the Mn surrounded 
by Cd. The initial activity of the Mn was meas. 
ured in each case under identical geometrical 
conditions by means of a thin-walled Geiger 
counter. 

(c) By mixing a small, known amount of U,0, 
with Mn powder and measuring its activity in 
the same way as was used for the activated Mn, 
the counting arrangement was calibrated. As. 
suming similar 8-ray spectra for Mn® and UX,, 
which we found by absorption experiments to be 
a reasonably good approximation, we obtained 
the total number of 6-rays emitted per minute 
from the saturated Mn sample. This number is 
equal to the number of neutrons absorbed by 
the Mn per minute. 

From these measurements we deduce 


\=7X10-" sec.-'+25 percent, 


or the half-lifetime of H*, 7=31-8 years. 

(2) Through the courtesy of Professor Cork, 
we were able to use for a second set of experi- 
ments Li bombarded, intermittently, during 
three weeks with slow neutrons from the Mich- 
igan cyclotron. As a ‘“‘neutron integrator” which 
allowed us to estimate the total number Of slow 
neutrons absorbed by the Li we used Sb, which 
is activated by slow neutrons with a period of 
60 days, and which was bombarded close together 
with the Li. The hydrogen extracted from this 
Li sample by the procedure described in (1) 
showed an activity, N’, about 100 times larger 
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than in the first experiment. By electrolyzing 
some of the remaining LiOH solution we found 
a small activity in the hydrogen developed. 
From this it followed that the efficiency with 
which H® was extracted from Li by the procedure 
described in (1) was 94+3 percent.’ 

The activity of Sb’* (60-day) was measured 
5 weeks after the last irradiation when Sb’ 
(2.6-day) had already decayed. It was compared 
with the activity of an Sb sample bombarded 
here (position A). If we assume the percentage 
contribution of C neutrons to the 60-day activity 
to be the same under the different conditions of 
bombardment employed in Michigan and here, 
we can estimate the number of H# nuclei formed 
in the Li sample. (This assumption, however, 
may introduce a serious error, as we found by a 
separate experiment (position A) that the con- 
tribution of neutrons above the Cd absorption 


*This high efficiency is not unexpected as a era 
high efficiency has been found for the production of H 
in the reaction LiH+D,0—HD+LiOD (Beutler, Brauer 
and Juenger, Naturwiss. 24, 347 (1936)). 


limit to the 60-day activity of Sb is considerable.) 
In this way we obtained a value of about 50 
years for the half-lifetime of H*. Though this 
experiment is less accurate than the first it 
yielded a valuable check on our method. 

We had observed previously‘ that the electrons 
emitted from H* have a range of about 0.5 mg 
Al/cm?*. This corresponds to an upper limit of 
the B-ray spectrum Eo>=15+3 kev. As the 
transition may be considered as allowed, we can 
now check the relation between Ep, and X, as 
given by Fermi’s theory of §-decay, for this 
rather extreme case. Using the value for the 
constant of proportionality as deduced by Grén- 
blom® from an analysis of the 8-ray data on light 
elements, we find good agreement between theory 
and experiment within the present accuracy. 

This research was aided by a grant from the 
Graduate School Research Fund of the Univer- 


sity of Illinois. 


' ous} D. O'Neal and M. Goldhaber, Phys. Rev. 57, 1086A 
+ B. O. Grénblom, Phys. Rev. 56, 508 (1939). 
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Connection Between the Meson Decay and the 
Beta-Decay 

According to Yukawa’s theory of the meson,! the beta- 

disintegration process can be described as the transforma- 

tion of a neutron (JN) in the nucleus into a proton (P) with 

the emission of a meson (Y~) and the disappearance of the 

meson with the simultaneous emission of an electron (e~) 

and an antineutrino (m’). Thus the two basic hypotheses 

of his theory introduce the following interactions between 
elementary particles: 

NeP+Y-, P2N+Y!*, (1) 

Yteet+n. (2) 


From the interaction of the latter type, connecting a meson 
with light particles, it follows immediately that the mesons 
are unstable and disintegrate spontaneously each into an 
electron and a neutrino. The proper lifetime, i.e., the mean 
lifetime at rest, of the meson calculated from the above 
theory was given by? 

to~(g*/he) -(m/m,)>-T, (3) 


where g*/hc is a dimensionless constant which characterizes 
the strength of the interaction (1) and which can be de- 
termined by the magnitude of the nuclear forces. In Eq. 
(3) m, and m denote masses of the meson and the electron, 
respectively; T is a time constant which appears in the 
theoretical formula for the beta-decay and which can be 
estimated by comparing the formula with lifetimes of light 
nuclei.* Since the most probable values for these constants 


are 
my~170m, T~105 sec., (4) 


we find from (3) 
sec. (5) 
On the other hand, a striking confirmation of the meson 
decay has been received from observations on the mass 
absorption anomaly for cosmic rays‘ and more direct 
evidence has also been reported by Williams and his 
collaborators.’ These observations all indicate that the 
proper lifetime of the meson is about 


ro~1 to 5X 10-6 sec. (6) 


Thus the theoretical estimate (5) is too short by a factor 
10 in comparison with the experimental value. 

We shall point out here that there are at least two ways 
of removing this discrepancy: 

(I) As was recently suggested by Mgller, Rosenfeld and 
Rozenthal,* we assume that there exist two kinds of mesons 
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and that while one of them has a lifetime of the order of 
10~* second, in accord with experiments, the lifetime of the 
other is so short that it can explain the sufficiently rapiq 


nuclear beta-decay. 
(II) We reverse the connection between the meson 


and the beta-decay. Namely, instead of (2) we sha 
introduce the original Fermi interaction 


() 
as a new basic hypothesis. If this is done, the meson decay 
takes place as a compound process of (1) and (7) in the 
following way. First a negative meson, for example, is 
absorbed by a proton, which is in the negative energy state 
and a pair of a neutron and an antiproton (P’) is created 
Then this pair disappears with the simultaneous emissigg 
of an electron and an antineutrino: 

Y-—>P’+ N-e~+n’. (8) 


From this theory, we obtain for the proper lifetime of the 
meson 


ro~(he/g*) -(m/m,)§-T, 


which is about (hc/g*)* times as long as (3) and whichis 
in good accord with the experimental value (6). 

It is to be noted that in this calculation there appears 
a diverging integral with respect to the momentum of the 
heavy particle which is virtually created in the inter. 
mediate state. This difficulty is quite similar to those which 
arise in calculations for the self-energy and the magnet 
moment of the heavy particle.? Hence we have cut off this 
integral at such a value of the momentum that it gives th 
right order of magnitude for the self-energy and th 
magnetic moment. 

A more detailed account will be published later in th 
Proceedings of the Physico-Mathematical Society of Japas 

In conclusion, I wish to express my cordial thanks t 


Professor H. Yukawa for his kind interest in this work. 
SHOICHI SAKATA 


Institute of Theoretical Physics, 
Kyoto Imperial University, 
Kyoto, Japan, 
July 30, 1940. 
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‘Mallen Rosenfeld and Rozenthal, Nature 144, 629 (1939). 
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The Nuclear Spin of C¥ 


Using Rabi’s magnetic-resonance molecular-beam 
method, Hay has observed both the sign and the magnitude 
of the gyromagnetic ratio of the C® nucleus.! By making 
the rather safe but still slightly theoretical assumptions 
that he is dealing with a *P state and that the spin part S 
of the “nuclear spin” J consists of neutron spin with no 
large admixture of proton spin, he is able to conclude from 
the sign of the gyromagnetic ratio that the nuclear spin 
js 1/2. This disagrees with the nuclear spin 3/2 inferred 
earlier from analysis of a band spectrum which could not 
be completely resolved? It is unfortunate, especially in 
view of such a disagreement, that an experimental de- 
termination should have to lean so heavily on theoretical 

eses. 

3" possible to reach Hay’s conclusion from a well- 
established empirical rule, instead of theoretical hypothesis, 
if we make use of the observed magnitude of the gyro- 
magnetic ratio as well as its sign. The rule is that all 
hitherto known magnetic moments of odd nuclei (there are 
more than fifty of them) lie between the limits given by the 
simple single-particle model.* The corresponding limits for 
the gyromagnetic ratio g in the odd-neutron case are 


—p,/(I+1), 


where uy = — 1.96, the lower limit of the observed magnetic 
moment of the neutron,‘ in nuclear magnetons. The 
measured gyromagnetic ratio, g= +1.401+0.004, does not 
lie strictly within these limits for any value of J. However, 
it exceeds the upper limit 1.30 for J=1/2 by only a few 
percent. For higher values of J, it exceeds the limit by 
about a factor two or more. 

In Fig. 1, nuclear magnetic moments are plotted against 


ee 

ee bd 

“t+ 3 


Fic. 1. The magnetic moments and nuclear spins of the odd-neutron 
nuclei. The triangles are possibilities for C™. 


nuclear spins, for odd-neutron nuclei only.? The upper 
line is the limit given by the single-particle model with 
I=L—1/2, and the lower line with J=L+1/2. The dots 
represent the previously known data for the odd-neutron 
nuclei. (The data for the odd-proton nuclei, not shown, are 
about twice as numerous.) The triangles show two possi- 
bilities for C'*. The previous data for small values of J tend 
to crowd close to the limit given by J=L—1/2, both in 
Fig. 1 and in the odd-proton case,? and the triangle for 
I=1/2, just above the limit, does not deviate appreciably 
from the regularity of the other data. The triangle for 
I=3/2 would be seriously discordant with the other data, 
and may therefore be rejected. 

One notes that the nuclear spin 1/2 for C¥, thus de- 
termined in another way from Hay’s data with a minimum 
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of theoretical description, is the value predicted both by 
the Hartree model® and by the alpha-model® on the basis 
of the Larmor-Thomas type of spin-orbit coupling. The 
magnetic moment is smaller than these models would 
suggest. 


D. R. INGLIs 
Johns Hopkins University, 
Baltimore, Maryland, 
August 16, 1940. 


1R. H. Hay, Phys. Rev. 58, 180 (1940). 
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3 Th. Schmidt, Zeits. f. Physik 106, 358 (1937). 

4L. W. Alvarez and F. Bloch, Phys. Rev. 57, 111 (1940). 

5M. E. Rose and H. A. Bethe, Phys. Rev. 51, 205 (1937). 
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Concentration of Radioactive Bromine Produced by 
Gamma-Radiation 


One of the difficulties encountered in the study of nuclear 
reactions of the (y, ) type is the small beta-ray activity 
ordinarily obtained, even with the strongest gamma-ray 
sources. The high energy gamma-rays used for such work 
(vig. those from the reaction Li’(p, y)Be*) are but weakly 
absorbed in matter, and this suggests that a concentration 
of the radioactive atoms produced in a large bulk of 
material should increase the activity obtained. Moreover, 
the (7, m) reaction, which is the only one thus far observed 
under gamma-ray bombardment, leads to an isotope of the 
irradiated element. Accordingly, we have made some ex- 
periments with a view to adapting the Szilard-Chalmers' 
technique to this problem. 

Bromine was irradiated in the form of ethyl bromide, 
500 cc of which were contained in a spherical vessel con- 
structed so as to surround almost completely the source of 
gamma-rays. The active bromine was separated as silver 
bromide in the usual way. A composite decay curve was 
obtained which was resolvable into periods of 6.4 minutes, 
18 minutes and 4.5 hours with relative intensities at 
saturation of 1 : 3 : 2, respectively. The 6.4-minute period 
is that of Br’ and the others belong to the isomeric Br™. 
These isotopes are formed from the two stable ones Br7* and 
Br*!, in each case by the loss of a neutron. It is noteworthy 
that the method of concentrating the activity also estab- 
lishes the chemical nature of the isotope. 

The radioactive sources prepared in this way were at 
least ten times as strong as those obtainable without con- 
centrating the activity. Beta-ray sources with initial ac- 
tivity of several thousand counts per minute were made. 
Contaminating activity produced by fast neutrons from 
the normal deuterium content of the ion beam striking the 
lithium target was shown to be negligible. 

WALTER H. BarxKas* 
R. CARLSON 
J. E. HENDERSON 
H. Moore 
PGniversity of Washington, 
Seattle, Washington, 
August 28, 1940. 


* Member, The Institute for Advanced Study, Princeton, New Jersey. 
1L. Szilard and T. A. Chalmers, Nature 134, 494 (1934). 
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Luminescence Excited by Exposure to Neutrons* 


Since many elements, during exposure to neutrons, re- 
lease energy in the form of gamma-rays, beta-rays and 
alpha-particles, it is reasonable to expect that these radia- 
tions may be used to excite luminescence. We have made 
- observations which show that certain luminescent ma- 
terials, which are known to be excited by x-rays, cathode 
rays, and radiations from radium, emit light as a result of 
exposure to neutrons. We have used photographic methods 
to detect this emission during exposure, and have made 
visual observations, after exposure, of thermoluminescence 
excited in many substances and of phosphorescence in 
luminescent ZnS. 

Exposures were made for varying lengths of time up to 
several hours at distances up to 50 cm from a neutron 
source with radium equivalent 100 kg and at a distance of 
a few cm from a source with radium equivalent about 1 kg. 
Specimens were simultaneously exposed in black paper, 
inside a solid box of paraffin 12 cm thick, and in a lead 
box 1 cm thick placed inside this same paraffin box. The 
effect was observed to be least in the specimens wrapped 
in black paper and exposed to fast neutrons. The two sets 
of specimens inside the paraffin box were observed to emit 
luminescence of practically the same intensity although 
one set was enclosed in a lead box and the other was not. 

Photographic observations.—Two identical specimens of 
luminescent CaWO,, one wrapped in black paper and the 
other in Cellophane, were placed in contact with the sensi- 
tive side of a panchromatic film and exposed to neutrons 
as described. The part of the film under the Cellophane 
was strongly blackened while the part under the black 
paper showed little or no blackening. A similar difference 
was observed when the CaWO, was painted on black card- 
board with a binder and part of it covered with black paper 
during exposure. Other materials were observed to show 
this same effect especially CdWO, activated by 1 percent 
Mn and phosphorescent ZnS, It is of interest to note that 
the three substances which we have found to be most 
strongly luminescent under neutron exposure are all known 
to show scintillations under alpha-particles from radium. 

Thermoluminescence.—After exposure to neutrons, about 
thirty substances were found to emit light when they were 
heated on a small electric stove to a temperature consider- 
ably below that of red heat. The strongest effect was shown 
in fluorites from certain localities, kunzite, ThO2 activated 
by Tb, and in CaSO,, CaF2, and CdSO, activated by Mn. 
The thermoluminescence is, in most cases, very dim and 
not comparable in intensity with that excited by x-rays. 
The intensity may be increased by mixing with the lumi- 
nescent material some substance which reacts strongly with 
neutrons such as GdSQ,. The effect is diminished by screen- 
ing the luminescent substance during exposure with boron 
carbide. The effect of surrounding some specimens such as 
fluorite and CaSO, excited by Mn, with cadmium metal 
during exposure is to increase the intensity of thermo- 
luminescence excited because, although cadmium absorbs 
slow neutrons, it emits, as a result of this absorption, 
radiations which are effective in excitation. 

Summary of results —(1) Luminescence is excited in some 
materials as a result of exposure to neutrons. In a material 
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for neutron excitation it is desirable that the lumi 
material contain some element which is highly reactive to 
neutrons and that it contain some heavy element capable 
of absorbing radiations such as gamma-rays, Materials 
excited by neutrons are excited also by x-rays, 
rays or the radiations from radium. (2) Slow neutrons are 
more effective as a source of excitation than fast neutrons, 
(3) The intensity of luminescence may be increased by 
mixing with the luminescent material a small quantity of 
some material containing an element which js 
reactive to neutrons. (4) Excitation of luminescence 
exposure to neutrons may be explained as follows: Neu. 
trons react with the nuclei of certain atoms which form a 
component part of the luminescent substance or with the 
nuclei of atoms of some other substance highly reactive to 
neutrons which may be mixed with the luminescent sub. 
stance or placed near it. As a result of this reaction, radia. 
tions are emitted which excite luminescence. Gamma-tays 
appear to be a factor in the excitation and secondary radia. 
tions and internal conversions may also play a part, 

The exposures to neutrons were made by the use of the 
cyclotrons at Columbia and Cornell Universities through 
the kindness of Dr. J. R. Dunning and Dr. R. F. Bacher 
whose assistance we greatly appreciate. 

FRANceEs G. Wick 
MABEL S. Vincent 


Department of Physics, 
Vassar College, 
Poughkeepsie, New York, 
August 15, 1940. 


* This research was made possible by a grant from the American 
Society. 
1 J. G. Hoffman and R. F. Bacher, Phys. Rev. 34, 644 (1938), 


On the Isotopic Weights of Chlorine, Argon and Iron 
by the Doublet Method 


With a mass spectrograph of Bainbridge-Jordan type 
whose resolving power is about 17,000, the isotopic weights 
of Cl*, Ci87, A*® and Fe®* were measured by the doublet 


method. 


The cylindrical discharge tube employed in this work’ 


was provided with an aluminum bar anode and an iron 
concave cathode. All photographic plates used were of the 
Schumann type prepared in our laboratory. 

The plate was calibrated with the bromine isotopes Br® 
and Br*, for the mass difference of which we accepted 2500 
packing fraction units as determined by Aston. 

(1) Chlorine.—Three doublets 
and C,H.—HCI*’ were obtained by electric discharge 
through the vapor mixture of normal hexane and carbon 
tetrachloride. 

These mass differences were as follows: 


C®,—HCl* =246.7+1.7, (AM.10~) 
C®,H —Cl** =421.7 
=419.8-41.1,f 
From these results and from the values H'=1,00812 
+0.00004 ! and C= 12.00394+0,00018,? the following is 
topic weights of Cl* and Cl*’ were obtained: 
Cl* = 34,97903 +0.00038, 
Cl? = 36.97786 +0.00036. 
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(2) Argon.—By introducing the mixture of normal 
hexane vapor and argon into the discharge tube, the two 
doublets C®,H,—A*® and C®,;H;—HA* were photo- 
graphed. The following mass differences of these doublets 


obtained, 
(1) 


HA“ =693.0+2.3. (2) 


From the results (1) and (2), the isotopic weight of 
A* is as follows: 


(1) A*®?=39.97637+0.00057, 
(2) A*®=39.97500+0.00662. 


At present, we have no interpretation of the discrepancy 
between the mass differences of the above two doublets. 
(3) Iron.—The wide and faint doublet C,H; — Fe® was 
photographed, from which we obtained 1235+17 as its 
mass difference. The ion of Fe** was considered to originate 
from the iron cathode. We obtain a provisional value 
55.9572+0.0012 as the isotopic weight of Fe**. Its packing 
fraction —7.7 is smaller than Dempster’s value —7.0+0.4.* 
T. Oxupa 
KK. Ocata 
K. AOKI 
Y. SUGAWARA 
Laboratory, 
Imperial University, 
Osaka, Japan, 
July 22,1940. 
1f, W. Aston, Proc. Roy. Soc. A163, 391 (1937) 


2 Asada, Okuda, Ogata and Yoshimoto, Nature 143, 797 (1939). 
2A, J. Dempster, Phys. Rev. 53, 64 (1938). 


Electrostatic Generator with Concentric Electrodes 


Our 2.6-million-volt generator has been described in 
previous issues of this journal.» * It was housed in a tank 
20 ft. long and 54 ft. in diameter, with a working pressure 
of 100 Ib. per sq. in. With only high pressure air in the tank 
this generator was limited to approximately 2.2 Mv by 
sparking directly across the air gap. As Freon was added 
to the high pressure air the direct spark-over voltage rose, 
but at 2.6 Mv sparking set in along the charging belt and 
additional Freon then gave no improvement. 
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This apparatus was dismantled, the tank was provided 
with a large, removable end plate, and new apparatus was 
installed as shown in Fig. 1. In the new installation the 
apparatus is supported from one end only by three Texto- 
lite tubes, T, of 3? inches outside diameter and } inch wall 
thickness. This arrangement gives a charging belt length 
and a minimum insulator length (except for insulator #) 
almost twice as great as in the previous installation. The 
new tube consists of 62 sections, each 2} inches long, 
compared to 53, 24-inch sections for the old tube. 

Electrodes B and C were expected to give improvement 
in the spark-over voltage for two reasons: First, they pro- 
vide a more uniform gradient between electrode A and the 
tank wall; and second, test work indicated that high pres- 
sure air withstands higher gradients in a short gap than in 
a long gap. 

Adjustable corona gaps provide potential distribution 
along the supporting structure and determine the poten- 
tials of electrodes B and C relative to the potential of A. 
Insulator K supports the hemispherical end of electrode C. 
Flexing of the structure by electrostatic forces would 
probably cause trouble if this electrode were supported 
only from its open end. 

The potential of electrode C is measured by means of a 
generating voltmeter. This measurement also gives a rough 
indication of the total voltage. In a recent test run the ion 
beam was brought out into the air through an aluminum 
foil with an air equivalent of approximately 12 mm and the 
range of the proton beam, measured visually, was used to 
determine the generator voltage. 

With an air pressure of 100 Ib./in*. in the tank, the 
maximum range of the proton beam was 18 cm (room tem- 
perature 27.2°C, barometric pressure 742 mm of Hg). After 
adding the air equivalent of the aluminum foil this gave a 
voltage of 3.5 Mv from the 1937 range curve of Bethe and 
Livingston. When 10 Ib. (weight) of Freon was added, the 
maximum range of the proton beam was 27 cm, giving a 
voltage of 4.3 Mv. The voltage was limited by direct spark- 
over, although insulator K gave some trouble. With 15 Ib. 
of Freon the maximum range was 29 cm (4.5 Mv). The volt- 
age was then limited by sparking along insulator K. In 
previous tests it withstood higher voltage, but in recent 
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Fic. 1. A, B, C; aluminum shells } inch thick. T; Textolite tubes 3} inch outside diameter and } inch wall thickness. Each Textolite tube is 


equipped with aluminum rings, made up from }-inch diameter rod, clamped 
cap H; aluminum hoops supported by studs with saddle-shaped ends. One stud projects from each ring. 


nches apart (not shown in the 


tightly to the Textolite and spaced 1 a past 
extolite tube equipped 


i 
s and a corona gap system for potential distribution. Porcelain rings for the accelerating tube are 2} inches long, corrugated inside 
make the seals. Springs 


and outside. The separating electrodes are steel spinnings. 


Gaskets of rubber dam 0.01 inch thick 
force of 5000 Ib. on the tube. Belt; cotton fabric, woven endless, 16 inches wide 


apply a compressional 
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exposure to the atmosphere during hot, humid weather it 
soaked up too much moisture and drying treatments at- 
tempted so far have not been entirely successful. 

The performance of the tube and charging belt indicate 
that they will be satisfactory for still higher voltage. 

We are indebted to Dr. J. L. McKibben for many 
valuable suggestions and for much help in the design work. 
Results of his test work were utilized and contributed 
greatly toward the success of the generator. N. D. Crane 
and A. O. Hanson gave valuable help with construction 
work and the Wisconsin Alumni Research Foundation 
provided financial support. 

R. G. HERB 

C. M. TURNER 

C. M. Hupson 

R. E. WARREN 
siversity of W leconsin, 


Madison, Wisconsin, 
August 15, 1940. 


Pm, Roa Herb, D. B. Parkinson and D. W. Kerst, Phys. Rev. 51, 
AN B. Parkinson, R. G. Herb, E. J. Bernet and J. L. McKibben, 
Phys. Rev. 53, 642 (1938). 


Internally Converted Gamma-Rays from Radioactive Gold 


The electron spectrum of radioactive gold chemically 
separated from platinum bombarded with 9.5-Mev deu- 
terons has been obtained with the magnetic spectrometer. 
A typical result with various peaks due to internally con- 
verted gamma-rays is shown in Fig. 1. It is evident that 
the spectrum of this element is very complex. The lower 
curve in this figure was taken 18 days after the upper one, 
and was obtained in order to evaluate the half-life of each 
electron group. 


Fic. 1. Electron spectra of radioactive gold[taken 18 days rt. 
The ordinates of the lower curve have by beter 
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The results so far obtained may be summarized as fol. 
lows: The activity of 164 days half-life shown! by a 
tion measurements to emit beta-particles of energy 0,45 
Mev and gamma-radiation of energy 0.11 Mev is too 
for analysis in the spectrometer. The activity of half-life 
5.6 days, which from absorption measurements was re 
ported to emit beta-particles of energy 0.36 Mey and 
gamma-radiation of energy 0.41 Mev, is now shown to be 
due to a partially converted gamma-ray of energy 35644 
kev. The responsible isotope may be either Au! or Ayis 
The other 4 peaks of Hp greater than 1800 oersted cm pa. 
be resolved into K- and L-conversion electron groups coy. 
responding to gamma-ray energies of 331+3 and 41044 
kev, respectively. The half-life of these gamma-rays ap- 
pears to be about 3 days which is somewhat longer than 
the value reported from an analysis of the decay curyes 
taken with an ionization chamber. Early observations jn 
the spectrometer confirm the existence of a positron emitter 
of short half-life. Further work with stronger samples must 
be carried out to analyze satisfactorily the electron groups 
of lower energy. 

This investigation was made possible by a grant from 
the Horace H. Rackham Fund. 

J. L. Lawsoy 


J. M. Corx 
Department of Physics, 
University of Michigan, 
Ann Arbor, Michigan, 
August 16, 1940. 


1J. M. Cork and J. Halpern, Phys. Rev. 58, 201 (1940). 


The Photographic Registration of Heavy Particles 
Emitted During Bombardment 


In order to observe the particles emitted from a target 
during its actual bombardment in the cyclotron, inde. 
pendent of the disturbing effects of the scattered primary 
particles and the neutron background, a magnetic spectro- 
graph has been devised. The arrangement employed is 
shown in Fig. 1. 

The primary beam of deuterons is incident upon a narrow 
strip of very thin foil of the element being studied, sus- 
pended at the center of a cavity in a cast lead block. 
Particles leaving the target at an angle of ninety degrees 
are collimated by a slit in the wall of the lead block. The 
strong field of the cyclotron together with an auxiliary field 
produced by a magnetic shunt on the gap of the large 
magnet, suffices to resolve the emergent particles. Their 
deflections are approximately proportional to ne/(mE) 
where ne is the charge, m the mass and E the energy of the 
particle. The photographic plate is held in a light-tight 
carriage that can be set at any desired position on an 
accurately milled track. By placing immediately in contact 
with the plate during bombardment, a stepped aluminum 
foil it is possible to also observe the range of each group of 
particles. This gives sufficient information to determine the 
nature of the particle and its energy. 

Preliminary observations have been made on aluminum, 
copper, palladium, gold and platinum, using foils about 
0.00004 in. thick. Considering the geometry of the sample 
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Fic. 1. Schematic diagram of magnetic spectrograph. 


and slit the lines are of satisfactory width. Platinum and 
gold show mainly only the scattered deuteron beam. In 
aluminum only protons are observed. The maximum range 
of the protons is 126 cm air equivalent or about 11.0 Mev 
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which is in agreement with the early observations of Mc- 
Millan and Lawrence,' considering their bombarding energy 
was 2.2 Mev. Palladium shows in addition to scattered 


_ deuterons, protons of low energy with something like a 


continuous distribution. With copper a proton group of 
energy about 8.5 Mev is observed together with the 
scattered deuterons. An additional line in the case of copper 
as yet not satisfactorily explained has a range of about 47 
cm air equivalent but shows a deflection even smaller than 
would be expected for alpha-particles of this energy 
(25 Mev). 

The apparatus with certain improvements should prove 
to be a valuable aid in the study of nuclear reactions. 

This investigation was made possible by a grant from 
the Horace H. Rackham Fund. 

J. M. Cork 


P. F. BARTUNEK 
Department of Physics, 
University of Michigan, 
Ann Arbor, Michigan, 
August 15, 1940. 


1E. McMillan and E. O. Lawrence, Phys. Rev. 47, 343 (1935). 
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